Knot Invariants from Four-Dimensional Gauge Theory 



Davide Gaiotto^ and Edward Witten^'^ 

^School of Natural Sciences, Institute for Advanced Study, 
1 Einstein Drive, Princeton, NJ 08540 USA 



^Department of Physics, Stanford University, Palo Alto CA 9^305 

Abstract: It has been argued based on electric-magnetic duality and other ingredients 
that the Jones polynomial of a knot in three dimensions can be computed by counting the 
solutions of certain gauge theory equations in four dimensions. Here, we attempt to verify 
this directly by analyzing the equations and counting their solutions, without reference to 
any quantum dualities. After suitably perturbing the equations to make their behavior 
more generic, we are able to get a fairly clear understanding of how the Jones polynomial 
emerges. The main ingredient in the argument is a link between the four-dimensional gauge 
theory equations in question and conformal blocks for degenerate representations of the 
Virasoro algebra in two dimensions. Along the way we get a better understanding of how 
our subject is related to a variety of new and old topics in mathematical physics, ranging 
from the Bethe ansatz for the Gaudin spin chain to the M-theory description of BPS 
monopoles and the relation between Chern-Simons gauge theory and Virasoro conformal 
blocks. 
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1 Introduction 

The Jones polynomial [1] is an invariant of knots that has multiple relations to many 
aspects of mathematical physics, including integrable lattice statistical mechanics, two- 
dimensional conformal field theory and associated representations of braid groups, and 
three-dimensional Chern-Simons gauge theory. Khovanov homology [2] is a more recent 
topological theory in four dimensions; in this theory, a knot is viewed as an object in three- 
dimensional space and the invariant associated to a knot is a vector space (of physical states) 
rather than a number. The relation between the two theories is that the four-dimensional 
theory associated to Khovanov homology, when compactified on a circle, reduces to the 
three-dimensional theory that yields the Jones polynomial. 

Khovanov homology has been interpreted physically [3] in terms of topological strings, 
building on earlier work on BPS states of open strings [4]. See [5]- [13] for a sampling 
of additional developments. An alternative but closely related physical interpretation of 
Khovanov homology has been given in [14], where more detailed references can be found 
concerning the Jones polynomial, Khovanov homology, and their relations to mathematical 
physics. 

According to this more recent proposal, the Jones polynomial can be computed by 
counting the solutions of certain elliptic partial differential equations in 4 dimensions, and 
Khovanov homology can then be constructed by counting the solutions of related equations 
in 4-t-l dimensions. The reasoning that led to this proposal relied on electric-magnetic 
duality of = 4 super Yang-Mills theory in four dimensions to transform one description 
that is rather "quantum" in nature (being closely related to Chern-Simons gauge theory 
on a bounding three-manifold) to another that is "semiclassical" in the sense that the 
partition function can be computed just by suitably counting the classical solutions of 
certain differential equations. 

Instead of relying on electric-magnetic duality to predict this perhaps mysterious result, 
can we understand it by a direct study of the equations? This is the goal of the present 
paper. We will gain a reasonable degree of understanding of the Jones polynomial and a 
good foundation for understanding Khovanov homology. 
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Figure 1. A knot has been placed at the boundary of the four-manifold M4 = W x 
1.1 A Brief Review 

The four-dimensional equations in question can be described as follows. The gauge group 
is a compact Lie group^ G. The fields in the equations are a gauge field A which is a 
connection on a G-bundle E — t- M4, with M4 an oriented Riemannian four-manifold, and 
another field (p that is a one-form valued in the adjoint representation of G. The equations, 
which were first studied in relation to the geometric Langlands correspondence [15], read 

{F - (p A(l) + tdA(l))^ = 
{F - (p A (p-t-'^d A(l)y = 

dA*<P = 0, (1.1) 

where the selfdual and anti-selfdual projections of a two- form b are denoted 6^; d^ = 
d -|- [A, ■ ] is the gauge-covariant exterior derivative; F = dA + A A A is the Yang-Mills 
field strength; ★ is the Hodge star operator; and t is a real parameter. (Actually, t takes 
values in MP^ = R U 00; for t ^ or t — 7- 00, one multiples the second equation by t or 
the first by t~^.) To study knot invariants, one specializes to M4 = W x M_|_, where W is 
a three-manifold and 1R+ is the half-line y > (fig. 1). 

The boundary condition at y = is slightly subtle but can be easily described in the 
absence of knots. Suppose first that t = 1. Since the boundary condition is local, we can 
specialize to Vl^ = in describing it. (In any case, that is the main example for the 
present paper.) Consider a classical solution that is invariant under translations along 
and such that A and the part of (f) normal to the boundary vanish. The equations then 
reduce to Nahm's equations for (p, the part of (p tangent to the boundary: 

^ + 0X(^=O. (1.2) 
dy 

These equations have a singular solution, first introduced by Nahm in his work on monopoles. 
Pick an embedding p : su(2) — ?• g (where su(2) and g are the Lie algebras of SU{2) and G, 

^We reverse notation from [14], writing for tiie gauge group in the Chern-Simons description and G 
for the gauge group in the dual "magnetic" description, on which we focus in this paper. 
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respectively), given by a triple of elements t £ q obeying [ti,t2] = *3, and cyclic permuta- 
tions. Then the solution is 

= -. (1.3) 

y 

Though any p gives a solution, the case we want is that p is a principal embedding. For 
G = SU{N), this means that p is an irreducible embedding of SU{2) in G; for any G, it 
means that the raising operator = ti + it2 is a "regular" element of the complexified Lie 
algebra qc (this means that the subalgebra of qc that commutes with has the minimum 
possible dimension) . Then one can define a boundary condition by allowing precisely those 
solutions of (1.1) that can be approximated for y — t- by the model solution (1.3) with the 
regular Nahm pole. This has an analog for any t; the starting point, as explained in an 
appendix, is to set the tangential part A of the gauge field to be a specific multiple of cp, 
so that the equations reduce again to Nahm's equations. 

When a link L G W is included, this boundary condition is modified along L. A link 
is simply the union L = UiKi of disjoint embedded circles Ki. The Ki are labeled by 
representations R'^ of the Langlands or GNO dual group to G, and in this description 
the knots enter the formalism only via the way they enter the boundary conditions. Roughly 
speaking, the modification is made by requiring the presence of singular BPS monopoles 
supported along the Ki with magnetic charges given by the i?^. 

The G-bundle E — )• M4 has an instanton number P defined in the usual way as a 
multiple of J^^^ Tr F A F. (The definition of P as a topological invariant involves some 
subtleties that are described in [14]; roughly speaking, the boundary conditions at the 
finite and infinite ends of M+ give suitable trivializations of E, enabling one to define the 
instanton number.^) For each value n of the instanton number, one defines an integer a„ 
by "counting" (with signs that are determined by the sign of the fermion determinant of 
A/" = 4 super Yang-Mills theory) the number of solutions of the super symmetric equations 
(1.1) with instanton number n. Then the partition function of a certain version of twisted 
A/" = 4 super Yang-Mills theory on M4 is 

Z(g) = ^a„(?-, (1.4) 

n 

where the definition of q in terms of parameters of A/" = 4 super Yang-Mills theory was 
explained in [14]. 

To get the Jones polynomial and its analogs for other groups and representations, one 
specializes to = and takes the boundary condition at y = 00 to be simply A,(l)^0. 
Then for example for = SU{2) and the two-dimensional representation of SU{2), 
Z{q) is supposed to become the Jones polynomial. Since Vl^ = is the case relevant to 
the Jones polynomial, it will be the main example in the present paper. However, many 
of our considerations apply also for W = M x C where C is a Riemann surface, so we will 
consider this case as well. 

■^For general W and a general choice of the boundary condition at y = 00, P takes values not in Z but 
in a certain coset of Z in Q. 
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A slight generalization of the above-described procedure is to modify the boundary 
condition at infinity so that A and (py vanish but 4> approaches, up to a gauge transfor- 
mation, a specified triple a of elements of t, the Lie algebra of a maximal torus T of G. 
Physically, this means that one takes the vacuum at infinity to be specified by a given 
point on the Coulomb branch. (In the presence of the Nahm pole boundary condition, 
turning on (py or the other two scalars of = 4 super Yang-Mills theory - called a, a in 
[14] - would break supersymmetry; so a are the only useful Coulomb branch parameters.) 
Continuously turning on Coulomb branch parameters should not affect the counting of 
solutions of an elliptic equation, so this procedure should give a slightly more general way 
to compute the Jones polynomial. To describe the basic solution of the equations (1.1) 
with a specified choice of a at infinity, one looks for a solution that still has A = (py = 
and is still invariant under translations along R^, but now obeys lim^^oo = g<ig~^, for 
some g £ G. The equations still reduce to Nahm's equations (1.3). A general theorem 
[16] says that for any simple Lie group G, and any specified choice of a, there is a unique 
solution of Nahm's equations with a regular Nahm pole at y = and the required behavior 
for y — 7- oo. This solution describes the ground state at the given point on the Coulomb 
branch in the absence of any 't Hooft operators on the boundary. 

1.1.1 Lift To Khovanov Homology 

Though our main focus will be to recover the Jones polynomial from this framework, we 
will also briefiy sketch how Khovanov homology is supposed to arise. A primary purpose 
of this is to explain the extent to which the particular values t = ±1 are or are not special, 
since this will be important later. 

To get Khovanov homology instead of the Jones polynomial, we are supposed to "cat- 
egorify" the above-described situation, which is just a fancy way to say that we must 
obtain everything that has been described so far from a theory in one dimension higher. 
For this, let x^,x'^,x^^ be local coordinates on W and decompose (p as (f) = (j) ■ dx + (pydy, 
where (py is the component of (p in the y direction. Categorification is accomplished by 
introducing a new time coordinate x^ and replacing (py by the covariant derivative D/Dx^. 
This replacement makes sense in that, since (py only appears in (1.1) inside commutators 
and covariant derivatives, the replacement does give a differential equation (rather than a 
differential operator), now on the five-manifold M x W x ]R_|_. Moreover this differential 
equation, whose details are described in section 5 of [14], is elliptic so problems of counting 
its solutions make sense. ^ 

This five-dimensional lift of the four-dimensional equations (1.1) also has a surpris- 
ing four-dimensional symmetry, provided we set t = ±1. The original four-dimensional 
symmetry relating the different directions in M4 = W x M+ has been spoiled by the re- 
placement (py — )• D/Dx^. But at t = ±1, the five-dimensional equations acquire a new 
four-dimensional symmetry: one can replace M x W hy a general oriented Riemannian 
four-manifold M, without additional structure, and formulate these equations on M x M_|_. 
For studying the Jones polynomial, the values t = ±1 are not particularly distinguished; 

^This equation has also been formulated and some basic properties described in [17]. 
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the counting of solutions of the elhptic equations (1.1) is independent of t. Moreover, cat- 
egorification - the substitution (l)y — D/Dx^ - is not hmited to t = ±1. What is special 
about t = ±1 is the four-dimensional symmetry of the categorified theory, which is likely 
to have important implications for Khovanov homology and its analogs on other manifolds. 

From a physical point of view, the five-dimensional lift of the equations (1.1) are BPS 
conditions of a certain twisted version of five-dimensional super Yang- Mills theory, formu- 
lated on M. X W X M_|_; they describe configurations that are invariant under one of the 
supercharges, which we will call Q. This operator obeys = 0, and the space of super- 
symmetric ground states is the same as the cohomology of Q. This is the candidate for 
Khovanov homology. Mathematically, the five-dimensional equations can be interpreted as 
Morse theory flow equations, and the space of supersymmetric ground states is the analog 
of Floer homology for this situation. Physically, to construct the space of supersymmetric 
ground states, one starts with time-independent solutions of the five-dimensional equa- 
tions - these are simply the solutions of the original uncategorified equations (1.1) in four 
dimensions. Expanding around any one of these solutions, one can construct an approx- 
imate supersymmetric state, and these furnish a basis for the space of supersymmetric 
states in the classical approximation. Then one computes quantum corrections by taking 
account of tunneling between classical vacua; the tunneling events are solutions of the full 
five-dimensional equations. 

From this point of view, the link between Khovanov homology and the Jones poly- 
nomial comes from the fact that the classical solutions that give a basis for the classical 
approximation to Khovanov homology are the same ones that must be counted to compute 
the Jones polynomial. 

1.2 Methods Used In This Paper 

A priori, to count the solutions of the nonlinear partial differential equations (1.1) is a 
daunting problem. Our attempts to simplify this problem are based on three ideas. 

The first is a standard idea in topological field theory. We consider knots in W = 
R X C, where C (which may be simply M?) is a two-manifold, and we parametrize M 
by x^. We stretch our knots in the x^ direction, so that except at a few exceptional 
values of x^ where the number of strands changes, the boundary conditions are nearly 
independent of x^ (fig. 2). We hope that, away from the exceptional values of x^, the 
solutions can be approximated by solutions that are independent of x^. Once one drops 
x^, the equations reduce to equations in three dimensions. The reduced equations preserve 
more supersymmetry and one may hope to understand their solutions. 

After finding the three-dimensional solutions, to recover a four-dimensional picture, we 
have to take into account an adiabatic variation of the parameters in the three-dimensional 
equations. This is because our knot, even after stretching, is not quite independent of x^. 
We also have to consider the jumping that occurs when the number of strands changes. 
Actually, there is an important special case in which one only has to consider the adiabatic 
variation of parameters. This is the case (fig. 3) that M x C is replaced by x C (or I xC 
where / is a closed interval, though this introduces questions about boundary conditions) 
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Figure 2. Stretching a knot in one direction - here taken to be the direction - to reduce to a 
situation that ahuost everywhere is nearly independent of one coordinate. After much stretching, 
the knot is everywhere nearly independent of x^, except near the finite set of critical values of x^ 
at which a pair of strands appears or disappears. (In the figure, these occur only at the top and 
bottom.) 



X 



\ 




Figure 3. A braid in / x C; by gluing together the top and bottom, one can make a closed braid 
in X C. After much stretching, a braid can be described by adiabatic evolution in x^, with no 
exceptional values where this description breaks down. 



and the link is replaced by a braid. In this situation, one would study not the Jones 
polynomial but its associated braid group representations, which are also of great interest. 

One important thing to mention about this program is that it is not guaranteed to 
work. As one stretches a knot in the direction, the solution might simultaneously 
"spread" in the y direction (fig. 4) so that even after stretching, the solution might not 
approach an x^-independent limit. In fact, we will find that this happens under some 
conditions. One of our main tasks will be to understand conditions under which the sort 
of behavior suggested in fig. 4 does not occur. 

In carrying out the program that we have just described, we start in section 2 at t = ±1 
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Figure 4. As a knot is stretclicd along tlie boundary, a solution of the supersymmetric equations 
might become delocalized in the y direction, normal to the boundary. This is schematically indicated 
here; the shaded region indicates the spatial extent of a solution - that is, of the region over which 
the chosen solution deviates significantly from the one that describes the vacuum in the absence of 
knots - and its thickness is proportional to the amount that the knot has been stretched. 

because these are special values for Khovanov homology (as we recalled in section 1.1.1 
above) and also because some simplifications in the three-dimensional equations at t = ±1 
were already found in section 3.6 of [14]. 

Because we encounter some puzzling phenomena (which we will ultimately understand 
along the lines of fig. 4), we look for some additional simplifications. In doing so, we 
primarily exploit two ideas. 

The first idea is to modify the boundary conditions to incorporate gauge symmetry 
breaking. The basic idea was already explained at the end of section 1.1: instead of asking 
for (j) to vanish at infinity, we ask for — )• gag~^, where g ^ G and the three components 
of a = (01,02,03) take values in a Cartan subalgebra of the Lie algebra g of G. (In the 
more general case W = M x G, we would similarly modify the boundary condition to 
require that the component of (j) in the M direction is in a specified conjugacy class at 
infinity.) Continuously changing the boundary conditions in this way should not change 
the counting of solutions that leads to the Jones polynomial. On the other hand, in such 
counting problems one often finds that perturbing to a more generic situation can make 
things easier. Moreover, in the present case, taking a to be generic reduces the nonabelian 
gauge theory that we are studying to an abelian theory at low energies. If we scale up our 
knots so that all relevant directions of the Ki are large compared to l/|a|, then we can 
reasonably hope to find some sort of effective abelian description of the relevant phenomena. 

The second idea that we exploit is perhaps even more obvious. Since the equations that 
arise at t = 1 with the Nahm pole boundary conditions described above are rather special, 
we perturb the value of t and/or the Nahm pole boundary conditions to something more 
generic. This proves to be very fruitful, especially when combined with gauge symmetry 
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Figure 5. In a time- independent situation, wc look for solutions on a three-manifold M3 = C x M_|_, 
where C (taken here to be a two-sphere) is a Riemann surface. Knots are placed at points on the 
boundary of M3, labeled here as z\ , . . . , Z4. 

breaking. 

1.3 Outline And Results 

In section 2, we analyze the three-dimensional reduction of equations (1.1) at t = 1. We 
get an interesting description in terms of Higgs bundles with some additional structure, 
but it becomes clear that the program suggested in fig. 2 will encounter some difficulties 
at t = 1. In section 3, we perturb the equations to t 7^ 1 and find that this offers a 
much more promising framework for understanding the Jones polynomial. The equations 
for generic t have surprising and useful relations to a variety of topics in mathematical 
physics, including the Bethe equations for an integrable spin system known as the Gaudin 
model, and certain special "degenerate" conformal blocks of the Virasoro algebra; the rest 
of the paper is based on these relations. In section 4, we discuss the general framework 
for constructing braid group representations from adiabatic evolution of the parameters 
governing time-independent solutions. The general framework is a little abstract, but in 
section 5, we show that in our particular problem, it can be made very concrete using the 
free field representation of certain Virasoro conformal blocks. In section 6, we implement 
that idea in detail. This finally enables us to understand how the Jones polynomial and 
the braid group representations associated to it can be recovered by counting solutions of 
the four-dimensional BPS equations (1.1). Section 7 is devoted to describing an effective 
superpotential for BPS monopoles that can be used to understand some of the subtle results 
of sections 2 and 3. In section 8, we place some structures encountered in this paper in a 
wider context of mathematical physics. Three appendices fill in details of the derivations. 

2 Analysis At t=l 
2.1 Some Preliminaries 

As explained in the introduction, after stretching a knot along the first factor of M4 = 
M X C X ]R_|_, we want to find the solutions that are independent of the first coordinate 
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(which we cah x^), that is the solutions that obey reduced equations on the three-manifold^ 
M3 = C X M_|_. Here C is a Riemann surface (which may be simply M?), and M+ is the 
half-line y > 0. If knots are present, we take their support to be of the time-independent 
form M.X Zi X {0}, where the Zi are points in C, and {0} is the endpoint y = of M+. The 
picture is sketched in fig. 5. 

Solutions that can be derived from three dimensions preserve more than the generic 
amount of supersymmetry - they preserve four supercharges, to be precise - and accord- 
ingly, as described in section 3.6 of [14], their structure simplifies. In all cases that we will 
encounter in the present paper, the reduced equations can be usefully described in terms 
of three differential operators Pj, i = 1, 2, 3. The equations say that the Di commute 

[Vi,V,]=0, (2.1) 

and obey a moment map condition 

3 

^[A,P|] = 0, (2.2) 

1=1 

where is the adjoint of 2?,, in a natural sense. The commutativity constraint (2.1) is in- 
variant under complex- valued gauge transformations — t- gVig'^, where 5 is a Gc-valued 
gauge transformation (Gq is the complexification of G), while the moment map condi- 
tion (2.2) is only invariant under G-valued gauge transformations. What will make our 
problem tractable is that solutions of the combined system of equations modulo G-valued 
gauge transformations are equivalent to solutions of just the commutativity constraint 
(2.1) modulo Gc-valued gauge transformations. But solutions of the commutativity con- 
straint modulo complex gauge transformations can be described in terms of holomorphic 
quantities, so it is possible to understand them. 

Different instances of this structure vary by the construction of the Pj, which depends 
on the choice of t, and also on the boundary conditions that we assume at y = and at 
y = 00. The most basic case considered in [14] is that t = 1 (or —1) and the boundary 
condition is given by a Nahm pole in the part of (p tangential to the boundary. That 
boundary condition sets to zero (py, the normal part of (j), at y = 0. If we also require 
(j)y to vanish for y — t- 00, then a simple vanishing argument shows that in a solution that 
is independent of x^, (py is identically zero; similarly, Ai, the component of A in the 
direction, vanishes in a three-dimensional solution. Once (py and Ai are set to zero, the 
equations can be put in the above-described form, as shown in detail in [14], section 3.6,^ 

To minimize confusion, we note the following. In this paper, we use two different decompositions of 
M4 = R X C X R+ as the product of a three-manifold and a one-manifold: we write M4 = W x R+ with 
= R X C, but also M4 = R x Afs with M3 = C x R+. 

""Our notation differs from the notation used there by a relabeling of the coordinates — >■ x*"*"^ (whose 
purpose is to make "room" for a new time coordinate x'^ upon categorification) . Also, for later convenience 
we permute the Di in an obvious way. 
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with 



D . D 

D 



1^3 = ^-i[(bi,-] (2.3) 



and the moment map condition 

3 

= A^] = F23 - [h,4>3] - Dyct>i. (2.4) 

i=l 

In writing these formulas, we have simply taken C = M? with coordinates x'^,x^. However, 
it is helpful to introduce a complex coordinate z = x'^ + ix^, and to write cl)2dx'^ + (p^dx^ = 
(pdz + Tpdl; also we introduce a complex connection Ay = Ay — i(f>i for parallel transport 
in the y direction and write Vy = dy + [Ay , ■]. Then we can write 

D 

2?i = 2— 

Dz 

V2 = 2[ip, •] 

= (2.5) 
Vy 

With this way of writing the Pj, they make sense on an arbitrary Riemann surface C, with 
D/D'z understood as the d operator and (/? as a (1, 0)-form on C. 

As one would expect in a geometry that preserves four supercharges, the commutativity 
constraint = can be derived from a superpotential, namely 



W = -^/ T.^Fy-,, (2.6) 



with Tyz = [T>y, Dz]. To be more exact, in varying W to derive the conditions = 0, 

one should require the variation of j4 at y = to vanish. Otherwise, the variation of W 
contains additional delta function terms at y = 0. 

The equations {Vi,Vj] = — have been called the extended Bogomolny 

equations in [15]; actually, these equations are a hybrid of the equations of Nahm, Hitchin, 
and Bogomolny. They reduce to Nahm's equations if we drop the dependence on z, to 
Hitchin's equations if we drop the dependence on y, and to the Bogomolny equations if we 
set (f = 0. This is not just an analogy: we can borrow standard strategies from the theory 
of moduli spaces of Nahm, Hitchin, or Bogomolny equations. 

The equation [T>i,T>2] = 0, taken for fixed y, defines a Higgs bundle {E,ip) in the 
sense of Hitchin. The fact that Vi and 2^2 commute with P3 simply means that the 
Higgs bundle is independent of y, up to a complex-valued gauge transformation. When 
we specialize to C = = C, we get a Higgs bundle on C that can be understood as a 
Higgs bundle on CP"*^ = C U 00, possibly with a singularity at infinity. This case will be 
considered in section 2.5. There is another natural way for singularities of the Higgs bundle 
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to arise. As in section 6 of [14] , one may include surface defects supported on codimension 
two submanifolds in M4. Taking these to be of the form R x g^j x M+, where the qi are 
points in C, the construction summarized above still applies and the Higgs bundles acquire 
singularities at the points q^. We consider this situation in section 8.4. 

Additional structure arises from boundary conditions at ?/ = and y = 00. We will 
discuss the consequences of the Nahm pole at y = in section 2.2. The analogy to Nahm's 
equations will be useful: we can extract some holomorphic data from the commuting pair 
1)2, for any given point in C. As 1)2,1^3 commute with Di, this data varies holomorphi- 
cally on C, or possibly meromorphically in the presence of singularities. 

Finally, an analogy with the Bogomolny equations will help us understand the physical 
content of our solutions, especially when we turn on gauge symmetry breaking for y — )• 00. 
Indeed, the holomorphic data in the commuting pair Di,!)^ is analogous to the data which 
specifies the position of BPS monopoles in a solution of the Bogomolny equations. 

2.2 The Boundary Condition 

Our next task is to analyze the boundary conditions, first in the absence of singular 
monopoles. The model solution (1.3) at t = 1 has a singularity in T>2 and T>3, but not in 
T>i, simply because T>i does not contain the scalar fields. For g = su(2), with a standard 
choice of the Lie algebra elements t, this solution is explicitly 



We consider the matrices here to act on the fiber of a trivial rank two complex vector 
bundle E ^ M? x M^. For G = SO{3), on a general Riemann surface, there could be a 
global obstruction to defining as a rank two bundle and one would then consider instead 
the corresponding adjoint bundle ad{E). Because our considerations will be local along the 
Riemann surface C, a reformulation in terms of the adjoint representation does not change 
much, and we will omit this. 

The solution (2.7) can be written 





(2.7) 



<^ = gvi9 



-1 




9 



-1 



(2.8) 



with 




(2.9) 



and g a Gc-valued gauge transformation that is singular at y = 




(2.10) 



In other words, the solution is obtained by the complex- valued gauge transformation g 
from a trivial solution = y?!, A-^ = Ay = 0. 




We want to consider solutions which look like the model solution near y = 0, up to 
a gauge transformation. We would like to express this constraint in terms of the Higgs 
bundle data {E, ip) away from the boundary. For that purpose, it is useful to consider the 
behavior of a local holomorphic section s of the gauge bundle E that is invariant under 
parallel transport in the y direction. We will first do the calculation very explicitly for the 
model solution, and then identify which features are valid more generally. 

Let s be a section of the gauge bundle E that obeys 

^ = 0. (2.11) 

Since T) /Vy = g{d/diy)g~^ , the solutions of this equation are of the form s = gsQ, where 
■So is independent of y. Thus a general solution takes the form 

(2.12) 

with constants a, h. In particular, a generic vector in E, when parallel transported in the y 
direction to ?/ = 0, will blow up as y~^l'^ ■ There is a one-dimensional subspace consisting 
of solutions of Vs/Vy = that actually vanish as y^^^ for y — t- 0. This subspace is simply 
characterized by the condition a = 0. 

We write Ey for the restriction of £^ — ?• C x M4. to C x {y} for any fixed y > 0. Parallel 
transport using T>y gives a natural identification of the Ey for all y, and we write simply E 
for Ey, regarded as a bundle over C. Similarly, the restriction of (/J to C x y is independent 
of y, up to parallel transport by Vy. So by restriction to C x y, we get a Higgs bundle 
{E,^)^C. 

The "small" sections o{ E - the solutions s of Vs/Vy = that vanish for y — ?• - 
generate a rank one sub-bundle L C E. In the model solution, it is simply the sub-bundle 

of sections of E of the form ( ^ | . L is a holomorphic sub-bundle of — )• C; a section of 

vv 

L is holomorphic if it is annihilated by Vi = 2D/Dz. Concretely, in the model solution, a 
section ( ^ | of L is holomorphic if 6 is a holomorphic function of z. 

The fundamental reason that L is holomorphic is that, as [P2!^3] = 0, we can ask for 
a small solution of V^s = to also be annihilated by V2 = 2D /D'z. However, we cannot 
also ask for a small solution to be annihilated by ip = V\/2. This is clear from the above 
formulas; in the model solution, p annihilates a vector if the bottom component vanishes, 
not if the top component vanishes. 

To measure the failure of L to be (^-invariant, we can proceed as follows. If L is a 
trivial line bundle, which will be the case in our applications, then we can pick a section s 

of L that is everywhere nonzero. For our model solution, we just pick s = . Then we 

define 

K = sAips. (2.13) 
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For the model solution, we see that k = 1, and in particular k is everywhere nonzero. 
Nonvanishing of k means that (ps is not a multiple of s, so L is not invariant under 
multiplication by (p. (If is replaced by a general Riemann surface C, it might be 
impossible to pick an s that is globally nonzero, but one can still pick a local section s of 
L and measure the failure of L to be (^-invariant by computing k = s A ips. Whether k 
vanishes does not depend on the choice of s as long as it is nonzero.) 

Three basic properties of s which held for the model solution remain true for any 
solution with a regular Nahm pole, since they are unaffected by the subleading behavior 
of Vy and (/9 as y — )• 0: 

• A generic vector in E, when parallel transported in the y direction to y = 0, will 
blow up as y~^^^- 

• The sections s that under parallel transport to y = actually vanish span a rank 1 
holomorphic sub-bundle L C E. 

• Finally, we cannot also ask for a small section to be annihilated by = T>i/2. On the 
contrary, in a solution that can be approximated by the model solution near y = 0, 
(p does not annihilate a small section at any point in C. 

Consequently, a solution of the full equations with a regular Nahm pole for y — )■ 
gives not just a Higgs bundle {E,Lp). Rather, there is some additional structure: E is 
endowed with a holomorphic line sub-bundle L which is nowhere stabilized by ip. On M^, L 
is inevitably trivial, so we can simply say that an everywhere non-zero holomorphic section 
s oi E exists with s A(ps = 1. Notice that this condition is far from sufficient to determine 

ip. If s = , then the condition s A (ps = 1 fixes the upper right matrix element of (p 

and puts no condition on the others. The reason that only one matrix element of (p is fixed 
in terms of s is that when we make a gauge transformation that behaves like (2.10) for 
y — 7- 0, acquires a singularity that only depends on its upper right matrix element. In an 
appropriate situation (on a Riemann surface C of higher genus, or on M? in the presence 
of singular monopoles, as introduced shortly), there can be a nontrivial moduli space of 
triples {E, ip, L) with E and L fixed and only ip varying. This is explained in section 8. 

For any triple {E,ip,L), we expect that the full system of equations can be solved 
by a complex gauge transformation. Suppose we are given a y-independent Higgs bundle 
{E, ip), and a holomorphic sub-bundle L C E which is nowhere stabilized by (p. In a basis 
of E given by tps and s, we make the gauge transformation (2.10). This will reproduce the 
Nahm pole singularity at y = 0, but generically it will not give a solution of the moment 
map condition (2.4). By further making a smooth complex gauge transformation, one can 
hope to get a solution of the moment map condition. 

2.2.1 Adding Singular Monopoles 

In order to add singular monopoles at the boundary, we need to replace the Nahm model 
solution with a more general singular solution, given in [14], section 3.6. Let us consider the 
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case of a single singular monopole, located at y = z = 0. The model solution has the same 
singularity as before for y — >• at 2; 7^ 0, but has a more complicated form near y = z = 0. 
It can be obtained by a complex- valued gauge transformation, described explicitly in [14], 
from a solution of the commutativity constraint with 

^=(00)' ^^ = ^y = ^- (2.14) 



To express this in our present language, L is still spanned by sections that in the gauge 

'0^ 



(2.14) are multiples of s = [] ', this is because g has similar behavior as before for y — )• 0. 



We now have 

sAips = z'^. (2.15) 

The zero of s Afs is interpreted as the position of the singular monopole, and its degree is 
the charge. This interpretation suggests immediately what a solution with several singular 
monopoles should mean. We consider a solution described by a Higgs bundle (E, ip) with 
a sub-bundle L C E that is generically not (/9-invariant. If 



s 



s A y^s = Hiz - Za)''\ (2.16) 

a=l 

then we say that the solution has singular monopoles of charges ka at the locations Za- One 
hopes to be able to prove that given such data, there is a unique solution whose singularity 
near each z = Za agrees with that of the singular model solution. 

Though we mainly consider su(2) in the present paper, we can readily generalize these 
statements to a more general Lie algebra. For simplicity, take G = SU{N) and view E as a 
complex vector bundle of rank N. Consider a Nahm pole based on the principal embedding 
of 5u(2) — )• 5u{N). The eigenspaces of (in the fundamental n-dimensional representation 
of 5u{N)) are one-dimensional, and we have a line bundle L defined by sections which 
decrease as fast as possible as y — t- 0. The only constraint on L is that, away from the 
positions of singular monopoles, 

L®ipLe---®ip''-'^L = E (2.17) 

When we specialize to C = = C, we can define the line sub-bundle L by a specific 
section s of E, defined up to rescaling, such that the sections s, (ps, • • • , ip^~^s are linearly 
independent. This constraint is relaxed at the location of singular monopoles, in a way 
which depends on their charges. An equivalent description in terms of zeroes of matrix 
elements of ip is given in [14], eqn. (3.59). 

This can be extended to the case of a general Nahm pole at ?/ = 0, not necessarily 
associated to a principal embedding of su(2). (This extension will not be studied in the 
present paper.) In general, for y — )■ 0, there are local sections growing as either integer 
or half-integer powers of y. We can define a fiag of holomorphic sub-bundles En of E by 
looking at local sections which grow at most as a given power y('^~^y^ as y — >■ 0. Clearly, 
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En+i C En, also Eq = E, and En = for large enough n. Upon rescaling by a 
generic vector in En has a finite limit as ?/ — ?• 0, and this limit is an eigenvector of with 
eigenvalue {N — n)/2. The kernel of this map is En+i] hence the y — t- limit identifies 
the quotient spaces En/En+i with the eigenspaces of t^. Multiplication by if gives maps 
En — )• En-2- As ~ the holomorphic maps (pn ■ En/En+i — )■ En-2/ En^i can be 
identified with the action of t+ on the eigenspaces of t^. 

2.3 Solutions Without Symmetry Breaking 

Now we want to use the ideas that have just been described to determine some moduli 
spaces of solutions of the supersymmetric equations (2.1), (2.2). We will do this for G = 
SO{3), so that the dual group whose representations label the singular monopoles is G'^ = 
SU{2). 

First we work at the origin of the Coulomb branch; this means that we consider 
solutions such that the scalar fields (j) vanish for ?/ — >• 00. In particular, if must vanish at 
infinity. Since the equation Vyip = means that the conjugacy class of ip is y-independent, 
if can only vanish at infinity if it is everywhere nilpotent. This means that by a complex 
gauge transformation, we can make (p upper triangular: 




with some polynomial p{z). We cannot, however, put s in a standard form at the same 
time. So we simply take 

with polynomials P, Q. Without changing 99, we can make an upper triangular gauge 
transformation, shifting P by a polynomial multiple of Q: 

P{z) ^ P{z) + U{z)Q{z). (2.20) 

This is the only freedom, apart from a rescaling of s by a complex constant. 

Now let us ask how we can pick p, Q, and P to describe a configuration with singular 
monopoles of charges ka at the points Za in the boundary. Since s/\ips = pQ^, the condition 
that we need is 

piz)Q{zf = Y[{z - za)'^ := K{z). (2.21) 

a 

If the singular monopoles all have minimal charge, fe^ = 1, we can only obey this with p = K 
and Q = ±1. Then we can set P{z) = by a transformation (2.20), and changing the sign 
of Q multiplies s by an inessential constant. So at the origin of the Coulomb branch, there 
is a unique solution with singular monopoles of specified locations and minimal magnetic 
charge. 

This fact is actually an obstruction to the program that was described in section 1.2 
(see fig. 2). The case of singular monopoles of minimal charge is supposed to be dual to 
the Jones polynomial, which is the invariant computed in Chern-Simons gauge theory for 
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= SU{2) with all knots labeled by the two-dimensional representation of SU{2). If 
there is only one solution, this means that the physical Hilbert space associated to this 
problem is one-dimensional. The Jones representations of the braid group would then be 
of rank 1. This is certainly not the case. We must be running into some version of the 
problem that was indicated in fig. 4. We will get a clearer picture of what is happening 
after including symmetry breaking in section 2.4, and after deforming to t 7^ 1 in section 
3. This will ultimately enable us to circumvent the obstacle just described. 

What happens if some singular monopoles have non-minimal charge? In this case, we 
encounter moduli spaces of solutions. In general we solve (2.21) by 

p{z) = \{{z-Zar\ Q = \{{z-Zar, (2.22) 

a a 

with 

ma + 2ra = ka. (2.23) 

For Q as in (2.22), we cannot set P to 1 by a transformation (2.20); rather, the values 
of P and its first — 1 derivatives are invariants at each point z = Za- So there are ra 
moduli associated to each Za, where the possible values of are {0, 1, ... , [A;a/2]}. We will 
argue beginning in section 2.4 and in most detail in section 7 that the moduli represent the 
positions in the y direction of some smooth BPS monopoles, together with some conjugate 
angles. 

In implementing the program described in section 1.2, if some knots are labeled by 
integers fca > 1 - in other words, by representations of G"^ = SU{2) of dimension /cq + I > 2 
- one would have to handle the evolution of the four-dimensional solution as a path in 
the moduli space of three-dimensional solutions. (This would be done by quantizing the 
moduli space to get an appropriate space of physical states, in which the evolution would 
take place.) Instead of following that route, we will perturb the equations that we have 
just analyzed to more generic ones, with symmetry breaking at infinity or with t 7^ 1. This 
will eventually reduce all of our moduli spaces to finite collections of points. 

2.4 "Real" Symmetry Breaking 

Our first approach to getting a clearer understanding of the solutions - and a more useful 
reduction to three dimensions - will be to move on the Coulomb branch. The basic idea 
was already explained at the end of section 1.1. We fix some constant, nonzero expectation 
values a for the tangential scalar fields (p, and consider only solutions of the supersymmetric 
equations (1.1) such that lim^^oo 4' = 90'9~^i for some g ^ G. 

From our present point of view, we want to analyze the reduced three-dimensional 
equations in the context of symmetry breaking. The reduction splits off the direction, 
so we write (p- dx = (pi dx^ + ipdz + Tpdz. The effects of an expectation value for (pi or for 
93 will be quite different. We first consider the case of turning on (pi only. As we will see, 
this has the effect of making the solutions somewhat more physically transparent, without 
adding or removing solutions. 
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For G = SO{3), we can pick a gauge where (pi is constant at large y. This means that 
for large y, Ay = Ay — approaches a constant matrix at infinity, say such that 

A - (: ij ■ 

with ai > 0. This expectation value breaks SO{'i) to U{1). We pick a generator of the 
unbroken normalized so that the off-diagonal components of an adjoint- valued field, 

which are the fields of minimum electric charge for G = S0{2>), have charges ±1: 

q4-T'^ (2.25) 




The choice of sign will be convenient. 

The equation Vyi^ = implies that the lower-triangular component of ip must be zero, 
since otherwise it would grow exponentially fast at infinity. The diagonal component must 
also be zero, since otherwise 93 would have a nonzero limit at infinity. Hence symmetry 
breaking provides a natural frame in which (p is strictly upper-triangular for — t- 00, just 
as the boundary condition did for ?/ — )• 0. But (p decays exponentially fast at large y. 

In the low energy effective U{1) theory, a field of electric charge 1 is a section of a line 
bundle that we will call M. The first Chern class of A^, integrated over the plane given 
y = UOj for some large constant yQ, is an integer that we will call the magnetic charge 
m. The upper triangular matrix element of 



99 ~ ^ (2.26) 




is a section of A4~^, and hence the magnetic charge, which is minus the first Chern class of 
is minus the number of zeroes of p. This is the same as minus the number of zeroes 
of p (defined in (2.18)), which is equivalent to by a complex gauge transformation. So in 
the notation of (2.23), the magnetic charge is 



m 



Y,ma = Y,(-''a + 2ra). (2.27) 



This means that for a given configuration of singular monopoles or 't Hooft operators at 
y = 0, the smallest possible value of m is — ka- Every time that we add a zero to Q, 
m increases by 2. 

We propose that this fact can be interpreted in terms of smooth BPS monopoles. If we 
simply set (/? = 0, the extended Bogomolny equations that we have been studying reduce 
to the usual Bogomolny equations. Far from the boundary, the Bogomolny equations are 
a very good approximation, since ip is so small. The Bogomolny equations on admit 
smooth monopole solutions of charge 2. (We measure magnetic charge in units such that the 
minimum magnetic charge allowed by Dirac quantization is ±1. This is also the magnetic 
charge of a minimum charge singular monopole or 't Hooft operator, but smooth monopoles 
have even charge.) Our proposal is that if Q has a zero of order at z = Za, then there 
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are smooth monopoles located at z = Za- This statement has a precise meaning only 
if the monopoles are located at very large y, where the extended Bogomolny equations 
reduce to the ordinary ones, and have smooth monopoles as solutions. However, we have 
found that if Q has a zero of order ra at z = Za, then there are precisely complex moduli 
associated to the point z = Za in the complex 2-plane. We propose that these moduli are 
the positions in the y direction of smooth monopoles that are located at z = Za, along 
with conjugate angles. Again, the precise meaning of this statement holds when the y 
positions in question are large. We will make this interpretation quantitative in section 7, 
but for now we consider qualitative arguments. 

A first motivation for this proposal comes from m; we interpret the formula (2.27) as 
the sum of the magnetic charges of the singular monopoles at the boundary (— ka) plus 
a contribution of 2 for every smooth monopole (contributing 2^^ra in toto). The fact 
that the charge of the singular monopoles is negative is worthy of note. Prior to symmetry 
breaking, the charge of the singular monopoles does not have a meaningful sign - it is dual 
to a representation of = SU{2) - but after symmetry breaking, this sign makes sense 
and with our normalization it is negative. 

To learn more, we recall some facts about the Bogomolny equations. For the Bogo- 
molny equations on M^, the holomorphic data are the commuting operators Vi and V^, 
defined exactly as in eqn. (2.5). Localized monopole solutions of the Bogomolny equa- 
tions manifest themselves in the holomorphic data through "bound states" in the parallel 
transport by P3, i.e. through normalizable solutions to 



Such solutions only appear at specific positions in the z-plane, which are interpreted as the 
z values of the monopole locations, because generically, if we pick s to decay exponentially 
for y — )• — 00, it will grow exponentially for y — )• +00. 

In our setup, we are limited to y > 0, and an analogous normalizable solution exists 
precisely if the "small" section s that vanishes for y — )■ also decays exponentially at large 
y. If (f also decays exponentially as well, then s A(ps vanishes for y — >• 00 and hence for all 
y. So sections s that vanish at both ends can arise only only at zeroes of K{z). Indeed, if 




(2.29) 



near infinity, then given the form of (2.24), vanishing of s for y — )• 00 is equivalent to 
Q = 0. Zeroes of Q are the same as zeroes of Q. (Indeed, s A ips is independent of y; 
its zeroes are the zeroes of Q if y is small or of Q if y is large.) So the relation to the 
ordinary Bogomolny equations does indeed suggest that there are smooth monopoles at 
each point z = Za at which Q vanishes. The position of the smooth monopoles in the y 
direction should be encoded in the values of P and its first — 1 derivatives. We will make 
this picture more precise in section 7; for now we simply observe that the normalizable 
small section s{za) behaves as P{za)e~^^y , which suggests that increasing |P(za)| moves 
the smooth monopoles towards large y. 
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The limit P{za) — ?• is somewhat singular, as s is supposed to be everywhere non- 
zero. The physical picture suggests that the limit will correspond to a monopole bubbling 
situation, where the smooth monopole is pushed to the boundary, and screens the singular 
monopole's charge. Monopole bubbling is a phenomenon (originally described in [18] and 
rediscovered in section 10.2 of [15]) in which an 't Hooft operator absorbs a smooth BPS 
monopole, lowering the magnitude of its magnetic charge. 

2.5 "Complex" Symmetry Breaking 

Now we consider the case of symmetry breaking in (p. Suppose, for G = SO{'i), that the 
eigenvalues of ip for y — t- oo are ita. Then, as commutes with T>y, its eigenvalues are ita 
everywhere. By a complex gauge transformation, we can reduce to the case that (p \s a, 
constant diagonal matrix: 

Writing as usual = for the small section, we find s f\Lps = 2aPQ. After putting if 

in the form (2.30), we can still make a gauge transformation by diag(A, A ^), mapping 

P^XP, Q^X-^Q. (2.31) 

In the presence of singular monopoles of charges ka located at z = Za (and y = 0), the 
condition we want to satisfy is s A ifs = Y[t=ii^ ~ ^aY°' '■— ^{^)- This becomes 

2aPQ = K{z). (2.32) 

Solutions of these equations are associated to factorizations of K[z) and (modulo a trans- 
formation (2.31)) have no moduli. 

In the case of d boundary 't Hooft operators that all have fca = 1, corresponding to 
minimum magnetic charge, the number of solutions is precisely 2'^. The solutions corre- 
spond simply to the possible ways to distribute the factors of K{z) between P and Q. 
The number 2'^ has a natural interpretation. On the Coulomb branch, we might expect a 
minimum charge 't Hooft operator to have two possible states, with its magnetic charge 
being aligned or anti-aligned with the symmetry breaking. The two states correspond to 
a zero in P or a zero in Q. In the dual description by Chern-Simons theory with Wilson 
operators, a minimum charge 't Hooft operator corresponds to a Wilson operator in the 
two-dimensional representation of SU{2); such an operator again represents two quantum 
states, with positive or negative electric charge along the axis of symmetry breaking. 

To confirm that the 2'^ solutions correspond to two possible choices of the magnetic 
charge for each 't Hooft operator, let us compute the magnetic charges of these solutions. 
Suppose that P is of degree di and Q of degree d2, where di + d2 = d. The ratio P/Q does 
not depend on the normalization of the small section s. This ratio has electric charge —1 in 
the low energy abelian gauge theory; it is a section of the line bundle On the other 

hand, concretely, P/Q has di zeroes and d2 poles on the z-plane. {P/Q has neither a pole 
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nor a zero at z = oo, when understood as a section oi M.^^. In fact, by a complex gauge 
transformation, the solution can be put near 2; = oo in the form of the original Nahm pole 
solution (1.3), and in particular is independent of z. This trivializes Ai near z = 00 and 
makes P/Q independent of z.) So the line bundle Ai has degree d2 — di and hence 

m = d2 — di. (2.33) 

From this, we see that the 't Hooft operator at z = Za contributes either —1 or +1 to m, 
depending on whether we place the factor oi z — Za in P or in Q. 

To get some insight into why there are more solutions for a 7^ than there are for 
a = 0, consider a slightly more general ansatz for ip: 



^= : '''' . (2.34) 




For a 7^ 0, the polynomial p{z) can be removed by an upper triangular gauge transforma- 
tion, but now there is a smooth limit for a = 0. The condition s A (ps = K becomes 

d 

2aPQ+pQ^ = Y[{z-Za). (2.35) 

a=l 

For a ^ 0, p is irrelevant, as it can be eliminated by P — ?■ P — pQ/2a. For a ^ 0, the 
equation (2.35) has 2'^ solutions (modulo a complex gauge transformation that preserves 
the form of (p), corresponding to factorizations of K{z) as {2aP + pQ)Q. But all of these 
solutions have P ~ 1/a except for the one solution with Q = 1 (and P = 0) that we found 
already in section 2.3. 

We interpret this as follows. In the presence of a minimum charge 't Hooft operator 
at a boundary point z = Za, there is always a magnetic charge —1 localized near the 
boundary. The magnetic charge that is localized near the boundary is the same for all 
solutions because the form of the solution near the boundary is always given by a standard 
model solution (the one described in section 3.6 of [14]), independent of everything else. 
In the case of a solution of (2.32) with Q{za) = and (therefore) P{za) 7^ 0, there is in 
addition a smooth BPS monopole, with magnetic charge 2, located at z = Za and at a 
value of y that depends on P{za). So the total magnetic charge associated to z = Za is 
— 1 + 2 = 1. As a — )• 0, P{za) — )• 00 and the smooth monopole disappears to y = 00. 
In this way, all of the 2'^ solutions become equivalent for a — ?■ 0, even though they are 
different for a 7^ 0. As always, such a description in terms of smooth BPS monopoles is 
only precise in the limit that the monopoles are located at large values of y, so that the 
Bogomolny equations are a good approximation near their positions. As we explain most 
fully in section 7, this is the case exactly when a is very small (compared to the inverse 
distances l/(-Za — -Zfe) between the 't Hooft operators) so in particular the description by 
BPS monopoles becomes precise for a — > 0. 

In the opposite case that a is large compared to the inverse distances, the symmetry 
breaking is strong and the different singular monopoles on the boundary are so far separated 
that they do not significantly influence each other. In this case, the essential statement 
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is simply that a single 't Hooft operator of minimal charge, in the presence of symmetry 
breaking with a 7^ 0, has two possible states, in which it looks like an 't Hooft operator 
of charge 1 or —1 in the effective low energy theory. Given this, a system of d widely 
separated 't Hooft operators of minimal charge naturally has 2*^ possible states. 

2.5.1 Implications 

The counting of solutions for a 7^ circumvents a difficulty that we encountered in section 
2.3. With 2*^ classical solutions, the physical Hilbert space in the presence of d 't Hooft 
operators of minimal charge will have dimension 2'^. These 2'^ states are potentially dual to 
the states of d Wilson operators labeled by the two-dimensional representation of SU{2). 
And the representations of the braid group on d strands that are associated to the Jones 
polynomial can certainly be realized in a vector space of dimension 2*^. 

Given this, it is reasonable to expect that the stretching strategy sketched in fig. 2 of 
section 1.2 can work if a is generic, but that for a = 0, an attempt to simplify a classical 
solution by stretching in one direction leads to behavior along the lines suggested in fig. 4. 

More generally, in the presence of an 't Hooft operator on the boundary of charge 
ka, we would hope to find ka + 1 states of magnetic charge —ka, —ka + 2, —ka + 4, . . . , ka- 
Indeed, such an 't Hooft operator in SO{3) gauge theory is dual to a Wilson operator in 
SU{2) gauge theory associated to the representation of spin ka/2; this representation has 
dimension fca + 1, and its weights are as indicated. We get the right number of solutions 
with the right magnetic charges if in solving eqn. (2.32), we allow arbitrary factorizations 
of K(z) with the zeroes split between P and Q in an arbitrary fashion. Unfortunately, we 
do not have a simple interpretation of the factorizations for which both P and Q vanish at 

z = Za- After all, s = is supposed to be everywhere nonzero. It seems possible that 

\QJ 

the factorizations in which P and Q have a common zero should be interpreted in terms 
of monopole bubbling, that is, as solutions of the equations in the presence of an 't Hooft 
operator of reduced charge ka — 2, ka — 4, etc. To develop this idea in detail is beyond the 
scope of the present paper. 

Going back to the simpler case that all ka are equal to 1, since the space of physical 
states has a promising dimension 2"^, the next step could be to try to compute the action 
of the braid group and to extract the Jones polynomial. Before attempting that, we will 
describe a further deformation of the problem, which turns out to be illuminating. Among 
other things, with this deformation, we will get a nice behavior regardless of the charges 
of the 't Hooft operators. 

3 Analysis At General t 
3.1 Some Basics 

Still searching for a useful description for 't Hooft operators of arbitrary magnetic charges, 
we consider the possibility of deforming the equation that we are trying to solve to one 
that might behave more conveniently upon stretching along one direction. 
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A strong hint that there is a useful deformation comes by considering the relation of the 
extended Bogomolny equations (2.3), (2.4) to Hitchin's equations [20]. Hitchin's equations 
are equations on an oriented two-manifold C for a connection A and an adjoint- valued 
one- form (f): 

F-(t)A(l) = 

dyl ★ = 

dA0 = O. (3.1) 

The moduli space of solutions of these equations, up to gauge transformation, is a hyper- 
Kahler manifold Mh{G,C)- As a hyper-Kahler manifold, Mh{G,C) has a family of 
complex structures parametrized by a copy of CP^. The complex structures of M.h{G-, C) 
have a simple interpretation. They correspond to ways of splitting the three real equations 
in (3.1) into two parts: two real equations, which can be combined to a single complex 
equation; and a third real equation that is "orthogonal" to the first two. The complex equa- 
tion then describes the holomorphic data parametrizing 7W//(G, C) in one of its complex 
structures, and the remaining real equation is a moment map condition. 

An exceptional split corresponds to the case that the complex equation is made by 
combining together the last two equations in (3.1) to get the Higgs bundle equation 

dAV = 0. (3.2) 

A more generic split involves a complex parameter C,. For any set 

(We write here (j)z and instead of ip and Tp.) For any ^, the equation 

[pS,p|] = (3.4) 

is equivalent to two real linear combinations of the Hitchin equations (3.1). The possible 
complex structures on A^h(G, C) are parametrized by Q, where we add a point at infinity 
to the Q plane to make CP^. For every Q, one defines a complex structure in which 
the equation (3.4) is regarded as an equation governing holomorphic data; in this complex 
structure, the holomorphic variables are = — C^^^t^z and = A- + (cj), and the 
equation (3.4) is holomorphic in those variables. The third linear combination of the 
equations is regarded in complex structure as a moment map condition. 

For generic the equation (3.3) simply says that the complex connection A"^ = 
Aidz + .A^-dz is flat. Once the equation (3.3) is expressed in terms of A'', it has no 
explicit dependence on ^. Thus, M.h{G,C) when regarded simply as a complex manifold 
in complex structure without worrying about its Kahler metric, is independent of Q for 
generic C- The exceptional values of C are and oo. For example, for <^ — ?• 0, we must 
multiply by whence it reduces to —[4>zi ■]• Meanwhile, for C = 0? T^z = Dz- So the 
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C — 7- limit of eqn. (3.4) is the Higgs bundle equation (3.2). By similar reasoning, the 
— 7- oo limit of (3.4) is the complex conjugate of (3.2). 

The Higgs bundle equation has played a prominent role in our analysis because in 
section (2.1), we derived the holomorphic data at y 7^ from the equation [2?i,P2] = 0, 
which was none other than the Higgs bundle equation. (The moment map condition in that 
analysis condition was not the usual moment map equation for Higgs bundles - namely the 
first equation in (3.1) - but rather it was the three-dimensional equation (2.4).) At this 
point, we are led to wonder whether we could modify the construction that led to (2.3) so 
that the holomorphic data at y 7^ will be given not by a Higgs bundle but by the solution 
of a different complex linear combination of Hitchin's equations. It is in fact possible to do 
so. Indeed the path to doing so is not quite uniquely determined. 

There is one particularly nice parameter by which we can vary the underlying four- 
dimensional equations. This is simply the parameter t in eqn. (1.1). In addition to changing 
the equations, we should consider the possibility of changing the boundary conditions. Of 
the six fields {A, (p) (that is, the components of A and (/> that are tangent to the boundary), 
only (/) has a singularity at y = in the basic Nahm pole solution (1.3). However, once one 
drops the dependence on the spatial coordinates x, = 4 super Yang-Mills theory has an 
50(6) symmetry that rotates the six components of A and (p. Therefore, one can obey the 
classical Yang-Mills equations with a solution obtained by applying an SO{6) rotation to 
the Nahm pole (such a rotated Nahm pole was studied in [19], section 4). Of course, the 
rotation in general will change the unbroken supersymmetry. It turns out, however (see 
Appendices A and B for more detail), that as long as the rotation matrix is contained in a 
certain 5^7(3) subgroup of 5*0(6), the unbroken supersymmetry is unchanged, and if it is 
contained in a certain U (3) subgroup, then the unbroken supersymmetry changes in a way 
that corresponds to a change in the parameter t in the four-dimensional equation (1.1). So 
there is considerable freedom in rotating the boundary condition, with or without a change 
in t. 

After making a choice along the lines just described, the deformed equations and 
boundary conditions have a three-dimensional reduction that takes the familiar form of 
(2.1) and (2.2), but with a different definition of the Pj. Deferring most of the details 
to the appendices, we will summarize some formulas that arise in an illuminating special 
case. If one wishes to preserve the 50(3) symmetry of rotations of the boundary, then the 
rotation of the Nahm pole can only depend on a single parameter: the polar part of A must 
be a constant times the polar part of (p. It turns out that such a boundary condition 
is compatible with the four-dimensional equations (1.1), with a modified value of t that 
depends on (. 

In the reduction to three dimensions, we must now impose Ai — d^cpi = = (py. The 
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(3.5) 



Here, as described in the appendix, Ay is an appropriate linear combination of Ay and (pi. 
The condition that the Di commute must be supplemented by a moment map condition, 
which is also described in the appendix. As one would anticipate for a geometry that 
preserves four supercharges, the commutativity constraint = can be derived 

from a superpotential, which in fact is a multiple of the Chern-Simons function: 



To be more exact, in varying W to derive the equations = 0, one imposes a 

constraint on the variation of ^ at y = 0, so as to avoid delta function terms in the 
variation of W. 

The holomorphic data away from y = are now easy to describe. We recover the 
picture studied in section 2 if = 0, but as soon as C is nonzero (and not infinite), the 
three commuting operators Di simply describe a complex flat connection on C x The 
covariant derivatives for this flat connection are Dz = Di/2, D-^ = — P2/2C1 T^y = ^3- So 
commuting operators T>i simply describe a complex flat connection. We will call the flat 
connection A, irrespective of how it was defined in terms of A and (j). A is constrained by 
a moment map condition, which does not quite coincide with the most commonly studied 
moment map condition for a complex flat connection [21], though it is qualitatively similar; 
we expect it to have a unique solution for any (. 

The holomorphic data away from y = are hence simply a complex flat connection 
on C X M^, or equivalently, since this space is contractible to C, a complex flat connection 
on C. We will mainly be interested in the case that C is or CP^. In either case, C is 
simply-connnected, so a complex flat connection on C is trivial. One may wonder therefore 
how anything of interest can happen. The answer is that most of the structure of interest 
will come from the boundary condition at y = 0. In the case of symmetry breaking, there 
is also some interesting structure in the behavior at y = 00. 

3.2 Nahm Poles and Opers 

The first point is to understand how a complex flat connection can have a Nahm pole. The 
answer is that the pole appears in Ay and Az, but not in Az- The model example of a flat 
connection with a Nahm pole is 




(3.6) 




y 

Az = o 




(3.7) 
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This describes a flat connection as long as [*3,*+] = We are interested in the case that 
t+ = *i + i*2, where the tj, i = 1,2,3 generate a principal 5u(2) subalgebra of g. 

Being flat, this connection can be described by a formula d + ^ = gdg~^. For example, 
for 5u(2), we can take explicitly 



yV2 



(3.8) 




which leads to 

Alternatively, the model solution can be generated from the non-singular flat connection 

A.= ^^^^, A- = Ay = (3.10) 
by the same singular gauge transformation as in (2.10): 

9=\ n ..1/2 • (3-11) 





We now proceed rather as we did in the Higgs bundle case to explain the condition 
that must be placed on a complex flat bundle E so that it can be placed in the form (3.9) 
near y = 0, modulo less singular terms. We write simply E for the restriction of E to 
C = C X y for some fixed y > 0. We consider solutions s of the equation VyS = that 
vanish as y^^^ for s — t- 0. Sections s obeying these conditions span a rank one sub-bundle 
L G E. In the case of the model solution, any such s is a multiple of 

(3.12) 

so as in the Higgs bundle case, L is simply spanned by sections whose upper component 
vanishes. 

Also as before, if we regard E as a flat bundle over C, then L is a holomorphic sub- 
bundle; indeed, the object s that we have just defined is a holomorphic section of L, since 
it is certainly annihilated by Dz- However, it is not true that s is annihilated by T>z- On 
the contrary, a look at the previous formulas shows at once that 

sAVzS = l. (3.13) 

This brings us to the mathematical notion of an "oper." (For an explanation of this 
notion as well as a review of many related ideas that will enter our story later, see [22] or 
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[23].) For G = SU{2), an oper is a flat rank two complex bundle E bundle over a Riemann 
surface C, with structure group SL{2,C), together with a holomorphic line sub-bundle 
L C E with the following property: L is nowhere invariant under parallel transport by 
"Dz ■ The last statement means the following. If s is a local nonzero holomorphic section of 
L, then "DzS, which will be valued since is a connection on E, is nowhere L- valued. 
An equivalent statement, since L is spanned by multiples of s, is that DzS is nowhere a 
multiple of s. Alternatively, s A P^s is everywhere nonzero. The last statement does not 
depend on the choice of the nonzero section s, since if we replace s by fs (where / is a 
nonzero local holomorphic function on C), we have s A VzS — )• /^s A VzS. 

We have extracted an oper structure from the Nahm pole boundary conditions; con- 
versely let us see that given an oper on C, that is a pair {E, L) obeying the conditions 
just described, we can construct a solution of the Nahm pole boundary conditions. Go to 
a gauge in which L is spanned by vectors whose upper component vanishes. Holomorphy 
of L means that Az is lower triangular in this gauge. The oper condition s A P^s 7^ for 
any nonzero local section of s implies that in this gauge, the upper right matrix element of 
Az is nonzero. By a further diagonal gauge transformation, we can set this matrix element 
to 1. Then we pull back the flat bundle E with connection A from C to C x ]R+ (to 
get a flat connection on C x with no dependence on y) and make the singular gauge 
transformation (3.11). Having an upper right matrix element of 1 means that after the 
singular gauge transformation, Az has the singular behavior of the model solution (3.9); 
being lower triangular. Ay acquires no singularity. Finally, the gauge transformation gives 
Ay precisely the form of the model solution. 

So we have shown, at least for the case that G has rank 1, that two-dimensional 
opers correspond precisely to solutions of the Nahm pole boundary conditions in three 
dimensions modulo less singular terms. Conjecturally, by a further smooth complex-valued 
gauge transformation, one can satisfy the moment map condition. 

3.2.1 Some Further Remarks 

We add the following technical remarks. Since we want to be able to consider 't Hooft 
operators of minimum charge, we will take the gauge group in the rank 1 case to be 
G = 50(3), rather than SU{2). Accordingly, we should restate the above derivation in 
terms of the adjoint bundle ad(-E) rather than E. Because our considerations have been 
local on C, rewriting the construction in terms of the adjoint bundle does not change very 
much and we will omit it. (The main difference is that what can be naturally defined 
globally is in general not L but L^, which is a sub-bundle of a.d{E).) 

Also, everything we have said for G of rank 1 has an analog for any semi-simple G, 
somewhat as we indicated in the Higgs bundle case at the end of section 2.2. For example, 
for G = SU (n), an oper is a flat complex bundle E ^ C oi rank n with SL{n, C)-valued 
holonomies together with a line sub-bundle L C E that is holomorphic and has the property 
that if s is a local nonzero holomorphic section of L, then s, VzS, . . . , 'D^~^s furnish a local 
trivialization of E. The equations [T>i,T>j] = together with the Nahm pole boundary 
condition determine such an oper structure, by arguments similar to those we have already 
given. 
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3.3 Opers With Singularities 

Now we would like to modify the Nahm pole boundary condition to incorporate additional 
singularities - which we will associate with 't Hooft operators - on the boundary at y = 0. 

The type of singularity that we want can be guessed by analogy with the discussion of 
Higgs bundles. We will still have a flat Gc bundle E CxM_|_, and near a generic boundary 
point, the flat connection will look like the model solution (3.7), up to a unitary (that is, 
G-valued rather than Gc-valued) gauge transformation. We can still define a holomorphic 
line sub-bundle L <Z Ehj considering sections s obeying DyS = and vanishing for y — t- 0; 
and we still require that if s is a local holomorphic section of s, then s f\T>zS is generically 
nonzero. 

The only difference is that now we assume the existence of exceptional points , a = 
1, . . . , d, at which s AV^s vanishes. In fact, we specify positive integers ka and require that 
s AVzS vanishes in order ka for z — )• z^: 



This is analogous to requiring s A 93s ~ (z — Za)'^" in the Higgs bundle case. 

To get a precise problem of classical or quantum gauge theory with this sort of bound- 
ary behavior, what remains is to specify precisely what sort of singularity a solution of the 
moment map condition X^J^i'^ll ~ ^ (or the four- or five-dimensional equations that can 
be dimensionally reduced to it) is supposed to have at 2; = Za- As in most such problems, 
to do this one finds a model solution for the case of only one singularity, at, say, z = (and 
y = 0) and with an arbitrary k. Then one asks that the singular behavior of the solution 
near each of the points z = Za, y = should coincide with that of the model solution, for 
k = ka- For the case of Higgs bundles, the appropriate model solutions were found (for G 
of rank 1) in section 3.6 of [14], but for the generalization considered here, at present we 
are only able to find the model solutions numerically. They are described in Appendix C. 

The objects that we have described so far correspond to solutions of the fiatness and 
moment map conditions on C x with boundary conditions associated to Nahm poles 
or opers, except at finitely many boundary points where the oper condition is corrected. 
In particular, as soon as one gets away from y = 0, one simply has a flat bundle (with 
a moment map condition). The monodromy of the fiat bundle around the points z = Za 
is therefore trivial, so the exceptional behavior at the points z = Za only affects the oper 
property of the pair (E, L), not the flatness of E. Singularities of this kind are called oper 
singularities with trivial monodromy. 

3.4 Oper Singularities And Bethe Equations 

Let us now make concrete (referring to [22] for much more detail) what sort of an object 
is an oper with monodromy-free singularities. We will make this analysis for the case that 
C is simply-connected, so that there are no moduli in the choice of the flat bundle E ^ C. 
So G will be either or CP^; that is, it will be the complex z-plane with or without an 
added point at infinity. There are two reasons for assuming C to be simply-connected: this 
is the most relevant case for understanding the Jones polynomial; and also, eliminating 




(3.14) 
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the choice of E from the discussion will make it easier to focus on the opers and their 
singularities. 

Since C is simply-connected, a flat bundle over C is trivial. So we can go to a gauge 
with Az = A.-Z = 0. The line sub-bundle L of S is inevitably trivial for C = M?, or trivial 
after omitting the point z = oo for C = CP^. So it is globally generated by a section s, 
but we cannot put s in a simple form while also setting Az = Az = 0. Instead, we take 




(3.15) 

with polynomials P and Q. P and Q are uniquely determined up to a linear transformation 

MgGL(2,C). (3.16) 




Now s A VzS reduces to PdzQ - QdzP. So if the polynomial K{z) = Y\a=ii^ = ^o)^" 
encodes the positions and charges of the oper singularities, then the equation we would like 
to solve is 

PdzQ - QdzP = K{z), (3.17) 

modulo the action of SL{2,C). (Choosing K to have leading coefficient 1 has reduced 
GL(2,C) to 5L(2,C).) It is convenient to fix the action of two of the three generators of 
SL{2, C) by requiring that the degree of the polynomial Q is less^ than the degree of P 
(if two polynomials have the same degree, a linear combination of them has smaller degree 
and we call this Q), and that Q has leading coefficient 1, 

q 

Qiz) = ll{z-Wi), (3.18) 
i=l 

for some Wi. These conditions leave only the freedom to add to P a multiple of Q. 
We can recast (3.17) as 

P _ K{z) 
Q ~ ~~Q2 

The left hand side of this equation has zero residues at the zeroes Wi of Q{z). The right 
hand side must also have zero residues. This gives the constraints 



5.7^ = (3.19) 



' ^ = l.--->?- (3-20) 



Wi - Za '—f. W; 



Vice- versa, given a solution of these equations, the residues of K/Q'^ are zero; hence 
/ K/Q^dz is a rational function P/Q, and P is fixed up to a constant multiple of Q, 
which is the expected indeterminacy. 

If C = M^, opers with the desired monodromy-free singularities simply correspond to 
the solutions of the equations (3.20). For C = C¥^, we must further ensure that the oper 



^Later on, in the presence of symmetry breaking, we will have to relax this condition. 
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does not have an additional singularity at infinity. The condition for this turns out to be 
that the degree q of the polynomial Q is just one-half of the degree of K: 



q = \, k = Y,ka. (3.21) 

a 

To determine whether the oper has a singularity at infinity, let p be the degree of P and 
fp\ fp\ 

define ~ = z ^ \ ^] ■ We view P, Q as polynomials inv = 1/z. The condition that the 
\QJ \QJ 

degree of P exceeds the degree q of Q implies that there is no cancellation of the leading 
power of z on the left hand side of (3.19) and hence that p + q = k + 1. Since q < p, it 
follows that 

q < I (3.22) 

in general. The condition that the oper has no singularity at z = oo or v = is that 
(PdyQ — QdvP] 7^ 0, and it is not hard to see that this coincides with (3.21). 

To get farther, we need the theory of integrable systems. Rather "miraculously," the 
equations (3.20) are the Bethe equations of an integrable model, which is the Gaudin model 
or a certain large impurity limit of the XXX spin chain. (The connection between opers 
with monodromy-free singularities and the Gaudin model is reviewed in [22], following 
earlier developments such as [24, 25]. For more on this, see section 8.6.) For a = 1,. . . ,d, 
let Ra be a copy of the representation of SU{2) of spin ja = ka/2, and let Ti. = ®'^^iRa- 
The Hamiltonians of the Gaudin model are the commuting operators on % given by 

H, = Y.^^. (3.23) 

Here Ta are the generators of su(2) acting on and Tq • Ti, is the inner product of Ta 
and Tf, (defined with the quadratic form such that Ta-Ta = ja{ja + !))• Actually, what we 
have written in (3.23) are the Hamiltonians for the Gaudin model for = SU{2). (We 
call this group as it is naturally dual to the gauge group G that appears in the rest 
of our analysis.) There is a Gaudin model for any G^, and it bears the same relation to 
opers that we are about to describe for SU{2), but if has rank bigger than 1, then the 
Ha are only part of a complete set of commuting Hamiltonians. 

Since the Gaudin Hamiltonians commute, they can be simultaneously diagonalized. 
Moreover, since they commute with the action of G^, their joint eigenvectors can be or- 
ganized in irreducible representations of G^. Because of the G^ action, to understand all 
of the joint eigenvectors of the Gaudin Hamiltonians, it suffices to understand those joint 
eigenvectors that are also highest weight vectors for the action of G^. 

In the theory of the Bethe ansatz for the Gaudin model of G^ = SU{2), it is shown 
that solutions of the Bethe equations (3.20) correspond to the joint eigenvectors that are 
also highest weight vectors for the action of G^. In this correspondence, the weight w is 
related to the degree g of Q by 

w=^-q. (3.24) 
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In particular, if we want G^-invariant joint eigenvectors of the Gaudin Hamiltonians, we 
need w = and q = k/2; but as we observed in (3.21), this is the condition that the 
corresponding oper extends over CP^ with no singularity at infinity. 

So the number of opers on CP^ with monodromy-free singularities is the same as the 
number of 5?7(2)-invariant joint eigenvectors of the Gaudin Hamiltonians. But the joint 
eigenvectors of the commuting Gaudin Hamiltonians are a basis for 7i, and similarly the 
G^-invariant joint eigenvectors are a basis for T-L^^ , the G^-invariant part of %. So the 
number of opers that obey the conditions that we have imposed is precisely the dimension 

This result is our first concrete success in comparing the counting of BPS solutions in G 
gauge theory to Chern-Simons theory with gauge group G^. Consider Chern-Simons theory 
on CP^ with charges in the representations a = 1, . . . ,d. We place these charges at 
points Za G CP^. In the classical limit, the space of physical states is just the G^-invariant 
part of "H = ®aRa] tlie restriction to G^-invariant states is the Gauss law constraint. This 
also gives the right answer for the dimension of the physical Hilbert space of Chern-Simons 
theory if the Chern-Simons coupling parameter is generic. On the other hand, in the 
dual description in which the Hilbert space is constructed starting with time-independent 
solutions in G gauge theory, the states should correspond,^ from the arguments we have 
given, to opers on CP^ with singularities of charge ka at the Za- Since the number of these 
opers is the same as the dimension of 'HP^ , we have at least succeeded in reconciling the 
dimensions of the spaces of physical states in the two descriptions. This gives an indication 
that with the help of the deformation that we have exploited in the present section to 
C 7^ 0, the program of section 1.2 based on stretching a knot in one direction can actually 
work. 

The counting of states is less transparent if we take G = rather than CP"*^. Quali- 
tatively, it is clear that the number of physical states in Chern-Simons theory is larger on 
than on CP"*^, because, as the flux can escape to infinity, a physical state need not be 
completely gauge-invariant. However, to understand the condition that should be imposed 
at infinity is rather delicate, and it is hard to understand in G^ Chern-Simons theory 
the result that seems to come from the opers: physical states on correspond to highest 
weight vectors in %. After incorporating symmetry breaking in section 3.5, the comparison 
between the two descriptions will be simpler. 

3.4.1 Relation To Conformal Field Theory 

In arriving at the Gaudin model, we have accomplished much more than simply getting 
a number of classical solutions that is reminiscent of known constructions of the Jones 
polynomial. The Jones representations of the braid group can be described [26] as the 
monodromy of the Knizhnik-Zamolodchikov equations [27]. These equations are as follows. 
Express the usual parameter q that enters the Jones polynomial as q = exp(27ri/(fc^ -|- 2)). 
Let 5S be the space of distinct d-plets zi^. . . ,Zd G C. And let %* be the trivial bundle 

^We explain in section 4.1 why time-dependent instanton corrections in the G gauge theory do not afTect 
this counting of states. 
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over *B with fiber T-L = (^'^^^Ra- The Knizhnik-Zamolodchikov equations are the following 
system of equations for a section G of %*: 



The Ha are the Gaudin Hamiltonians (3.23). The Knizhnik-Zamolodchikov equations 
describe parallel transport of the section of H* with respect to a certain flat connection, 
which is implicitly defined in (3.25). To verify flatness of the connection, one uses the fact 
that the Ha commute and also the relation dH^/dza = dHa/dzi,. The solutions Q of the 
Knizhnik-Zamolodchikov equation are conformal blocks of two-dimensional current algebra 
with symmetry group C^; they are important in two-dimensional conformal field theory. 

Since opers with monodromy-free singularities correspond to a basis for H, we will, in 
our approach to the Jones polynomial, eventually be using gauge theory to construct a flat 
connection on the bundle H*; moreover, as we hope to recover the Jones representations 
of the braid group, this flat connection should be gauge-equivalent to the one deflned by 
the Knizhnik-Zamolodchikov equations. Actually, this tempting-sounding route is not the 
one we will follow. The very same Jones representations of the braid group have another 
realization in conformal fleld theory in terms of Virasoro conformal blocks for correlators 
of a product of degenerate fields and this will prove more useful. 

3.5 Symmetry Breaking Again 

Just as in section 2.4, we can gain some further clarity by moving away from the origin of the 
Coulomb branch. As always, we do so by turning on constant and commuting expectation 
values for cj) near y = oo. In the present context, this means that the connection form A 
does not vanish at infinity, but is a one- form with constant coefficients; moreover, these 
coefficients commute with each other. 

3.5.1 "Real" Symmetry Breaking 

First we consider the case that only Ay has an expectation value at infinity. This expecta- 
tion value arises from the value of (pi at infinity, and so just as in (2.24) we have 



where we can take ai > 0. 

The condition that Az and A-z should have no exponential growth at infinity tells us 
that they must be upper triangular in a gauge in which Ay looks like (3.26) for y — t- oo. 
Thus, near y = oo, the real symmetry breaking gives us a natural way to put the whole 
flat connection in a triangular form. 

An oper that is endowed with a covariantly constant reduction of its structure group 
to the group of upper triangular matrices - that is, to a Borel subgroup - is called a Miura 
oper. This notion is described in detail in [22]. Any oper bundle without monodromy 
can be given a Miura oper structure; in fact, there is a one-parameter family of ways to 




(3.25) 




(3.26) 
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do so. Concretely, if the rank two flat bundle E ^ C has trivial monodromy, then the 
associated bundle of CP^'s (whose fibers are obtained by projectivizing the fibers of E) 
also carries a flat connection without monodromy. Let us call this bundle B. Picking an 
arbitrary section of B over some given point p £ C and parallel transporting it, we get 
a covariantly constant section of B which turns the underlying oper into a Miura oper. 
This procedure introduces one complex modulus - the choice of a point in the flber of B 
over the starting point p. This means that a Miura oper without monodromy depends on 
a complex modulus. (When - as in section 3.5.2 - we introduce symmetry breaking in a 
complex direction, this modulus will disappear, because there will be no freedom to make 
a gauge rotation of Ay at infinity relative to Az-) 

The Bethe roots have a particularly nice interpretation in the case of a Miura oper. 
To explain this most simply, let us go back to the case that C = M? and use a gauge with 
Az = Az = 0. The behavior for y — )• oo singles out a sub-bundle L of the rank 2 bundle 
E that is invariant under parallel transport; we may call it a flat sub-bundle. In a gauge 
with the asymptotic behavior (3.26), L is generated by a covariantly constant section s 
that vanishes for y — ?• oo. After a complex gauge transformation to set Az = A- = 0, J is 
simply constant; we may as well take 




(3.27) 



On the other hand, the behavior for y — >• determines a holomorphic (not flat) sub-bundle 
L C E, generated as before (in a gauge Az = A- = 0) by 




(3.28) 



The choice of the Miura structure s gives a way to pick a natural linear combination of the 
components of s, namely 

sAs = Q. (3.29) 

The Q determined this way is not necessarily the one that we used in section 3.4, where 
we took Q to be a linear combination of components of s with minimum degree; now Q 
is simply determined by the Miura structure. With our new choice, the zeroes of Q have 
a simple interpretation: they are the points at which s is a multiple of s. In other words, 
the sub-bundle determined by the behavior for ?/ — )• oo is generically different from the 
sub-bundle determined by the behavior for y — )■ 0. The zeroes oi Q - which are called 
Bethe roots - are precisely the points at which these coincide. Another way to say the 
same thing is that the Bethe roots are the values of z at which there is a solution of DyS = 
that vanishes for both y — )• and y — )• oo. As we discussed in section 2.4, in the context 
of the Bogomolny equations one would say that there are smooth BPS monopoles at those 
values of z (and some values of y). 

For a concrete example, suppose that there are no 't Hooft operators at all. The flat 
bundle ii^ — t- is completely trivial and it has up to isomorphism a unique oper structure 
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with 

(3.30) 

In the absence of symmetry breaking, our convention that Q is the hnear combination of 
components of s with smaller degree leads to Q = 1, and hence (up to the freedom of 
adding to P a multiple of Q and rescaling it) P = z. 

In the presence of real symmetry breaking, the Miura structure gives a distinguished 
choice (3.29) of Q which has no reason to be a constant. If Q is not constant, then by 
adding to P a multiple of Q and rescaling it, we can set P = 1, so 

(3.31) 

for some w. The polynomial K = PQ' — QP' is 1, consistent with the absence of any 't 
Hooft operators. We have found, in the presence of real symmetry breaking, holomorphic 
data corresponding to a one-parameter family of solutions depending on the choice of a 
point € R^. This is a solution with no 't Hooft operator and a single Bethe root. We 
will call it a bare Miura oper. In the right context, when lifted back to four dimensions, 
we will interpret this solution later as a "string" that is localized at z = w and at a value 
of y that depends on C,. 

In general, the degree of K is at most one less than the degree of Q (this bound is 
achieved precisely if P = 1), so the number of Bethe roots is at most one more than the 
degreee k = J2a^a. of K. 

3.5.2 "Complex" Symmetry Breaking 

If we give expectation values at infinity to all components of (j), while requiring A to vanish 
at infinity, then the complex connection A is constant and diagonal for y — t- oo, 




(3.32) 

where a is a complex number. The factors of (" arise in the change of variables from ip to 
A. 

At infinity in the z direction, for any y, the solution reduces to the unique solution 
[18] of Nahm's equations which has a Nahm pole at y = and behaves as (3.32) at y = oo. 
In particular, the connection form is constant for z — )• oo with fixed y. For this form of the 
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connection, we can write the small holomorphic section s as 



(3.33) 



where a' (which equals a for y — ?• oo) is the constant value of A-z/C at large z with fixed y, 
and So has a finite limit at large z and fixed y. 

In order to analyze the holomorphic data in such solutions, it is unnatural to gauge the 
connection away. The information we want would be hidden in the behavior of the necessary 
gauge transformation at infinity. Rather, we will pick a complex gauge transformation 
which brings the connection exactly (not just asymptotically) to the form 



As we will now see, this can be done by a relatively simple type of gauge transformation. 
Consider any gauge field on x R+, such as one that has the asymptotic form discussed 
above, such that the connection form approaches a nonzero constant for z — )• oo with 
fixed y. If we were to compactify to CP^, we would say that such a connection has an 
irregular singularity at z = oo with a double pole. Connections with irregular singularities 
have Stokes phenomena. If one only considers gauge transformations with no essential 
singularities (that is, with polynomial growth only at infinity), the Stokes data is gauge- 
invariant. Two fiat connections with an irregular singularity are gauge-equivalent by a 
gauge transformation with only polynomial growth at infinity if and only if they have the 
same monodromy and the same Stokes data. 

As our connection has only has a double pole at z = oo, it has only two Stokes sectors. 
Stokes theory tells us that the monodromy around z = oo can be decomposed into the 
product of a diagonal formal monodromy matrix, and a sequence of Stokes matrices, which 
are alternatingly upper and lower triangular with ones on the diagonal. With only two 
Stokes sectors, the expression for the monodromy is 



with constants 6,6, and fi. As we are on M^, which is simply-connected, the monodromy 
at infinity must be M = 1. This together with the form (3.35) of the monodromy implies 
that the three factors in (3.35) - the Stokes matrices and the formal monodromy - must 
all equal 1. Hence we can bring our connection to the constant diagonal form (3.34) by a 
gauge transformation which grows only polynomially at infinity. 




(3.34) 




(3.35) 
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In particular, s will now take the form 



(3.36) 



with polynomials P and Q. The equation s A P^s = K{z) gives 



{Pd^Q - Qd^P) - = K{z) 



(3.37) 



which can be converted to the requirement that Ke^""^ jQ'^ has no residues at the zeroes 
Wi of Q{z). This becomes 



We are interested in solutions of (3.37) modulo a rescaling P — > AP, Q — > A~^Q (which 
corresponds to an automorphism of the diagonal flat connection), so we actually only care 
about the zeroes of Q{z). 

These are, again, Bethe equations, this time for the Gaudin model with an irregular 
singularity at z = oo [28-30] . It is shown in [30] that solutions of these Bethe equations are 
in one-to-one correspondence with eigenvectors of the spin chain, making in all na(^a + ^) 
solutions. It may also be possible to extract this result from the theory of the XXX 
spin chain, which is more familiar than the Gaudin model. The Bethe equations of the 
Gaudin model arise from a specific "large impurity limit" of those of the XXX chain. The 
deformation parameter a maps to the twist of the XXX chain, which is known to break the 
global S\J (2) symmetry to \J (1), and simplify the counting of solutions to Bethe equations: 
instead of a solution for each eigenvector of the spin chain Hamiltonian that is an SU (2) 
highest weight, one gets a solution for each eigenvector. 

As a simple example, consider the case of a single 't Hooft operator of charge 2, so 
K(z) = z^ . Just as in section 2.5.1, there is a solution of (3.37) with P = \ and Q a 
quadratic polynomial in z, and a solution with Q = 1 and P a quadratic polynomial. 
What happens if P and Q are both linear in z? Again as in section 2.5.1, we can solve 

(3.37) with P = Q = z^/— C/2a, but this solution does not correspond to an oper, since P 
and Q have a common zero. The novelty is that there is also an acceptable solution with 
P = — (C/2a)(z — C/^*)) Q = z + (/a, corresponding to a solution of the Bethe equations 

(3.38) with a single Bethe root. This gives a total of 2 -|- 1 = 3 opers obeying the necessary 
conditions. The last solution disappears if we take a = and has a common zero for P 
and Q if we take — t- 0. 

3.6 Opers And Stress Tensors 

In our analysis of opers with trivial monodromy, we have used a gauge in which the con- 
nection is trivial, Az = Az = 0. Correspondingly, we had to make a general ansatz for the 
small section s that generates the holomorphic sub-bundle L G E: 




(3.38) 




(3.39) 
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Here we will make a gauge transformation to put s in the standard form with upper 
component vanishing, and see what we can say about Az- This will have two benefits. We 
will begin to understand the relation of opers to conformal field theory. And we will get 
a description that is more general, not limited to the case (which however is particularly 
important in the present paper) of oper bundles of trivial monodromy. 

To put s in a simple form by a smooth gauge transformation would make Az and A- 
both nonzero. It turns out to be more helpful to keep A^ = 0. To also keep Az regular 
would require that our gauge transformation should be holomorphic, which is too restrictive 
a condition. Instead we will consider meromorphic gauge transformations, which will keep 
Az = 0, put s in a standard form, and generate poles in Az- 

The most obvious meromorphic gauge transformation that puts s in a standard form 

is 




^ = I " ^-1 1 • (3-40) 

We have 





hdzh-' = dz+{ , (3.41) 



where as before K = PQ' — QP'. It turns out to be more convenient to go to a gauge in 
which the upper right matrix element of Az is —1, at the cost of mapping s to a multiple 
of itself (this leaves unchanged the line bundle generated by s). We make a further gauge 
transformation by 

H=(''^ %-). (3.42) 



VK 

The possible double-valuedness of y/K is of no concern, for the following reason. If = 
50(3), so that G = SU{2), then i^T is a perfect square as all its zeroes are of even degree. 
If instead = SU{2), then G = 50(3), and we should really be writing all formulas 
in the adjoint representation, rather than the two-dimensional representation; accordingly, 
the sign of a gauge transformation is irrelevant. After a gauge transformation by h, Az 
takes the form 

where we have set 

K' Q' ^ ka/2 



V = 



y^^+y^— . (3.44) 

a I 

In the last step, we used K = Yiai^ ~ ^a)''" , Q = Yiii^ ~ ""^i)- Finally, a lower triangular 
gauge transformation 

/ 1 n\ 

(3.45) 
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leads to our final result for Az- 

A = r I , (3.46) 




with 

t = -v'-v'^. (3.47) 
In general, t has poles at both the Za and the Wi. Near z = Za, 

^ ' ^"^^^ + . . . , ja = ka/2. (3.48) 



{z - ZaY Z- Za 



Near z = Wi, 



t 



— ihj \ ^ — ^ — ^ — ^ inj — in- I 



Z - Wi \^ Wi - Za Wi-Wj 



Thus, t has no singularity at z = Wi if and only if the Bethe equations (3.20) are satisfied. 
To get some more insight, set = dz + [Az, ■] and look for a flat section, that is a 

ff\ 

holomorphic solution of Vz ( ~ 1 =0. We find that f = f and 

(^ + *(^))/ = 0. (3.50) 

We have carried out this derivation using a particular local coordinate z, but the notion 
that we started with - a flat bundle with an oper structure - did not depend on the local 
coordinate. So eqn. (3.50) must be covariant under a change of the local coordinate, 
with a suitable transformation for t. A short calculation (or a more careful study of 
the above derivation) shows that under a change of local coordinate z — )• z , and a suitable 
transformation oft, the object t transforms like a stress tensor in two-dimensional conformal 
field theory. In other words, it transforms not as a quadratic differential, as one might 
naively think from its pairing with the second derivative d^/dz'^ in (3.50), but with an 
"anomalous" term involving the Schwarzian derivative ^{z, i}. The double pole in t at 
z = Za is as if there is a primary field inserted at Za- 

If we had started from Az = C~^diag(a, —a) and used the same sequence of gauge 
transformations, with K = PQ' - QP' - 2aC~^PQ, we would have arrived to the same 
formulas, but with an extra constant term in v. 

^ = + (3.51) 

C^Z-Za^Z-Wi 

a, I 

Now t{z) has a pole of order four at z = oo. This would correspond in conformal field 
theory to the insertion at infinity of a somewhat unusual operator [36]. 
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Finally, the Bethe equations have the following interesting property. They describe 
stationary points at fixed Za (and a) of a Yang- Yang function: 

W{Wi,Za) = -^log((wi - Wjf) + ^^ka\og{Wi - Za) - ^^^kakb\og{{Za - Zbf) 
i<j i,a a<b 



2a -r-^ a s—^ 



(3.52) 

In other words, the Bethe equations can be written as 

^ = 0. (3.53) 

OWi 

The Yang- Yang function has another interesting property: the coefficients of the single 
poles in t at z = Za - sometimes called accessory parameters - are given by 



dza 



(3.54) 



Some terms in W which are independent of the Wi have been included to insure that this 
relation is satisfied. A similar relation holds for adaW. 

3.7 Opers and Virasoro Conformal Blocks 

In this section, we will show how these formulae arise naturally in the semiclassical limit 
of Virasoro conformal blocks. The semiclassical limit is defined as a limit in which the 
central charge c of the Virasoro algebra goes to infinity, while the conformal dimensions of 
operators also scale in the same way as c. 

There is a useful way to parametrize the central charge: c = 1 + 6Q^, with Q = 
b + b^^. In this parametrization, we take 6 — )• to get a semiclassical limit. The conformal 
dimensions of operators are conveniently parametrized as A = a{Q — a). The parameter 
a is often referred to as "momentum." We keep ba = rj fixed as 6 — )■ 0. Then the insertion 
of an energy- momentum tensor T in a correlation function scales as and we can define 
the finite limit t = b^T. We propose to identify this t, inserted in certain conformal blocks, 
with the t of section 3.6. 

The quantum stress-tensor T{z) has an anomaly under conformal transformations; it 
shifts by a multiple -(^{z, z} of the Schwartzian derivative. Hence t = b'^T has a conformal 
anomaly that is independent of b for 6 — ?• 0. The behavior of t near z = Za 'va the previous 
section corresponds to the behavior near a Virasoro primary field with = — These 
operators are very special: correlation functions and conformal blocks which involve only 
operators of this type can be described very easily by a free-field realization. We can 
describe such a realization in close parallel to the discussion in the previous section. 

Let X be a two-dimensional free field with two-point function {x{z)x{z')) = — ^ ln(z — 
z'). A standard way to construct an energy-momentum tensor of central charge c = l + GQ^ 
is to take 

T = - : dxdx ■■ +Qd\. (3.55) 
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If we define v = —bdx, this definition reduces to (3.47) in the hmit 6 — )• 0. This is the first 
hint that a free-field realization can be useful for us. Operators of dimension A = a{Q — a) 
can be readily described as normal-ordered exponentials of the free boson, 

Va{z) =: e^"'^^") : . (3.56) 

A second hint comes from (3.44): the quantity v is the semiclassical limit of the 
expectation value of —bdx iii the presence of chiral vertex operators of momenta —ka/2b 
at z = Za and of momenta 1/b at z = Wi. The operators Vi/!,{wi) have dimension 1. They 
are usually called screening operators in free-field realizations [32, 33], and are naturally 
integrated over curves. The singular part of the stress tensor near the location of a screening 
operator 

T{z)Vi/b{wi) ~ ^ .2 ^i/b{wi) + — duj,Vi/b{wi) H = dy,^ — Vi/b{wi) 

\Z Wi) Z Wi \Z Wi 

(3.57) 

is a total derivative, and drops off upon integrating over the position of the screening 
operator. 

We can easily compute the following free field correlation function: 

(U^^M^^^U^-kaM^a)] = Uiw^ - W,)-^ Hiw, - Z^p Yi^Za " Z^Y^'"^''' 
\ i Ob I free i<i i,a a<h 



exp ( ^y\>{w,z\ 



(3.58) 



The exponent on the right is the Yang- Yang function (3.52)! (For the moment, the terms 
proportional to a are absent as we have not included symmetry breaking.) 
Consider the integral 

lj\y~k^l2b{zS) = f (ll^l/biw^)llV_,^/2biZa)^ Hdw,, (3.59) 
\ a, / r ^ \ i a I free « 

where T is any integration cycle for which the integral converges. Because of (3.58), this 
is equivalent to 

(^y-K/2b{^a)^ = j^^My^{Za,W,)lb')\{dw,. (3.60) 

A Virasoro conformal block for the expectation value of a product of primary fields is 
a candidate correlation function that is compatible with the Virasoro Ward identity: 

(n.)Y[V_,^„A.^^ =('^_A_ + _2_|_)(nv_,./,,(.j) (3.61) 



r 



The functions defined in (3.59) or (3.60) have this property for an arbitrary choice of 
the number of w^s and the integration cycle F; this is proved using the definition (3.59) 
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and the fact that the screening charges are primary fields of dimension 1. (The function 
{'^{■^)Y\a^-ka/2b{za))-p defined by the integral (3.59) with an insertion of T{z) in the 
free field correlation function on the right hand side.) What we have just described is the 
free- field realization of the conformal blocks [32, 33]. 

The space of possible integration cycles T for the integral (3.60) has a natural fiat 
connection (the Gauss-Manin connection) as the points Za, a = 1,. . . ,d vary. Hence the 
functions (Ha ^-ka/2biza))r - for any fixed number of u;'s - furnish a representation of the 
braid group on d strands. It is known [34, 35] that these are precisely the representations 
of the braid group that are associated to the Jones polynomial and its generalizations. 
This relation of the Jones polynomial to conformal field theory will be more useful for 
the present paper than the relation via the Knizhnik-Zamolodchikov equation, which was 
noted in section 3.4.1. 

The functions (Jla ^-fca/2fe(-^a))r are not all possible conformal blocks for a product 
of primary fields with the dimensions of the V-ka/2b: rather, they are all such conformal 
blocks if the operators l^_fc„/2fe are degenerate primary fields in the sense introduced in [31]. 
Alternatively, these functions are all possible conformal blocks if the oper derived from the 
small b limit oit = b^T is supposed to have trivial monodromy at the points z = Za- These 
concepts and the relation between them are described in section 3.7.2. 

3.7.1 The Irregular Case 

Now let us incorporate the complex symmetry breaking parameter a in the above discus- 
sion. We can certainly in the integral (3.60) over the w's modify the exponent to include 
the terms proportional to a in the Yang- Yang function (3.52). But what does this mean 
in conformal field theory? We need to replace the free field correlation function (3.58) by 

(ll^l/b{m)llV^,^/2tiZa)) 

\ i a I free 

= n(^^ - ^^O"^ n(^^ - ^«)^ W^'^ - z,)-i^''^^^e^^'^^'"^-^^'^\ (3.62) 

i^i i,a aj^b 

What is the conformal field theory interpretation of this formula? 

Almost by construction, the right hand side is the free-field correlation function of 
the given product of fields with peculiar boundary conditions for x at infinity, x ~ az/bQ. 
Alternatively, we have inserted at infinity an "irregular vertex operator" , i.e. the L — )• oo 
limit of 

exp {-'^-^L^dxiL)^ ■ (3.63) 

In the presence of such an irregular vertex operator, the stress-tensor has the expected 
degree four pole. Nothing changes in the above formulae, except that the choice of possible 
integration contours is enlarged. The result of the integral is a conformal block with an 
irregular puncture at infinity, as defined in [36], in addition to the standard punctures 
of momenta —ka/2b. In the context of free fermions, operators associated to irregular 
singularities were originally defined in [37, 38] . 
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3.7.2 Degenerate Primary Fields And Trivial Monodromy 

Now we will review some standard facts about representations of the Virasoro algebra. 
This will enable us to explain what is special about the particular conformal blocks that 
come from the free field representation. 

For a generic value of the conformal dimension A, the Verma module defined as the 
span of all possible Virasoro descendants of a highest weight vector of dimension A is 
irreducible. For a set of special values 

(r-l)b s-1 , , 

a = ar,s = -^^ r,s = 1,2,3,... (3.64) 

or a = Q — Or^s, this is not true: a certain descendant at level rs is again a highest weight 
vector, and has zero norm. In a unitary conformal field theory, the descendant in question 
will vanish, and even in a non-unitary theory, it might vanish. The primary field whose 
descendant vanishes is called a degenerate primary field. We write ^r,s for such a field. 

The vanishing descendant of a degenerate primary field will certainly decouple in cor- 
relation functions. We call conformal blocks obeying such a relation degenerate conformal 
blocks. They satisfy a condition known as a "degenerate fusion rule." (It can be proved 
using the differential equation that we mention shortly.) In the OPE of an operator ^r,s 
and an operator of momentum a, only operators of momentum — — appear, with 

r' = r — l,r — 3,-- - ,1— r s' = s — 1, s — 2 • • • ,1 — s. (3.65) 

A conformal block with an insertion of momentum a^^s satisfies the null-vector decou- 
pling condition if and only if its correlation functions satisfy a certain differential equation 
of order rs. An important special case is r = 2, s = 1; the equation is 

92$2,i(^) + 6^ : r(z)$2,i(^) := 0. (3.66) 

This reduces in the semiclassical limit to the differential equation (3.50) associated 
to an oper. (For 6 — )• 0, the normal ordering in (3.66) is irrelevant; the only part of 

that survives for 6 — )• is the response to the "heavy" fields with momenta of order 
1/6.) Moreover, in the semiclassical limit, the expectation value oi t = b'^T will always be 
such that the monodromy of the differential equation around points with additional 
insertions is trivial. Indeed, the degenerate fusion rule implies that the OPE of $2,1 and 

contains only one primary field $2,s; from this, it follows that the monodromy of $2,1 
around a <I>i^s puncture is trivial. This is a quantum version of the trivial monodromy 
condition on the differential equation (3.50) associated to the oper. On the other hand, 
the singularity of t near a ^i^s insertion is precisely that which we have exhibited in (3.48) 
(with j = s/2). The upshot of this is that the semiclassical limit of a conformal block 
for a correlation function ^0(1=1 ^i.feaC-^a)^ determines an oper with precisely the sort of 
monodromy-free singularities that we extracted from three-dimensional gauge theory in 
section 3.3. 

The conformal blocks constructed from the free field formula (3.59) describe correla- 
tion functions of degenerate primary fields, simply because the free field vertex operators 
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Figure 6. A snapshot at fixed time of a time- independent situation. In the three-manifold M3 = 
C X R+, knots are present at boundary points zi, . . . , ^4. 

of momenta — 1| are degenerate. Related to this, it is possible to show that the confor- 
mal blocks which are given by the free-field realization do satisfy the degenerate fusion 
constraints. The corresponding differential equations are equivalent to the Picard-Fuchs 
equations satisfied by the free- field integrals, or to the natural flat connection on the space 
of integration cycles. They are a close analogue to the Knizhnik-Zamolodchikov equa- 
tions. The free field realization gives all the conformal blocks for the correlation function 
^0(2=1 *^i,fca(-^a)^ that are allowed by the fusion rules, so there are no more to be had. 

The interpretation of opers with monodromy-free singularities in terms of correlation 
functions of degenerate conformal fields gives an intuitive explanation to the Bethe equa- 
tions. Naturally, t should have no poles at the points Wi, because no conformal fields are 
inserted there. 

4 Four-Dimensional Solutions and Parallel Transport 
4.1 Introduction 

We now turn to the problem of analyzing time-dependent solutions of the original BPS 
equations (1.1). Even though this is a problem of classical partial differential equations, a 
quantum mechanical view is helpful. 

We start with a time-independent situation - twisted = 4 super Yang-Mills on a 
four-manifold M4 = M x Afa, where M3 is a three-manifold and we think of M as the time 
direction. In our application, M3 = C x with C a Riemann surface. If knots are 
present, we assume initially that they are time- independent. In this situation, which is 
depicted again for convenience in fig. 6, we want to find the BPS states - quantum ground 
states. 

The first approximation, already analyzed in sections 2 and 3, is to find time-independent 
classical solutions. In going from classical solutions to BPS states, we will ignore the non- 
compactness of M3. This means that we will ignore the existence of a continuum of non-BPS 
excitations. 
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If there is only a finite set of classical solutions of the BPS equations (and they are 
nondegenerate ~ there are no zero modes in expanding around such a solution), then the 
classical approximation to the space of BPS states is very simple. Let X be the set of 
classical solutions. Then for each I G X, there is in perturbation theory a quantum ground 
state '(1)1 that is localized near /. In perturbation theory, the form a basis for the space 
U of BPS states. 

Nonperturbatively, in problems of this general type, instanton effects might lift some 
of these approximate ground states away from zero energy. However, we are actually 
here dealing with a problem in which this does not occur. This is because in the time- 
independent case, even with knots present, our problem has four supercharges, and an 
instanton (that is, a classical solution with non-trivial time dependence) violates at least 
two of them. This leads to the existence of two fermion zero modes in an instanton 
background, which is one too many to contribute to a matrix element of the supercharge 
Q between approximate ground states V'/ and V'j- 

4.1.1 Relation To Morse Theory 

A more explicit understanding of why instantons do not lift the classical vacuum degenera- 
cies comes from the relation of this problem to Morse theory.^ Supersymmetric quantum 
mechanics related to Morse theory [42] is, in general, a theory of maps from M to ZY, where 
is a Riemannian manifold with metric tensor g endowed with a real-valued function h 
that we call the superpotential. For generic U. and /i, the model has two supercharges, one 
of which is conjugate to the exterior derivative: 

Q = e^de-''. (4.1) 

The classical vacua correspond to critical points of h. If /i is a Morse function - that is, 
its critical points are all isolated and nondegenerate ~ then in perturbation theory, each 
critical point I corresponds to an approximate quantum ground state ■0/- The fermion 
number qj of ipj is equal to the Morse index of the critical point I (the number of negative 
eigenvalues of the matrix of second derivatives of the function h at /). Since Q increases 
the fermion number by one unit, quantum corrections inducing non-zero matrix elements 
i'^jlQl'^i) are possible only if 

qj = qi + 1- (4.2) 

Such nonzero matrix elements can be computed by counting, in a suitable sense, the instan- 
ton solutions that interpolate between the critical point / in the far past and the critical 
point J in the far future. The relevant "instanton" equations, in other words the conditions 
for a map M — >■ Z// to be Q-invariant, are the gradient fiow equations of Morse theory: 

The problem we are studying of twisted A/" = 4 super Yang-Mills theory on M x M3 
(we primarily take M3 = C x M_|_ but the following remarks are more general) is an infinite- 



'For a relatively accessible introduction to the relevant aspects of Morse theory, see [39]. 
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dimensional problem of this sort,^ with lA being the space of complex- valued connections 
on M3. We view the Chern-Simons function 

W=-^ [ Tr ( ^ A d^ + A ^ A ^ ) (4.4) 
4vri J Ms \ 3 J 

as a holomorphic function on the complex manifold U. Holomorphy means that the one- 
dimensional sigma model with target U and superpotential W has four supercharges (this 
actually depends on the fact that the metric of U is Kahler and is also true in the gauge- 
invariant case mentioned in footnote 9). We actually want to study this model in the 
context of a twisting of A/" = 4 super Yang-Mills theory in which a particular supercharge 
Q is distinguished. This supercharge depends on a twisting parameter t [15] and is an 
infinite-dimensional version of Q = e^de~^ where h is the ordinary Morse function 

/i = Re(e^"W), (4.5) 

and 



1 — sin a 
t = 

cos a 

For a Morse function of this type, the gradient flow equation becomes 



(4.6) 



where the are local holomorphic coordinates on U. 

A down-to-earth manifestation of the relation of our problem to Morse theory is that 
the underlying four-dimensional supersymmetric equations (1.1) are the gradient flow equa- 
tions (4.3) for the Morse function h on the infinite-dimensional manifold lA. (This is one 
of the main ideas in [40, 41], where the gauge invariance mentioned in footnote 9 has been 
taken into account.) 

Now we can give a more explicit explanation of why nonperturbative effects in our 
problem will not spoil the supersymmetry of any of the approximate quantum ground 
states In general, for a Morse function that is the real part of a holomorphic function, 
isolated critical points all have the same (middle-dimensional) Morse index and thus the 
same value of the fermion number q. Hence the condition (4.2) is never satisfied. 

Another route to the same result is as follows. In general, gradient flow for a Morse 
function such as h that is the real part of a holomorphic function has a conserved quantity, 
namely the imaginary part of the relevant holomorphic function, in our case 

i = Im(e*"W). (4.8) 



® To be more precise, our problem is a gauge-invariant version of a such a problem - corresponding to 
a supersymmetric sigma-model with target Li coupled to gauge fields that gauge a symmetry of U. The 
gauge group in our case is the group of maps from A/3 to the finite-dimensional group G, while U is the 
space of complex- valued connections on A/3. However, in our problem the gauge group acts freely on W; 
this is ensured by the Nahm pole boundary condition. As a result, the gauge-invariance will not play a 
major role. In effect, for our purposes, we can replace U by its quotient by the group of complex gauge 
transformations and reduce to the case that there are no gauge fields. 
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For generic a, all critical points have distinct values of j and hence distinct critical points 
cannot be connected by a solution of the gradient flow equation. Hence there are no 
instantons that might spoil the supersymmetry of the states ^l'J. 

4.1.2 Time-Dependence 

As explained in section 1.2, we do not literally want to consider a time-independent situa- 
tion; rather, we want to allow for a slow time-dependence of the positions of the knots. Let 
*B be the space of distinct points zi, . . . , £ C. The space Ti of BPS states is the fiber of a 
bundle T-L over *B. This bundle carries a natural flat connection. This is a general property 
of topological fleld theory, but the Morse theory interpretation leads to a particularly nice 
description. 

In supersymmetric quantum mechanics related to Morse theory, since the supercharge 
Q = e^de~^ is conjugate to the exterior derivative d, the ground states ipj associated 
to critical points must have an interpretation in terms of the cohomology or dually the 
homology of U. There is a standard way to understand this in Morse theory. To a critical 
point /, one associates the downward-flowing cycle consisting of all points in U that can 
be reached by gradient flow starting at /. In other words, one considers solutions of the 
gradient flow equation on a half-line (—00, 0], with initial conditions that the flow starts at 
/ at t = —00. parametrizes the values at t = of such flows. 

The case that h is the real part of a holomorphic function has special features, and 
is particularly simple, so let us focus on that case. We make the further simplifying 
assumption that the critical points of h are isolated and irreducible. Vl is inevitably not 
compact (or it would not admit a non-constant holomorphic function). In this very special 
situation, the J} are called Lefschetz thimbles. The thimbles JT/ are closed if the angle a 
used in deflning h is sufficiently generic (to prevent the existence of gradient flows between 
distinct critical points), but they are not compact. So they do not represent classes in the 
ordinary homology ofVl. However, as J^j is deflned by downward gradient flow with respect 
to h, one has h — )■ —00 at infinity along J'j. As a result, the J^i are elements of a certain 
relative homology group - the homology H{U,U<:) of V( relative to the region ZY< where 
h goes to —00. (In the notation, we do not indicate the dimension of a homology cycle, 
because this relative homology is nonzero only in the middle dimension. That is related to 
the fact that the critical points all have a middle-dimensional Morse index.) 

The space Ti of supersymmetric ground states can be identified with the relative ho- 
mology H{U,V(<:). In this correspondence, the quantum ground state ipj associated to a 
critical point / maps to the element J'j of H(l/(,U<^). For an explanation of this from a 
physical point of view (in the context of supersymmetric quantum mechanics related to 
Morse theory), see [43] or [41]. 

The interpretation in terms of relative homology means that Ti has an integral structure 
and hence a natural flat connection. To give it a fancy name, this flat connection is the 
Gauss-Manin connection on the relative homology. This connection is trivial for generic 
values of a and the z^: the ipj are flat sections, and the connection on T-L is fully described 
by the smooth evolution of the classical critical points and corresponding thimbles. (In 
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transporting the thimbles, one must keep track of their orientations; the sign of the relative 
homology class associated to Ji depends on the orientation of J/.) 

Crucially, there are codimension one walls in the space 5^ x 55 of parameters q and Zi 
where the thimbles fail to be closed (and so do not define elements of the relative homology), 
and the map from the critical points to quantum states jumps discontinuously. This can 
occur if there are gradient flow lines from / to J; in this case, Ji is not closed as it contains 
points arbitrarily close to J, but not J itself. 

We write i\ for the locus in S"^ x 5S on which the following necessary conditions are 
obeyed for flows from / to J: the value of the conserved quantity j is equal at / and J, 
while > /i(J). The first condition is a single real condition, while the second is just 
an inequality. So Ij is of real codimension 1, and we call it a Stokes wall. 

In crossing a Stokes wall Ij, only the thimble Ji becomes ill-defined. It jumps by a 
multiple of Jj: 

Ji^ Ji + mijJj, (4.9) 

where m/j is the "number" of gradient fiow lines from I to J counted in an appropriate 
sense. A given line contributes 1 or —1 to the sign depending on the direction in which the 
difference between the values of j = lm{e^'^yV{w, z)) at / and J passes through zero. For 
an elementary explanation of such matters, see section 2 of [40]. 

The correspondence between states ^7 and thimbles J^j means that the V'/ have the 
same jumping in crossing Stokes walls. In the ipj basis, the connection is trivial except 
across the Stokes walls, where the transport matrix is a triangular "Stokes factor" 

S[ii] = l+mije-^i (4.10) 

Here e"^/ is the matrix whose only non-zero element is 1 at position J, /. 

In particular, the parallel transport along a path P in 5^ x *B is a path-ordered product 
of factors of the following kind: (a) between two Stokes walls, one has only the "formal 
monodromy" which expresses the permutations of the classical critical points, with minus 
signs that keep track of the orientations of the thimbles; (b) every time one crosses a Stokes 
wall, the monodromy acquires a corresponding Stokes factor. 

One can visualize the matrix elements nj[V] of the transport matrix for the path V 
as counting paths from a critical point / to a critical point J, where away from Stokes 
walls, one has to follow a critical point continuously, but in crossing a wall, one is allowed 
to "jump" along a gradient fiow trajectory from one critical point to the next. A matrix 
element of the transport matrix is given by a sum of contributions of hybrid paths of this 
type, with each path contributing 1 or —1 depending on how the orientation of a thimble 
evolves along the given path. 

Incidentally, the fact that the monodromy for parallel transport along a path V depends 
only on the homotopy class of V implies wall-crossing formulas for the numbers m/j that 
control the jumping. More generally, the whole picture can be interpreted in terms of BPS 
states in an LG model based on W, but we will not pursue that interpretation in this 
paper. 

Let us collect a few properties of the matrix elements nj[V]: 
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• [P] only depends on the homotopy class of V; 

• The composition of paths maps to matrix multiplication: nf\PioV2] = Xlj t^/ ['Pi]tIj'['P2]; 

• For an infinitesimal path 5V from z to z + 5z ii is almost always true that nj[5V] = 5j 
- this fails only in crossing a Stokes wall; 

• In crossing a Stokes wall, nj[5V] — 6j = m/j, where m/j is computed by a count of 
flow lines. 

4.1.3 The Dual Basis 

If r is any cycle in the relative homology H{U,U<:), then as the thimbles are a basis for this 
relative homology, T is equivalent in relative homology to a linear combination of thimbles, 

r = ^c^j/. (4.11) 
I 

How can the coefficients be determined? If we had, in some sense, a dual basis 
the j7/'s with pairings 

then we would identify the coefficients in (4.11) as 

c^ = (/C^r). (4.13) 

The Poincare dual of the relative homology H(U,U^) is the opposite relative homology 
H{U,U^) of U relative to the region Uy with h — >• +oo. A natural basis of H{U,Uy) is 
given by the upward-flowing thimble K,^ associated to critical points. For each critical 
point /, K,^ is defined as the boundary values at t = of solutions of the gradient flow 
equation on the half-line [0,oo) that approach the point / for t — )■ -|-oo. Since h decreases 
along gradient flow lines, the smallest value it assumes along such a flow is its value at 
t = oo, which is its value at the critical point /. So h is bounded below along /C^, but 
possibly not bounded above, and K,^ takes values in H{U,U-^). As for the pairing (4.12), 
from the definitions of the /C's and J^s, {K,^ , Ji) counts flows on the full real line (— oo, oo) 
that start at J at t = — oo and end at / at t = -|-oo. For J ^ I there are no such flows 
(for generic a where the thimbles are well-defined). For J = I, since h strictly decreases 
along a nonconstant flow, the only flow is the constant one that sits at J for all times. For 
a certain natural relative orientation of the ^T's and /C's, the contribution of the constant 
flow to the pairing is -|-1. 

4.1.4 Non-Single Valued Superpotentials and (/-Grading 

In the framework that we have presented so far, the matrices which represent the action 
of the braid group on Ti have integer-valued entries, since the relative homology has an 
integral structure. Now we want to consider a situation where W is not single- valued. To 
be more precise, we will consider a holomorphic function like the Chern-Simons functional 
W, whose real part is single- valued, but whose imaginary part is well defined only modulo 
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27rZ. We can reduce to the framework which we have employed so far by replacing U by 
the smallest cover lA on which W is single- valued. 

Passing to lA comes at the cost that now the number of critical points will be infinite, 
since each critical point in U has infinitely many preimages vnU. A critical point / in U is 
the same as a critical point I in U together with a choice of a branch of W at /. Locally, 
we can pick an arbitrary preimage Iq of a critical point /, and denote as In the critical 
point for which the value of W is shifted by — 27rin compared to the value at Iq. 

The Stokes factors and transport matrices are now matrices of infinite size, but their 
matrix elements can be computed by the techniques we have described, and are integer- 
valued. As the calculations only depend on the gradient of W, these matrices commute with 
the deck transformation /„ — )• /„+i. It is convenient to introduce a variable q taking values 
in C*. Then if we simply write q'^ipi as a symbolic shorthand for ih-j- , we can replace infinite- 
dimensional matrices whose entries are integers with finite-dimensional matrices whose 
entries are Laurent polynomials in q with integer coefficients. We denote such a matrix as 
S[£j; q]. In order to keep track of the lift of W as we move in 03, the formal monodromy 
matrices which encode the permutation of critical points and the Stokes matrices associated 
to gradient flows are valued in powers of q, to keep track of the change of W along a path. 
The matrix elements of the transport matrix are now obtained by summing over hybrid 
paths (continuous evolution away from Stokes walls and gradient flow across Stokes walls) 
with weight =tq", where — 27rm is the change of superpotential along the path and as usual 
the sign involves the orientation of the thimbles. 

An alternative description of all this is as follows. Since lA admits the non-single- valued 
superpotential W, its first Betti number is positive and one can introduce a "theta-angle" 
6 in the super symmetric quantum mechanics, weighting by e*"^ a path in which W jumps 
by — 27rin. This has the effect of replacing the relative homology of U with a twisted 
version of the relative homology, valued in a flat line bundle of monodromy q = e^^. Then 
we consider the Gauss-Manin connection for homology twisted by this flat bundle. The 
holonomy matrices for this connection have entries that are Laurent polynomials in q with 
integer coefficients and they can be computed as just described. 

4.2 Classical Description of Counting of Four-Dimensional Solutions 

Now, let us consider the problem that we are really interested in - time-dependent solutions 
of the super symmetric equations (1.1) on R x C x M+. The solutions will be time-dependent 
because the boundary conditions are time-dependent - we allow the positions Zi of singular 
monopoles on the boundary of M3 = C x M4. to vary with time. Although we typically 
assume an adiabatic evolution of the monopole positions, the counting of four-dimensional 
solutions is topological, and the adiabatic assumption is not necessary. 

In the simplest setup, with singular monopole strands at the boundary braided in 
time, the superpotential depends holomorphically on some parameters (the positions of 
the strands in C), and the parameters evolve in time. Schematically, the equations take 
the form of a "forced gradient fiow" 



dt 



dwj 



(4.14) 
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We suppose that the singular monopoles begin at positions Zi = {zi, . . . , z^) near time 
t = — oo and end at positions Zf = {z[, . . . , z'^.) near t = +00. In fact, we can assume that 
the positions Zi{t) of the singular monopoles are constant except in a bounded interval 
—T < t < T, for some T, during which they follow a path V in their parameter space 
5S. In such a situation, we can look for solutions of the forced gradient flow equation that 
begin at a specified critical point / of >V(w,ii), and end at a specified critical point J 
of yV{w,Zf). The "number" of such solutions, with each solution weighted by the sign of 
the fermion determinant, is a topological invariant - unchanged under deformations of the 
path V or the metric on U. We will call this invariant 91/ . (For a reason that will be clear 
momentarily, we really only want to define if Zi and Zf are not on Stokes walls.) 

In section 4.1.2, we already associated an integer invariant nj[P] to this situation. 
nj[V] was a matrix element of the Gauss-Manin connection for transport along the path 
V from Zi to zj. One can think of it/[7^] as the expansion coefficients when a thimble 
J'j in the relative homology H{U,hl<^^^ is transported along the path V using the Gauss- 
Manin connection, and then expressed in terms of the thimbles Jj that furnish a basis of 

Ji = ^njj!j. (4.15) 
J 

We claim that in fact 

mj = n/. (4.16) 

The importance of this relation is that 91/ is what we want, the counting of time-dependent 
solutions, while nj is more easily computed, since this requires only the study of time- 
independent problems. 

One explanation of (4.16) is as follows. By definition, J'i[zi] is the set of points in U 
which can be reached by flows 

d^ ^ ^i^ j^ dWjw, Zi) 
dt dwj 

which asymptote to / in the past. Here the Zi are regarded as constants. Now consider the 
equation (4.14) for forced gradient flow on the semi-infinite interval (—00, to]- For ^ —T, 
the values at to of a solution of this equation parametrize Jilzi], but for to > —T, they 
parametrize a to-dependent continuous deformation of this space that we will call Ji{zi;to\. 
We saw in section 4.1.3 that that for any given cycle F in Hm(piM<)zj^ the coefficients 
in the expansion 

V = Y,c'jj[zf] (4.18) 
J 

count, in the sense of an index, the number of fiows 

^ = -e-V^- ^^i"'"-^ (4.19) 
dt owj 

which start from F and asymptote to the critical point J in the future. 
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Hence the number of four-dimensional solutions which flow from Ji[zi]tQ\ at some given 
time to after the braiding occurs to the critical point J in the future are the coefficients 
9Tj in the expansion 

Ji[z^-M=Y,'^iJj[^f]- (4-20) 
J 

But since i7/[ii;to] parametrizes flows on the interval (—00, tg] that start at /, a flow from 
Ji[zi;tQ\ to J on the interval [tojCo) is equivalent to a flow from I to J defined on the 
whole real line. So the 91/ are the same as the desired invariants 91/ : 

Ji[zuto]=^^-}jj[zf] (4.21) 
J 

The relative homology is defined over Z, and an integral relative homology class such as 
i7/ has no continuous deformations. So clearly, as long as the continuous deformation 
from J'i[zi\ to Ji[z;to] induced by the fiow equations (4.14) lives at any given time in 
Hmipl M<)z(t)^ it coincides with the natural transport along V by the Gauss-Manin con- 
nection. In this case, (4.21) is equivalent to the desired result 9T/ = n/. 

To show that J7[z;to] lies in Hm(ptM<)z{to) for any t^, we are supposed to prove that 
ReyV(w;, z{to)) goes to —00 at infinity along Ji[to\. Indeed, if a sequence of forced gradient 
flows on the semi-infinite interval (— oo,to] goes to infinity, it does so by diverging for 
t — )• toi in which case Re(W(7y, -2(^0)) (whose gradient drives the fiow for t — >■ to) must go 
to —00. 

An alternative approach to (4.16) is the following. Suppose that Zf = Zi and V is the 
trivial path between them. Then 9t/ = nj = 6j. As we vary zj, both 9X/ and nj may jump 
in and only in crossing Stokes walls; they jump in exactly the same way, so they remain 
equal. We have already describing the jumping of nj. The jumping of 91/ occurs because 
in crossing a Stokes wall, a time-dependent solution may disappear to infinity, as follows. 
Suppose that, for Zf on some Stokes wall, there is a jump in n/, resulting from an ordinary 
gradient fiow from some critical point J' of }V{w, Zf) to J. Such a fiow produces a jump 

n/^n/±n/', (4.22) 

where the sign depends on the direction in which one crosses the Stokes wall. To see a 
corresponding jump in 9l/, one looks for forced gradient trajectories from I to J that 
consist of a forced gradient trajectory from / to J' followed, at some time very far in the 
future, by the same ordinary gradient trajectory from J' to J that causes the jump of n/. 
A two-step forced trajectory of this kind exists if Zf is near the Stokes wall and on the 
proper side of it; the time at which the second step of the fiow occurs diverges as Zf crosses 
the Stokes wall. This leads to the disappearance of the two-step solution and a jump of 
9T/ that just matches the jump of nj. 

5 Prom Braiding Of Thimbles To Free Field Integrals 

According to the reasoning in section 4, to understand the braid group representations 
associated to the Jones polynomial, we are supposed to compute a natural monodromy 
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action on the middle-dimensional relative homology of an infinite-dimensional space U of 
connections on a three-manifold M3 = C x M+. This may sound hopelessly abstract. We 
will now show how it can be turned into something concrete and calculable. 

5.1 From Thimbles to Integrals 

A convenient way to describe the evolution of states in including the g-grading, is to 
view the homology cycles as integration cycles. Instead of looking at the evolution of the 
homology cycles, it is equivalent to look at the evolution of the integrals. Of course, we 
need an integrand which can be integrated on the thimbles J7/, which are not compact. A 
function that fills the bill is e^^"^, where e is chosen in a suitable half-plane. Since 

/i = Re (e'"W) (5.1) 

goes to minus infinity along a thimble, the condition we want is Re (e*°e) > 0. This will 
ensure the convergence of the integrals 

Xv = j e^('")/MO (5.2) 

where dO is a holomorphic volume form on U (which will be kept fixed in what follows) and 
the integration cycle T is a thimble, or more generally any cycle in the relative homology 
H{UM<)- 

The thimbles are particularly nice integration cycles, because the e — t- limit of the 
integral over a thimble is very simple. On a thimble Ji defined by gradient flow from a 
critical point /, the function h has a unique maximum, namely the critical point. So for 
e — )• 0, the integral over a thimble is 

Xi := j e^^"")/" dO ~ exp(W//e) {e~'^'''^^'/^co + ...), (5.3) 

where W/ is the value of W at the critical point /. 

This formula is valid throughout the half-plane Re (e*°e) > 0, but actually as long as 
a is not on a Stokes wall, this asymptotics holds in a slightly larger sector in the complex 
plane. ^'^ This property uniquely characterized the basis of integrals X/ among all the 
possible Iy'- if we were to take a linear combination of several X/, the asymptotics would 
fail at some ray in the extended half-plane where two critical points exchange dominance. 
This characterization is familiar in Stokes theory, and motivated the terminology "Stokes 
walls." 

As we vary the parameters of W, the homology H{l/(,U^) will vary continuously. If 
we vary the integration cycle T continuously, the integral will vary holomorphically in the 

^°For the stated asymptotics to break down, the first step is to cross a Stokes wall, so that the thimble we 
started with evolves into a linear combination of thimbles with at least two terms. Initially, the asymptotics 
(5.3) remain valid, as any extra thimbles that appear at the Stokes wall initially make exponentially small 
contributions. If one varies a further, one of the extra thimbles may eventually become dominant. The 
combined process always involves varying a by an angle strictly greater than n/2 from its initial value. To 
show this, one just compares the values of W at the two critical points; these values have equal imaginary 
parts at the Stokes wall, and equal real parts when the two critical points exchange dominance. 
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parameters of W. The monodromy of the cycles T is the same as the monodromy of the 
integrals Ir- 

5.2 Two Chern-Simons Theories 

In our present context, the thimble integral (5.2) is a Chern-Simons path integral on the 
three manifold M3 = C x M+, albeit on an unusual integration cycle. Such an integral 
can be concretely expressed in terms of A/" = 4 gauge theory on M3 x . (In this section 
only, we write M or for the direction to distinguish it from the y direction M_|_.) This 
statement was one of the main conclusions of [40, 41]. In our case, since M3 = C x M_|_, 
the four manifold is M4 = C x IR+ x IR+. 

This is the second Chern-Simons theory to appear in this paper, at least implicitly. 
Our whole analysis concerns the calculation of the Jones polynomial, in a gauge theory 
setup which is S'-dual to a setup which computes the Jones polynomial by Chern-Simons 
theory on W = C x M. In that "original" Chern-Simons theory, the knot is a Wilson 
loop, the gauge group is G^, and the coupling parameter is k"^ . The Jones polynomial is a 
Laurent polynomial in 

g = exp(27rV(fc^ + 2)), (5.4) 

where 2 is the dual Coxeter number of SU{2). 

As explained in [14], Chern-Simons theory on a three-manifold W can be computed 
via topologically twisted = 4 super Yang-Mills theory on x if one relates to 
the twisting parameter of the AA = 4 theory by 

^'^ = ±{k^ + 2). (5.5) 

In this description one uses a D3-NS5 boundary condition at the origin in M+. The sign it 
depends on the relative choice of orientation between M4 and W . 

One can also apply S-duality, converting the gauge group from G"^ = SU (2) to G = 
5*0(3) and converting the D3-NS5 boundary condition to a D3-D5 boundary condition; this 
boundary condition involves a Nahm pole, as explained in [19]. In this new description, 
which has been the starting point of the present paper, the twisting parameter is ^' = 
— 1/^^. Because this dual description is difficult, we have tried to simplify it, as first 
explained in section 1.2, by "stretching" W in one direction. Thus we approximated W by 
M X C, and looked for solutions of the BPS equations onWx M_|_ that are "pulled back" from 
M3 = CxM+. If this dual description is formulated as supersymmetric quantum mechanics, 
with M as the time direction and the field variables being gauge fields on M3 = C x M_|_, 
then the superpotential is ^ (not ^'^) times the Chern-Simons function. 

We then relate the braiding of solutions on M3 to the braiding of thimbles for Chern- 
Simons theory on M3, which we express by AA = 4 theory on M3 x M_|_. In this description, 
we impose D3-NS5 boundary conditions at the origin of ]R_|_ (and D3-D5 at the origin of 
M+). The level k of this Chern-Simons description is related to the twisting parameter of 
the A/" = 4 theory by 

^ = T{k + 2). (5.6) 
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The two Chern-Simons descriptions are related by S-duality, together with the exchange 
of the roles of M+ and M+. The opposite sign in (5.5) and (5.6) is due to the fact that 
exchanging M+ and IR+ reverses the orientation of M4. 

Combining this with ^ = —1/^"^ and '^^ = k'^ + 2, we find that the relation between 
the level A;^ in the Chern-Simons description that is related to the Jones polynomial in 
the traditional way and the level k in the Chern-Simons description that relates the Jones 
polynomial to Nahm poles and opers is 

k + 2 = l/{k'^ + 2). (5.7) 

This four-dimensional setup, with the A/" = 4 theory on a manifold M4 = C x M+ x M+ 
with a "corner," is rather interesting, and we believe it deserves to be explored further. 
We will not do so in this paper. 

5.3 Prom Chern-Simons to Conformal Blocks 

For this paper, more useful than the relation of Chern-Simons theory on M3 = C x to 
four dimensions is its relation to conformal blocks on C. The most familiar version of this 
statement [44] is that a path integral with Az fixed and Az varying gives a WZW conformal 
block (that is, a conformal block in two-dimensional current algebra, with the symmetry 
group G being the same as the gauge group of the Chern-Simons theory). Local variations 
of the fixed value of A- insert a holomorphic current J{z) in the conformal block: 

6logZ = J Tt{Jz)5A- (5.8) 

This statement has an analog [45] that leads to Virasoro conformal blocks in the case 
of SU{2) or 5*0(3) gauge theory, or to more general VF-algebra conformal blocks in the case 
of gauge groups of higher rank [46, 47]. This analog involves a different boundary condition 
in which, at the price of breaking some gauge symmetry at the boundary, one fixes some 
parts of Az, and some parts of A-z. For gauge group SU{2), the boundary condition which 
leads to Virasoro conformal blocks is simply stated: 





^^=11 (5.9) 



Here we denote as * the elements which are free to vary, and as x elements which are fixed. 
The connection just described is an oper! (Since A-z is lower-triangular, a bundle with 

this connection has a holomorphic sub-bundle L whose sections are of the form , and 

because the upper right matrix element of nowhere vanishes, this sub-bundle is nowhere 
preserved by T>z = dz + [Az, ■]■) So the complex boundary condition (5.9) is the one that 
is induced by the Nahm pole boundary condition studied in section 3. To explain the 
relation to Virasoro conformal blocks, we note that the Nahm pole boundary condition 
depends on a choice of complex structure on C. Once a complex structure is picked with 
a local complex coordinate z, nearby complex structures can be described by a Beltrami 
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differential Tlie relation is that in the new complex structure, the holomorphic fields 
on the space of complex connections are not As but Az — fJ^Az- Making this deformation 
is equivalent to replacing the "0" in the boundary condition for Az in (5.9) with fj^ (so 
that now A- — ^Az is lower triangular). Hence a local variation of the boundary condition 
associated to a change in complex structure inserts a holomorphic stress tensor 



and this leads to the relation between Chern-Simons theory with the oper boundary con- 
dition and Virasoro conformal blocks. 

In the semiclassical limit, the operator Tzz reduces to the classical stress tensor t(2;)[^] 
of the oper. More precisely, the identification of parameters from [45] is that if we define 



with k the Chern-Simons level, then the stress-tensor has central charge c = 1 -|- 6(6+ 6^-*^)^ 
and in the semiclassical limit, 6^T(z) — t- t(2;)[^]. 

If we combine (5.11) with (5.7) and (5.4), we find the relationship between the variable 
q usually used in describing the Jones polynomial and the parameter h used in describing 
Virasoro conformal blocks: 



We have here assumed that G is SU{2) or 50(3). For general gauge group G, both 
Nahm poles and M^-algebras are labeled by an 5u(2) embedding in the Lie algebra g of 
G. Inspection confirms that the boundary conditions used to define a general VF-algebra 
conformal block are induced by the corresponding Nahm pole. 

5.3.1 Analog For Liouville 

Though we will not need this fact in the present paper, we should remark that the relation 
between Virasoro conformal blocks and Chern-Simons theory has a simple extension to a 
relation beween Liouville theory and Chern-Simons theory. To do Liouville theory on a 
Riemann surface C, one considers Chern-Simons theory on C x / where I is a unit interval. 
At one end of /, one imposes the Nahm pole boundary condition and at the other end, one 
imposes a variant of the Nahm pole boundary condition with z and 1 exchanged. Liouville 
partition functions and correlation functions are built by combining holomorphic and anti- 
holomorphic Virasoro conformal blocks, which arise naturally from Chern-Simons on C x I 
with boundary conditions just stated. (For the case of a compact symmetry group, it is 
already known that Chern-Simons on C x / reproduces the WZW model on C.) By slightly 
extending arguments that we present presently, light degenerate fields and generic primary 
fields of Liouville theory, inserted at a point p £ C, correspond to Wilson operators or 
monodromy defects on p x I. 

In the classical limit, the correspondence between Chern-Simons and Liouville theory 
means the following. A classical solution of Chern-Simons theory on C x / with boundary 
conditions as above is a flat bundle on C whose holomorphic and antiholomorphic structures 




(5.10) 



- 6-2 = fc + 2 



(5.11) 



q = exp(— 27ri/6^). 



(5.12) 
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both obey the oper condition. Indeed, a classical solution of Liouville theory corresponds 
to a metric on C of constant negative curvature. If u and e are the vierbein and spin 
connection of this metric, then we can define a corresponding SL(2, M) fiat connection 
A = ujt^ + 6^4+ + Czt-. With a standard representation of the ti, Az is upper triangular 
with an upper right matrix element that is everywhere nonzero, and Az is lower triangular 
with a lower left matrix element that is everywhere nonzero. So both the holomorphic 
and antiholomorphic structures defined by this fiat connection satisfy oper conditions, as 
expected in the Chern-Simons description. (The antiholomorphic oper structure is defined 
with the roles of "upper triangular" and "lower triangular" matrices reversed.) So this 
gives the mapping between the two theories at the classical level. 

5.4 Wilson Line Operators 

To further understand the mapping from three-dimensional Chern-Simons theory with oper 
boundary conditions to Virasoro conformal blocks in two dimensions, we will explore the 
interpretation of Wilson line operators. First let us recall what happens if one uses standard 
boundary conditions that relate Chern-Simons theory to current algebra. In this case, a 
Wilson line operator ending on the boundary of a three-manifold M3 represents insertion 
of a conformal primary field at that boundary point in the WZW conformal block. If 
the Wilson line operator transforms in a finite-dimensional representation R of G, then 
the corresponding conformal primary field transforms in the same representation. This is 
consistent with the fact that a Wilson line operator ending on the boundary is not gauge- 
invariant, but transforms in the representation R, just like the corresponding primary field 
of the WZW model. 

What is the analogous interpretation of a Wilson line operator that ends on a boundary 
at which one imposes Nahm pole boundary conditions? The Nahm pole breaks the gauge 
symmetry at the boundary, so we have to pick a component of the Wilson line operator. 
As we discussed in section 3.2, a generic vector diverges when parallel transported to the 
boundary. Given a Wilson line operator ending at y = 0, the most easily defined gauge 
invariant information is the coefficient of the most negative power of y. This is extracted 
simply by contracting with an appropriate power of the small section s. Actually, we will 
find useful a rescaled version of s, namely s = K(z)~^^'^s{z), which satisfies 

{Vl + t{z))s = VzS = Vys = (5.13) 

and has definite conformal dimension —1/2. In the gauge (3.46), 's = 

In the two-dimensional representation of SU{2), we would consider an operator 

Pexp j A^ s{z), (5.14) 

where 7 is a path ending on the boundary at y = 0. In the classical limit, under conformal 
transformations of the boundary, this has the same conformal dimension as 's, i.e. —1/2. 
For a spin k/2 representation, one must contract with k powers of s~and the classical limit 
of the dimension is —k/2. 
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We want to argue now that a Chern-Simons path integral with a spin kjl Wilson 
operator ending on the boundary gives a Virasoro conformal block with the insertion of a 
"light" degenerate field of Liouville momentum —hk/2. Such a field has the correct classical 
dimension —k/2 in the & — )• limit, and furthermore the exact formula for its quantum 
dimension 



is the sum of the classical dimension of 's^ and the dimension of a spin k/2 operator in 
a WZW model (of level fc + 2 = —1/6^ as in (5.11)). Furthermore, the operator (5.14) 
satisfies classically the correct differential equation: the k = 1 operator is annihilated by 
+ t(z), etc. 

Part of what makes possible the correspondence between spin k/2 Wilson lines and 
degenerate conformal fields possible is that the degenerate fields satisfy fusion rules which 
coincide with the fusion rules of spin k/2 operators in the WZW model. This last fact is 
part of the input in the statement that the braid group representations associated to the 
Jones polynomial can be computed by the braiding of either primary fields of the WZW 
model or degenerate conformal fields of the Virasoro algebra. 

5.5 Singular Monopoles and "Heavy" Degenerate Fields 

Next, we would like to identify in the Chern-Simons description of Virasoro conformal 
blocks the "heavy" degenerate fields of Liouville momentum —k/2b. These are the degen- 
erate fields whose conformal dimension diverges for 6 — )■ 0. We claim that they correspond 
to the insertion of singular monopoles at the boundary. 

The main insight of section 3.3 was that at a Nahm boundary with singular monopoles 
(and generic C); the connection is an oper with singularities of trivial monodromy. We 
observed that the classical stress tensor of such an oper has poles that agree with the 
semiclassical limit of the quantum stress tensor in the presence of a heavy degenerate field. 
Moreover, the trivial monodromy condition holds quantum-mechanically as well: a light 
degenerate field has no monodromy around a heavy degenerate field, as they fuse in a 
unique channel. 

In the Chern-Simons setup, the classical trivial monodromy condition follows naturally 
from the fact that the singular monopole does not extend in the bulk. This was part of 
our derivation in section 3.3. Quantum mechanically, we need to consider the behavior 
when a light degenerate field - represented in three dimensions by an expression such as 
P exp ^— ?- approaches the singular monopole. The Wilson loop itself is topological, 

and the small section s' has trivial monodromy around the singular monopole. 

5.6 Putting The Pieces Together 

We can now finally establish a link between the solutions of the four-dimensional BPS 
equations (1.1) that we started with and the braid group representations associated to the 
Jones polynomial. 

The time-independent solutions of the BPS equations correspond to opers with trivial 
monodromy. We have identified the braiding of the corresponding quantum states with the 




(5.15) 
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braiding of complex integration cycles for Chern-Simons theory, and then with the braiding 
of degenerate Virasoro conformal blocks. These are known [34, 35] to be the braid group 
representations associated to Jones polynomials. So we have arrived at our goal, though 
in a form that may sound a little abstract. 

We can put this result in a perfectly concrete form using the free field representation of 
the conformal blocks. Opers with trivial monodromy are also associated to critical points 
of the Yang- Yang function for the Bethe equations; this is the logarithm of the integrand 
in the free field realization of conformal blocks. We can derive a degenerate Virasoro 
conformal block either from an infinite-dimensional thimble associated to an oper with 
trivial monodromy, or from a finite-dimensional thimble associated to a critical point of 
the Yang- Yang function. Either way, we get a conformal block with definite and uniform 
semiclassical limit in a sector of angular width greater than ir in the 6^ plane. As those are 
unique, the two bases of conformal blocks must coincide. 

Hence the braiding representations associated to the four-dimensional gauge theory 
coincide with the braiding representations of integration cycles in the space of Bethe pa- 
rameters Wi. This is not as surprising as it may seem if we turn on a symmetry breaking 
parameter: then we have interpreted the Wi as positions of bulk BPS monopoles, and our 
claim possibly amounts to the statement that the four-dimensional nonabelian gauge the- 
ory on the Coulomb branch reduces to a theory of massive monopoles and abelian gauge 
fields. We will develop this point of view further in section 7. 

6 Braiding Representations of Integration Cycles 
6.1 Overview 

A highlight of what we have learned so far is the existence of a natural map from the braid 
group representations derived from the four-dimensional gauge theory equations (1.1) to 
the braid group representations associated to correlation functions of Virasoro degenerate 
fields. Those braid group representations can be effectively studied using the free field 
representation, which we reviewed in section 3.7. Making this explicit will be our goal 
here. 



We consider a degenerate correlation function (^00=1 ^-ka/2b(.Za)^- We assume that 
the Za are distinct points in C. To represent conformal blocks, we introduce q variables 
Wi £ C, which we assume to be distinct from each other and from the Za- The allowed values 
of q have been analyzed in section 3. Moreover, we consider the Wi to be indistinguishable, 
in the sense that configurations that differ by permuting them are equivalent. We write Ai 
for the space of such distinct and indistinguishable variables Wi € C\{zi, . . . , z^}- We also 
write M for the smallest cover of M on which the Yang- Yang function W of eqn. (3.52) 
is single- valued. 

In the free field representation, degenerate conformal blocks are written in the form 
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r is a middle-dimensional cycle in M, chosen so that the integral converges. 

Morse theory offers a systematic way to produce all such integration cycles: a basis 
of integration cycles is given by the thimbles associated to critical points of W. Cycles 
of this kind are never compact; they have noncompact ends on which the Morse function 
h = Re W goes to — oo. In our problem, this happens when one of the Wi either approaches 
one of the Za or, in the presence of symmetry breaking, goes to infinity in the correct 
direction. In simple situations, instead of using Morse theory, one can describe integration 
cycles by hand. In constructing a cycle T by hand, one has to make sure that the Morse 
function really goes to — oo at infinity along T. For example, this will fail if too many Wi 
approach simultaneously the same Za- 

With symmetry breaking, some of the important integration cycles have ends at = oo 
and the use of noncompact integration cycles is unavoidable. However, in the absence of 
symmetry breaking, the noncompact integration cycles produced by Morse theory have 
their ends at Wi — )• Za, for various i and a, and are equivalent in the appropriate twisted 
relative homology to compact cycles in which the Wi wrap around the Za in a suitable 
fashion. (For an example, see fig. 7 below.) In the extensive literature on integration 
cycles in free- field realizations of conformal blocks [32, 33] and their application to the 
Jones polynomial [34, 35] , compact integration cycles are often used. Symmetry breaking, 
or in other words the introduction of an irregular singularity at infinity, has not been 
considered in this context, as far as we know. 

The use of Morse theory has advantages and disadvantages. The main disadvantage is 
that the thimbles do not correspond to a standard BPZ basis of conformal blocks defined 
by fusing the degenerate fields in specific channels. The main advantage is that in the 
basis of thimbles, the braid group is manifestly represented by matrices whose entries are 
Laurent polynomials in q with integer coefficients. This property, which was explained in 
section 4.1.4, is important vis-a-vis the Jones polynomial and Khovanov homology. 

In what follows, we will first analyze a few important examples with a small number 
of degenerate insertions Za and Bethe roots Wi, with or without symmetry breaking. Then 
in section 6.5, we analyze the general case in the presence of symmetry breaking. From 
that analysis, we get the experience we need to deduce a general description of the Jones 
polynomial - not just the associated braid group representations. This is presented in 
section 6.7. 

6.1.1 SU{2) Versus S0{3) 

We pause for a technical remark concerning the assertion that the entries of the braiding 
matrices are Laurent polynomials in q. 

^^As explained in [34], the cycle F should actually be odd under the exchange of any pair of w's, to 
compensate for the sign change of the differential form dwi A ■ ■ • A dwq under permutations. This means 
that the appropriate relative homology is actually the part that is antisymmetric under permutations of the 
w's. To be concrete, suppose that there are two w's and we find a solution of the Bethe equations at which 
the w's equal a and /3 up to permutation. Then we can define a cycle C' associated to the critical point 
wi = a, W2 — /3, and a cycle C" associated to the critical point wi = (3, W2 = a. The difference C' — C" 
is an element of the antisymmetric part of the homology. In practice, we can omit to explicitly form such 
differences and also ignore minus signs arising from permutations of the factors in dwi A • • • A dw„. 
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If the gauge group is G = SU{2), meaning that the dual group is = 50(3) and 
the charges ka of the singular monopoles are all even, then the Yang- Yang function W 
as defined in (3.52) is well-defined mod 27ri. Hence a change of branch of W multiplies 
exp(W/6^) by an integer power oiq = exjp{—2iri/b^), and the braiding matrices are Laurent 
polynomials in q with integer coefficients. 

If instead G = 50(3), = SU{2), then some of the ka may be odd. (Indeed, we 
will do our detailed computations for the case that all ka are 1.) Then W is well-defined 
mod 2111 if the Wi are varied for fixed positions Za of the knots, but is only well-defined 
mod 27ri/4 when the Za are varied. Consequently, for G = 50(3), the braid matrices will 
actually be Laurent polynomials in q^^^. 

Of course, we could eliminate this by writing the formulas in terms of q = q^^^, but 
we prefer not to do so since q as we define it is the natural instanton counting parameter 
in four dimensions. The underlying reason for the difference between G = SU(2) and 
G = 50(3) is that the Chern-Simons function, normalized as we have done in (3.6), is 
gauge-invariant mod 27ri in SU{2) gauge theory, but gauge-invariant mod 27ri/4 in 50(3) 
gauge theory. The last statement holds on any sufficiently rich three-manifold M3, or on 
any M3 if singular monopoles of odd charge are present. In the latter case, one considers 
only gauge transformations that are trivial at the position of the singular monopole. 

6.2 A Single Critical Point 

We begin with the two examples in which the Yang- Yang function has only a single critical 
point. 

6.2.1 Braiding of Two Primaries With Minimal Charge 

The first example arises in the absence of symmetry breaking, with two Za of charge k = 1 
and one w. An obvious integration cycle C12 is a segment joining the two Za- Now, let us 
compare it with the thimble. The Bethe equation 

+ = (6.2) 

W — Zl W — Z2 

has a unique solution w = ^{zi + Z2)- The thimble flows down from w to the Zi along a 
straight line, and coincides with the obvious cycle Ci2- The critical value of the Yang- Yang 
function 

1 3 
W = -- log(zi - Z2) + log(tt; - Zl) + log{w - Z2) = - \-Og{zi - Z2) + const. (6.3) 

is such that exp(W/6^) agrees with the OPE coefficient of the two degenerate fields in the 

3 

identity channel V^i/2b{zi)V^i/2b{z2) ~ {zi — Z2y^. (Recall that for small 6, the dimension 
of y_i/2fe is —3/46^ + 0(1).) This is expected, because with two Za and one w the oper has 
no singularity at infinity, in view of the discussion of eqn. (3.22), and hence describes the 
fusing of two fields to the identity. 

In order to fuse the two degenerate fields in a channel of momentum —b we would 
consider something even simpler: a case without any w's. In this channel, the conformal 
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Figure 7. The thimble C12 (top) compared to a closed integration contour (bottom), which is 
equivalent to (1 — q^^)Ci2. 



Z W 

m • 



Figure 8. This ray parallel to the real axis is the Lefschetz thimble for the case of one primary 
field of minimal charge with symmetry breaking. 

block is simply the free field correlation function {V^i/2b{zi)V_i/2b{z2)) with no integral 
at all. With no w''s, the discussion of (3.22) shows that the oper does have a singularity 
at infinity - corresponding to fusion of the two k = 1 degenerate primary fields to a fc = 2 
degenerate primary field. 

To be precise, in defining an integration cycle such as C12, we should specify a choice 
of branch of W. (Differently put, the cycle is supposed to be defined in the covering space 
M..) If we braid zi around Z2, Cu evolves continuously, but we may end up with a different 
branch of W. As the charges at zi and Z2 are identical, a basic braiding move is to exchange 

the position of zi and Z2, and then relabel them. In the absence of w, we would only have 

]_ 

the factor {zi — Z2) coming from the part of W which only depends on z^, so a braiding 

_ 1 

which exchanges zi and Z2 counterclockwise would give a factor of e SP' = gi. In the 
presence of one w, we can get the result by following the value of W at the midpoint of 
C12: 

C12 ^ -q-'^Ci2. (6.4) 

The minus sign follows from the change in orientation of Ci2- 

It is interesting to compare the thimble with closed contours which are commonly 
used in order to describe BPZ conformal blocks, as depicted in fig. 7. As illustrated in 
the picture, these contours are equivalent in homology to the thimble times a Laurent 
polynomial in q. From the point of view of the Stokes matrices, they are not as elementary 
as the thimble. 
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Figure 9. The cycles Ci, C2 and C12 
6.2.2 One Primary Field With Symmetry Breaking 

The other basic example with one critical point occurs in the presence of symmetry breaking 
with a single degenerate field of /c = 1, and a single Bethe root. The Bethe equation reads 

1 2a , , 

(6.5) 



w — z C 

So w = z + K For convenience, we will take the constant c to be real and positive. 
The thimble again coincides (fig. 8) with the most natural integration cycle C, along a 
ray starting at z and parallel to the positive real axis, passing through w. This example 
illustrates that an integration cycle in the presence of symmetry breaking may end at 
infinity. The cycle that we have just described cannot be replaced with an equivalent 
compact cycle. 

6.3 Two Critical Points 

Next we can consider examples with only one Bethe root w, but two critical points. This 
occurs with two primaries of A; = 1 in the presence of symmetry breaking, or with three 
such primaries in the absence of symmetry breaking. 

6.3.1 The First Example 

We consider first the case of two primaries with symmetry breaking. We denote the posi- 
tions of the primaries as zi and Z2 and we continue to assume that the symmetry breaking 
parameter c is positive. As long as zi — Z2 is not real, there is a symmetric choice of basic 
integration cycles Ca, a = 1,2: rays which start at Za and are parallel to the positive real 
axis. The difference Ci — C2 = C12 is a segment from zi to Z2. See fig. 9. Any two of 
these three cycles can be thimbles, depending on the relative values of c and zi — Z2- Since 
there are always only two thimbles, it is never the case that all three of Ci, C2, and C12 are 
thimbles. 

In this simple example, we already have the basic ingredients of the general braid 
group representation. The elementary move is to exchange the two Za, either clockwise or 
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Figure 10. The cycles Ci, C2 and C12 after a braiding operation. 



counterclockwise. We start with a configuration in which the Za have distinct imaginary 
parts, so that the Ca are well-defined. We may as well take the real parts of the Za to be 
zero. We want to define the branches of W along Ci, C2, and C12 so that it is true that 
C12 = Ci — C2 . Picking any branch of W on C12 , we define W on Ci and on C2 so that at the 
unique point where Ci intersects C12 or where C2 intersects C12, the definitions agree. This 
will ensure that C12 = C\ — €2- Although Ci and C2 do not intersect, we can deform them 
slightly so that they meet at a reference point far to the right, and then the two values of 
W will agree at this reference point. (This is ensured by the fact that one can define W to 
be single- valued in the semi-infinite rectangle bounded by Ci, C2, and €12-) 

Suppose that Im (zi — Z2) > 0. Then we can cross to Im (zi — Z2) < in two ways, 
with zi passing either to the left or to the right of Z2- The two operations are inverses, 
so it will suffice to consider one in detail. If zi passes to the right of Z2, then Ci evolves 
continuously as a ray parallel to the positive real axis. On the other hand, C2 does not. 
We can use instead C12, defined continuously as a segment from zi to Z2- In the basis of 
Ci and C12, this braiding move is diagonal: we only have to keep track of the branches of 
W and orientation of cycles. We want to express the final result in the basis of Ci and C2. 

In the half-braiding, Ci is multiplied by a factor of q~^: z\ is transported clockwise 
around Z2 and the reference point does not move significantly. After the half-braiding, we 
rename C\ as C2: 



On the other hand, C2 becomes the cycle in fig. 10. It can be deformed to the sum 
of three pieces that zig-zag to and from Rez = 00, as in fig. 11. Each of the three 
pieces is equivalent to a power of q times one of the original cycles C\ and C2. Just as in 



pieces are images oi —q 4C2 and q -iCi under a deck transformation. Hence the braiding 




(6.6) 



(6.6), one piece (after again exchanging the labels of Ci and C2) is q ^€2- The other two 
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Figure 11. The cycle C2, after the braiding operation, has been copied from fig. 10: it starts at 
^2, curves below zi to the left and then goes to Re z = 00. It is equivalent in homology to a zig-zag 
cycle, labeled C2 in the figure, which, starting at Z2 , heads directly to Re z = 00 before doubling 
back around zi and returning to Re z = 00. Thus C2 is the sum of three pieces, each of which heads 
to or from Rez = 00; each piece is equivalent to a power of q times an elementary cycle Ci or C2 (a 
ray starting at zi or Z2 and parallel to the positive z axis). 

transformation is 

C2^q~^C2-qiC2+q~^Ci. (6.7) 
Notice that with these transformation rules, 

C12 ^ -q-^Ci2- (6.8) 

This is the same result that we found in the eqn. (6.4) (the sign of the exponent is reversed 
because in deriving (6.8), we braided zi clockwise around Z2); symmetry breaking does not 
affect the fact that C12 represents the conformal block in which the two degenerate fields 

3 1 

fuse to the identity. The braiding matrix has eigenvalues —q* and q~*, which correspond 
to the two possible fusion channels. The linear combination Ci + q^^C2 transforms as 

Ci + q-^C2 ^ q~HCi + g-^Ca). (6.9) 

and hence it represents the fusion in the channel of momentum —b. 

Now, other choices of normalization of Ci and C2 may occur more naturally in various 
situations. If we change the relative normalization between Ci and C2, say by setting 
Ci = q~^l'^C\ and C2 = q^/'^C2, we can get braiding formulae 

B12: Ci^q-'--iC2 (g"^ -g3)C2 + g"+3Ci. (6.10) 

We will find the choice s = — | to be useful momentarily, so whenever we write Ca we 
assume that choice of s. 

The behavior of the Ca as a function of c{zi — Z2) is depicted in fig. 12. The Ca fail to 
be well-defined when c{zi — Z2) is real. In this very simple example, the integral over the 
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Figure 12. The cycles (Ci,C2) (defined as rays in the direction of Recz that start at z = zi or Z2) 
provide a local basis for homology, but the definition of the cycles jumps by the braiding matrix 
B12 across the line on which Z ~ c{zi — Z2) is real. 

Ca can be expressed explicitly in terms of familiar functions: 

/ 1 1 _ 1 cz-i CZ2 

la = f {w — Zi)b^ {w - Z2)'^ {Zl - Z2) ^e~^^ ^'^dw 



Ca 

-^-^ i iW - Z)^W^ Z^'^e^-^dW (6.11) 

Ca 



2b^ 



Here we defined W = c{w — Z2) and Z = c{zi—Z2). This integral can be explicitly written in 
terms of Bessel functions. The integrals over Ci and C2 give a basis of Bessel functions with 
specific asymptotic behavior at large Z, and the braiding matrix B captures the Stokes 
phenomena of Bessel functions: 

1 1 +1 / 1 \ f Z 



r- 1 



z^+^rfi + lWi^, f-l) (6.12) 



In the large Z limit, the integral along Ca is controlled by the region where w is close 
to Za- It is useful to pick a branch of the logarithms such that for w near zi (or Z2), the 
sum of the logarithms in the superpotential approaches ^ log(2;i — ^2)^. This is the same as 
the change in normalization between Ca and Ca for s = — |. This basis of Ca will be useful 
whenever we are at strong symmetry breaking. In this basis 

B12 ■■= Ci ^ q-iC2 C2 ^ (^q~~* - q^^ C2 + q-*Ci. (6.13) 

So far, we have analyzed this problem using cycles Ci, C2, and C12 that are visible by 
inspection. Let us compare these cycles with the thimbles. The Bethe equation with two 
z's and one is ^ ^ 

+ = c (6.14) 



W — Zl W — Z2 

and it has two solutions. 

There are two regimes of interest. If \Z\ >> 1, then the critical points are approxi- 
mately w = Za + l/c, a = 1,2. For each of the two critical points, assuming that c > and 
Re Z = 0, both w — zi and w — Z2 have positive real part, and it is natural in defining W 
to pick branches of log{w — zi) and log(w — 22) such that the imaginary parts are bounded 
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Figure 13. The pattern of Stokes walls and the bases of thhubles in the c(zi — Z2) plane 



Figure 14. The two Morse flows from the critical point (empty dot) near Z2 to the critical point 
(empty dot) near zi. The flows occur at slightly different values of ImZ. 

by ib7r/2. The thimbles defined this way coincide with the Ca we defined above (and not 
with the Ca)- The advantage of this choice is that it extends naturally to the general case 
of many fields with symmetry breaking, which we will treat in section 6.5. If \Z\ « 1, 
then we have approximate critical points w = ^{zi + Z2) and w = 2/c. The first critical 
point sits between the two Za, and the associated thimble is C12. The second critical point 
is associated to Ci if ReZ > 0, C2 otherwise. 

We depict the Stokes walls for the system of thimbles in fig. 13. There are regions 
in parameter space where the basis of thimbles consists of any two of Ci, C2 and C12, up 
to a choice of branches of the superpotential. Pairs of regions meet along Stokes walls, 
and triples of regions meet at the points Z = ±2? where the two critical points coincide. 
This is a general feature; in a generic problem of this type, there will always be loci of 
complex codimension 1 where two critical points coincide. In a plane transverse to such a 
locus, yy can be modeled by a simple cubic function W = + Sw of one variable w, with 
a parameter 6; this function has two critical points that coincide for 5 = 0, where three 
Stokes walls meet. The local Stokes behavior is universal; it corresponds to the behavior 
of the Airy function. 

As we explained in our general discussion of section 5.1, the asymptotic behavior of 
integrals for 6 — >• is clearest in the basis of thimbles. For the present problem, fig. 13 
captures the relevant information. There are four regions: the upper and lower regions 
correspond to the system of thimbles we saw in the \Z\ >> 1 limit, the two intermediate 
regions to the system of thimbles which we saw in the \Z\ « 1 limit. 

The formulas (6.10) or (6.13) describe the braiding matrix B12 that compares the 
region at the top of fig. 13, where the thimbles are Ci and C2, to the region at the bottom, 
where again the thimbles are Ci and €2- We did this computation by inspection, not by 
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analyzing the Stokes lines. In deriving (6.13), we assumed that zi moves half-way around 
Z2 in a clockwise direction; this is equivalent to moving from the top to the bottom of fig. 
13 with ReZ > 0. From fig. 13, we see that in this process, we will cross two Stokes lines. 
This means that in the basis of thimbles, the braiding matrix B12 "decomposes" into a 
sequence of two elementary moves, each associated to one Stokes line. Each elementary 
move involves a gradient flow in which the Bethe root w flows from a critical point just to 
the right of Z2 to a critical point just to the right of zi. There are two such flows, differing 
by whether w passes above or below zi. The two flows, which occur at slightly different 
values of Im Z, are sketched in flg. 14. In the formula B12C2 = {q"^'^ - q^/^)C2 + q^/^Ci 
of eqn. (6.13), the term q^/'^Ci on the right hand side is the contribution of the formal 
monodromy alone, while the other two terms are contributions from the two gradient flow 
solutions. 

In the small Z region, braiding around Z = is represented in the basis of thimbles 
by a triangular matrix: C12 is an eigenvector. There is a general reason for this fact; the 
integral over C12 is smaller than the integral over any other cycle either in the Z — t- 
limit, or in the 6—7-0 limit. In either limit, the integral over C12 is controlled by the usual 
saddle point approximation, and the critical point associated to C12 is the one at which the 
Morse function is the smallest. So C12 must be an eigenvector of the monodromy. With 
any number of z's and u)'s, the set of thimbles which has one w in between a given pair 
of very close Za span the "small" eigenspace of conformal blocks where the two degenerate 
flelds of momentum — ^ fuse to the identity. 

6.3.2 The Second Example 

The second example with two critical points arises if there are three singular monopoles of 
minimum charge and a single Bethe root w. Placing the singular monopoles at zi,Z2,Z2,, 
there are three obvious possible integration cycles: a straight path connecting Za to Zh for 
any a, h. We can pick the branch of the superpotential in such a way that these three cycles 
add to zero. 

This definition makes sense if the three points are not aligned. The parameter space of 
Za is then split into two halves: either zi, Z2 and 2:3 form a triangle with positive orientation, 
or they form a triangle with negative orientation. We will denote the three natural cycles 
in either case as C^. So 

C+ + C2+3 + C3+ = Cf2 + C23 + C3I = (6.15) 

but the two bases are related in an interesting way across the loci where the Za are collinear. 

There are three such loci, where one of the three Za passes between the other two. For 
example, if Z2 passes between zi and Z3, CJ*^ and will be related to Cj^ and simply 
by a change of branch of W, while the transformation of Cjtj then follows from (6.15). 

The Bethe equation 

111 , , 
+ + = 6.16 

W — Zl W — Z2 W — Z3 

is equivalent to a quadratic equation for w, so it has two solutions, corresponding to two 
independent thimbles. The thimbles are equivalent to two of the paths joining a pair of z's. 
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Figure 15. The pattern of Stokes walls in the space of shapes of the triangle with vertices zi, Z2, 

but which pairs appear depends on the choice of the z's. If the z's are colhnear, then the 
Bethe roots are located in the segments between adjacent z's, and the thimbles coincide 
with those segments. For example, if Z2 is between zi and 2:3, the two thimbles are and 
C^, up to powers of q. 

As usual, the thimble joining Za and zj, corresponds to the conformal block where the 
corresponding two degenerate fields fuse to the identity. The relations (6.15) correspond 
to the elementary "skein relation" between the three different ways to fuse two of the three 
primary fields to the identity. In other words, whenever Za and z^ are close together, there 
is a thimble which joins them and is an eigenvector of the braiding of Za and Zb, and that 
braiding matrix is triangular. 

We depict the Stokes walls in fig. 15. The integrals lab can be readily evaluated in 
terms of hyper geometric functions. 

6.4 A Final Example 

There is one more example that is both instructive and relevant to understanding the 
general picture. This is the case of two Za of charge 1, accompanied by two Wi. This can 
only happen in the presence of complex symmetry breaking. It is easy to see that the 
Bethe equations (3.37) have only a unique solution. Q and K are both of degree 2, so P 
is of degree and can be set to 1. This leads to linear equations that uniquely determine 
the coefficients in Q. Since the solution of the Bethe equations is unique, there is only one 
conformal block and the monodromy in braiding the Za can only be multiplication by a 
function of q. 

If zi, and Z2 are well-separated, i.e. c{zi — Z2) has large absolute value, then the 
solution of the Bethe equations is easily described: wi ^ zi + l/c and ^2 ~ -22 + l/c There 
is an obvious integration cycle, with wi and W2 integrated respectively over the rays from 
zi and Z2 to infinity in the c direction; these rays were labeled Ci and C2 in fig. 9. We can 
denote this cycle as C = Ci x C2. 
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The unique thimble for the problem is a small deformation of C if the imaginary part 
oi Z = c{zi — Z2) is large. As we deform Z, the unique thimble will deform continuously. 
Here Morse theory is rather useful: it is rather tricky to verify by hand that C goes back to 
itself under braiding of zi and Z2 , as it requires a contour deformation which does not keep 
C in a simple product form. But the evolution of the thimble provides us implicitly with 
such a deformation. It is likewise far from obvious at first sight that the conformal block 
produced by the contour integral over wi and W2 will be a simple function with abelian 
monodromy in the Z-plane. Nevertheless, this must be the case. 

The effect of braiding is easily computed at large c, by looking at the saddle point 
estimate for the integral. At the saddle, tfi ~ zi + 1/c and W2 ~ 2:2 + 1/0 remains true 
along the whole braiding, as long as it is executed at large \Z\. The Hessian of W is close 
to a large, Z-independent, multiple (~ |cp) of the identity matrix, so the braiding phase at 
large Z is controlled by the value of exp(W/6^) at the critical point. The relevant factors 
are 

{Zl - Z2)-'~"/\wi - Z2f~\w2 - Zi)'-\wi - W2y~" ~ (zi " ^2)"'"'/^ (6.17) 

Hence we recover the same braiding phase as in a setup with two Za of charge 1 and no Wi. 

This last statement will be part of the input for constructing an effective abelian 
description in section 6.7. In an abelian theory, one would expect that the braiding of two 
objects depends only on the product of their charges. So braiding of two objects both of 
charge —1 (z's unaccompanied by w's) or two objects both of charge 1 (z's accompanied 
by w^s) gives the same result. Braiding of an object of charge 1 with an object of charge 
— 1 is less simple, since non-trivial gradient flows enter the picture, as we saw in section 
6.3.1. 

6.5 General Picture With Symmetry Breaking 

In this section, we develop a general picture with symmetry breaking. We take the sym- 
metry breaking parameter c to be real and positive; we consider any number of singular 
monopoles at positions Za, a = 1, . . . ,d, and any number of Bethe roots Wi, i = 1, . . . ,q. 
At first we will set all the charges ka to 1. 

We can produce a basis of integration cycles by hand. We assume that the imaginary 
parts of the Za are distinct, and ordered so that lioa{za — Za') > if a' > a. We define the 
rays Ca which start at Za and are parallel to the positive real axis. 

Our integration cycles will be products Ca^aj - a, := Ca^ x x • • • x Ca, for every subset 
of q distinct z's. This gives ('^) integration cycles, which in the region of large c can be 
interpreted as thimbles. Indeed the Bethe equations 

V^^ = c + V^— (6.18) 

^ Wi- Za ^.Wi- Wj 

have approximate solutions for large c with each wi equal approximately to + 1/ c, with 
Oj 7^ Qj for i / j. The corresponding thimbles are precisely the Ca^a2...aq- Summing over 
all g, we get 2'^ critical points or conformal blocks in all, as expected. 



-68- 



Figure 16. An integration cycle for one singular monopole of non-minimal charge. 



z < 
1 





Figure 17. An integration cycle for several singular monopoles of non-minimal charge. To Za we 
attach Qa of the w's, for suitable qa- In the picture, the qa are 1,3, and 2. 

The braid group representation is rather simple in this basis. We let the Za move 
around in the complex plane. The Stokes walls are approximately at the locus where the 
real parts of two Za coincide. Morse flows will occur if and only if a Za unaccompanied 
by a Wi passes to the right of a Za' accompanied hy a Wi. The resulting behavior involves 
only the two z's that are crossing and is precisely what we analyzed in detail in section 
6.3.1. The braiding matrix when Za crosses Za' acts non-trivially only on Ca and Ca' (and 
any cycle Caia2...aq that contains one or both of these), and takes the same form as (6.13). 
As in fig. 14, the braiding matrix for this process can be interpreted as resulting from a 
pair of Morse theory flows. It is tedious, but elementary, to check that the braid group 
relations are satisfied, as they should be. 

6.5.1 Degenerate Fields Of Any Charge 

Now, still with complex symmetry breaking, we will relax the constraint that all ka equal 1. 
First, we can consider a single z of charge k, with q Bethe roots. The Bethe equations only 
have solutions if k > q. (This is clear from eqn. (3.37) for the opers.) Actually, according 
to the theory of Bethe equations, the solution is unique, for given q, and corresponds to a 
unique thimble or integration cycle. 

Indeed, it is easy to describe this unique possible integration cycle C^'^^: one integrates 
each Wi from z to +oo, along distinct, non-intersecting paths, as in fig. 16. 
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We can then immediately describe a basis of integration cycles in the most general 
case, with any number of Za of charge ka, and the number of Wi being q < ka- We get 
a unique cycle for each way to decompose q = J2a la '^ith qa < k^. for each a, we integrate 
qa of the Wi from Za to +00, along distinct, non-intersecting paths, as in fig. 17. This gives 
an integration cycle Yla ^a^°^ ■ 

At strong symmetry breaking (or equivalently, if the Za are well separated along the 
imaginary axis), these integration cycles correspond to the thimbles associated to the 
unique solutions of the Bethe equations with qa of the Wi near Za- There is no obstruction, 
in principle, to derive the braid group representation in this basis. 

For a homological approach - essentially corresponding to the free field realization - 
to the construction of braid group representations for any ka, though without symmetry 
breaking, see [48]. 

6.6 Turning OfT Symmetry Breaking 

Since symmetry breaking is so useful in simplifying our analysis, the question arises of 
verifying that symmetry breaking does not affect the values of knot or link invariants. 

The description of integration cycles in the previous section is applicable as long as the 
symmetry breaking parameter c is non-zero. But for c — )■ 0, we lose some integration cycles, 
as the Wi cannot go to infinity any more. In terms of critical points of W, the behavior for 
small c is easy to describe. For small c, for each solution of the Bethe equations, the Bethe 
roots split naturally into two subsets. Some number go of Bethe roots, which we call wf'\ 
i = 1, . . . ,(7o, remain of order 1 in the limit c — )• 0, while the remaining q'oo Bethe roots, 
which we call Wj^\ j = ^, ■ ■ ■ ,Qoo, are of order c~^. 

In this situation, the Bethe equations for the wf^'' are well approximated by the Bethe 
equations for that number w^s in the absence of symmetry breaking. On the other hand, 
the Bethe equations for the Wj°°^ are well approximated by the Bethe equations for that 
number of w^s, in the presence of a single z of charge kes = ka — 2qQ. Eqn. (3.22) 
ensures that kes is non-negative, so in fact 



As we have discussed in section 6.5.1, the Bethe equations for the w^ have a single 
solution if k^s > ^oo; otherwise, they have no solutions. Summing over all decompositions 
Q = Qo + Qoo and all solutions of the Bethe equations for the wf'\ and finally over all 
possible values of q, one gets the expected number ria=i(^a + ^) °f solutions of the Bethe 
equations. 

The counting can be carried out as follows, in terms of Wilson operators of the dual 
SU (2) gauge theory. The singular monopoles of charge ka correspond to representations Ra 
of SU (2) of spin ka/2 and dimension ka + 1. The tensor product R = <i$aRa has dimension 
na(^a + This representation can be decomposed as a direct sum of SU{2) modules 
of spin kcs/2 (where feeff is bounded by < k^s < J2a^o., just as in (6.19)). The states 
that transform with this spin are as numerous as the solutions of the Bethe equations with 

ka — 2qo = kcs/2 (and all possible values of q). 




(6.19) 



a 
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Now consider braiding of the Za- For sufficiently small c, when one crosses a Stokes wall, 
there are Morse flows in which Qoo becomes smaller, but no such flows in which goo becomes 
larger. The reason for this is that for c — t- 0, the values of the Morse function h = Re W at 
a critical point are greater the greater is goo- (This is because some contributions to h are 
of order ln(l/|c|) for c — )• 0. With the help of (6.19), one can show that the coefficient of 
ln(l/|c|) is an increasing function of goo-) As a result, decomposing the space V. of physical 
states according to the value of goo, the monodromy matrix is block triangular: 



B 



* * * *\ 

* * * 
* * * 

yo */ 



(6.20) 



The diagonal blocks (which are of rank 2, 2, and 1 in the example given) are the monodromy 
representations that one would have in the absence of symmetry breaking for given go = 
q - Qoo- 

To the extent that one can compute knot or link invariants by taking traces of braid 
group representations, the off-diagonal blocks in (6.20) are not important as they do not 
contribute to traces. Actually, to compute the Jones polynomial and related invariants of 
knots and links, one needs in addition to the braid group representations an additional 
"fusion" operation in which a pair of Za of the same charge is created or annihilated. 
The additional information that we need to ensure that knot invariants are unaffected by 
symmetry breaking and do not change upon setting c = is that fusion never involves 
creating or annihilating any w's at inflnity. This is natural because of the local nature of 
the fusion operation. 

6.6.1 A Clarification 

A careful reader might notice a small sleight of hand in this derivation. The inequality 
^eflF ^ was deduced from (3.22), but the original derivation of (3.22) was based on picking 
Q to have a smaller degree than P. The Bethe equations at c = certainly have solutions 
in which this is not the case. Why are we entitled to restrict to this case? 
Suppose that at c = 0, we find a pair Pq, Qq obeying the oper condition 

pf._^Q = K{z). (6.21) 

az az 

Now suppose that we turn on very weak symmetry breaking. The oper equation becomes 

^dQ _ dP^ _ ^ ^ 
az az 

We hope that as c is turned on, there is a pair (P(z; c), Q(z; c)) obeying (6.22) and such 
that Q has an expansion 

Q{z;c) =Qo + cQi + c^Q2 + .... (6.23) 
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This will ensure that the Q{z;c) has roots that approach the roots of Qo as c — )■ 0, plus 
possible additional roots that go to infinity for c — t- 0. One might expect that P would 
have an expansion of the same form, but this is not the case. The degrees p and q of 
polynomials P, Q obeying (6.21) satisfy p + q = k + 1, but as soon as c 7^ 0, the relation 
becomes p + q = k. So the degree of P must drop as soon as c 7^ 0. The way that this 
happens is that the expansion for P is actually 



where we write Pq for the coefficient of c , as this polynomial does not coincide with Pq. 
Plugging (6.23) and (6.24) in (6.22), we learn from the term of order in the equation 
that P-i is a multiple of Qq, and this multiple must be nonzero or else the term of order 
c" in the equation would force p + q > k + 1. So p > qo, and this together with p + q = k 
and q > qo implies that qo < k/2, as desired. 

The moral of the story is that solutions of the Bethe equations for c = with q > k/2 
do exist, but they are unstable to symmetry breaking. Various forms of this statement are 
known in the literature on integrable systems. 

6.7 Three-Dimensional Interpretation 

To apply our results to knots and not just to braids, it will help to understand the three- 
dimensional interpretation of what we have computed so far. We consider knots in M^, so 
the four-manifold on which we are trying to count solutions of eqns. (1.1) is M4 = x R_|_. 
We describe with Euclidean coordinates . The adiabatic evolution considered 

so far has been in the direction, while we have combined the other coordinates to a 
complex variable z = + ix^. As usual, we take the gauge group to be G = 5(9(3). 

We will focus on the case of strong symmetry breaking. The symmetry breaking 
involves the choice of an expectation value (j) = diag(a, —a), where a is a vector in R'^. As 
long as the complex symmetry breaking is nonzero, this vector does not point in the x^ 
direction, that is, the direction that we chose for the adiabatic evolution. Topologically, if 
the directions are not the same, we may as well think of them as orthogonal: we consider 
adiabatic evolution in the x^ direction, and symmetry breaking with c = —2a/C, pointing 
in the positive direction. This will be strong complex symmetry breaking with real, 
positive c, in the terminology that we have used so far. 

We will concentrate on the case that the strands have minimum magnetic charge only, 
and thus are dual to the two-dimensional representation of SU{2). At a generic time, each 
strand has two possible states: it is or it is not accompanied by a Bethe root Wi. In the low 
energy effective abelian gauge theory, the strand has magnetic charge 1 if accompanied by 
a Bethe root, and otherwise —1. 

The magnetic charge of a given strand changes when the Bethe root accompanying 
that strand moves to another strand. In the context of adiabatic evolution, this results 
from a Morse theory flow in which a Bethe root moves from one strand to another. This 
happens at a value of x^ at which one crosses a Stokes wall. The lesson of section 6.3.1 is 
that (in the limit of strong symmetry breaking) one crosses a Stokes wall at a time (that is 
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Figure 18. A three-dimensional picture of the process that leads to a Morse theory flow. Just two 
strands are pictured here. The direction is plotted vertically and the direction runs into the 
paper. We are looking at the picture from along the negative x'^ axis. The coordinates x'^ and x^ 
combine to a complex variable z ~ x'^ + ix'^ . For a particular choice of the direction of complex 
symmetry breaking, a non-trivial Morse theory flow can occur only at values of x^ at which the 
two strands have the same value of x^ and thus project to the same point in the x^ — x'^ plane. 
In the language of knot theory, we make a two-dimensional picture by projecting a knot or link to 
the x^ — x^ plane. In this projection, there are crossing points, and these are the points at which a 
non-trivial Morse flow may occcur. 

a value of x^) at which two strands have the same value of = Imz. So the crossing of a 
Stokes wall occurs when two strands have common values of and x^, and thus differ only 
in x^. Differently put, this happens when the two strands are separated in the direction of 
symmetry breaking. 

A three-dimensional picture clarifies things (fig. 18). When two strands align along the 
direction, a Morse flow can occur. The Morse flow occurs on a time scale fast compared 
to the adiabatic evolution, so in the adiabatic picture it is essentially instantaneous. The 
flow involves a Bethe root moving towards the positive direction. The flow can only 
occur between strands of opposite charge, and will allow positive charge to move towards 
positive x^ only. 

In general, given any knot, we can usefully project it to the x^ — x'^ plane, and look at it 
from the negative x^ direction. We suppose that the embedding of the knot in is generic 
enough so that its tangent vector always has a non-zero projection to the x^ — x^ plane, 
and moreover so that the projection of the knot to the plane has only simple crossings; 
finally we will assume that the function x^ has only simple maxima and minima along the 
knot. A simple example of a knot projection is given in fig. 19. In such a knot projection, 
the low energy abelian description is valid away from crossings, so away from crossings and 
local maxima and minima, which we discuss in section 6.7.1, each strand can be labeled 
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Figure 19. A simple example of a knot projection, with only simple crossings and simple maxima 
and minima of the height function. This figure also illustrates the fact that the projection to a plane 
of an oriented knot allows one to define an integer invariant p = (l/27r) f dsdO/ds that equals the 
total change in moving around the knot of the angle 9 defined by the tangent vector to the knot. 
For the example shown, p — 2. p is the only invariant of a knot projection that can be written as a 
local integral along the knot. It depends on the choice of projection and is not an invariant of the 
knot per se. 




+ - - + + - --f 



Figure 20. The vertex model assigns the indicated factors to every crossing of two strands. The + 
and — signs labeling the strands express upward flow of magnetic charge +1 or —1; if one turns the 
picture upside down, the weights remain unchanged, provided one exchanges all + and — labels. 
Charge can be exchanged between strands, but only when a positive charge flows in from the bottom 
above a negative charge. If one reflects the picture from left to right, while also replacing q with 
, the weights remain invariant; this reflects the behavior of Chern-Simons theory under reversal 
of orientation. 

by its magnetic charge 1 or —1. The charges are unchanged at crossings, unless a strand 
of positive charge passes over a strand of negative charge, in which case a charge exchange 
process (corresponding in Morse theory to a non-trivial gradient flow) is possible. 

Our previous calculations assign a weight to each possible crossing. These weights 
are just g^^/^ if the charges are unchanged at the crossing. When the charge jumps, the 
weight is q^^^'^ — q^'^^'^. The weights are summarized in fig. 20. We have arrived at a 
known vertex model representation of the braid group representations associated to the 
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Jones polynomials. See for example the i2-matrix on page 125 of [49] (where A is our q^^^) 
or see [50], especially pp. 1777-8, or fig. 10 of [51]. 

6.7.1 Creation And Annihilation Of Strands 

In order to reproduce the knot invariants, we need to understand the loci where the adia- 
batic approximation is invalid, because two strands are created or annihilated. Although 
the adiabatic approximation is invalid near such points, the low energy abelian description 
remains valid. As an immediate consequence, conservation of charge in the abelian theory 
makes it clear that only pairs of strands with opposite charge can be created or annihilated. 

The map from line operators in a microscopic theory to line operators in an effective 
low energy description is akin to the corresponding map for local operators, but it has 
a little twist: the coefficients are not c-numbers, but rather quantum mechanical vector 
spaces that have to be transported along the line. In the present case, the vector spaces 
are one-dimensional (in the U{1) theory, an 't Hooft operator has no structure except its 
charge), but we can still get an overall factor from parallel transport. This factor has to 
be written locally along the loop, and must also be consistent with topological invariance. 
For a knot without any additional structure there is no topological invariant that can be 
written as a local integral along the knot, but once one is given a projection of the knot 
to a plane - in our case the — plane - there is precisely one such invariant, the 
total winding number of the tangent vector to the knot (fig. 19). This can be written as 
{1/2tt) <f dsdO/ds, where the knot is parametrized by a variable s, and 9 is the angular 
direction of the tangent vector to the knot in the x^ — x^ plane. To define the sign of this 
invariant, one needs an orientation of the knot, which in our case comes from the direction 
of flow of magnetic charge. When the microscopic SU{2) theory is approximated at long 
distances as an effective U{1) theory, the effective action for an 't Hooft operator may 
acquire a term —ir]<fdsd9/ds, with a universal coefficient ij. If the tangent direction to a 
knot changes by an angle A9, this will contribute a factor 

exp{ir]A9) (6.25) 

In creation or annihilation of a pair of strands, the change in the tangent angle is = vr 
or — vr, depending on whether the positive charge bends to the left or to the right. This 
effect will associate a universal factor exp(ibz7rr/) to each creation or annihilation event, 
depending on the direction of flow of charge. 

Up to sign, there is a unique choice of r] that leads to a knot invariant, namely 
exp(iTn]) = =Fi(7^/^. The choice of sign does not matter, since every knot has an even 
number of creation and annihilation events; we will take exp(z7r?7) = —iq^^^. The weights 
for creation and annihilation events with this value of r] are shown in fig. 21. The value of 
rj could possibly be computed directly by studying the four-dimensional gauge theory BPS 
equations near the abelian limit. The factor of q^^^ should express the difference between 
instanton number computed microscopically in the 5*0(3) theory and instanton number 
computed in the low energy U{1) theory. The factor of =F^ should involve a comparison 
between fermion determinants for S0{3) and for U{1). Globally, the contribution of a 
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Figure 21. The weights of the vertex model for creation or annihilation of a pair of strands. As 
in fig. 20, the weights arc invariant under rotating the picture upside down if one exchanges -f and 
— labels, or under a reflection from left to right if one replaces q by q^^ ■ 



given classical solution to the Jones polynomial is proportional to the sign of the fermion 
determinant, a subtle invariant that may receive contributions from charged modes in the 
microscopic 5'0(3) theory; to write this sign as a product of local factors, one apparently 
must use factors of i^^, with overall signs that depend on how one trivializes the determi- 
nant line bundle. 

The value of r/ can actually be deduced by combining the information which is available 
in the abelian description and information available from the conformal block description. 
When a pair of strands is created or annihilated, we expect them to be fused to the identity. 
This means that the two nearby strands, located at say zi and Z2, are accompanied by a 
Bethe root w, and that the integration cycle for this Bethe root is the thimble C12 connecting 
zi and Z2- On the other hand, in section 6.3.1, we also defined integration cycles Ci and 
C2 with the property that the strands at zi and at Z2 have definite magnetic charges. 
(For example, in Ci, the Bethe root accompanies zi, so the charges of the two strands are 
respectively 1 and —1. We order zi and Z2 in order of decreasing = Imz.) The relation 
among these cycles turned out to be 



C 



12 



- q^l^C2. 



(6.26) 



The ratio of the ampitude to create a pair of charges (1,-1) to the amplitude to create 
a pair of charges (—1,1) is the ratio of the coefficients on the right hand side of (6.26), 
or — On the other hand, in the abelian description, this ratio is exp(27ri7y). So 
exp(i7rr7) = =Fiq~^/^, as shown in fig. 21. 

We have arrived to what is essentially a standard vertex model algorithm for calculating 
the Jones polynomial: given a knot or a link, pick a projection to the x^ — x^ plane such that 
there are only simple crossings and the function x^ only has simple maxima and minima. 
Divide the link into segments separated by the maxima, minima and crossings. Label 
the segments by it and sum over all labelings, weighting each labeling with the product 
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Figure 22. An unknot projected to the plane in tlie most obvious way. Tiiere are no crossings, 
but there is a creation event and an annihilation event. In the vertex model, the invariant for 
the unknot is computed by summing over all ways to label the two sides of the knot (that is, the 
segments between crossing, creation, and annihilation events) by charges + or — . Each labeling 
is weighted by the product of the appropriate local factors. In the present example, only the two 
choices shown make nonzero contributions, leading at once to the result —ci^^^ — q~^l'^ . 

of the local weights at crossings, maxima, and minima. Some simple examples are given 
momentarily. 

6.7.2 Some Examples And Some Topological Details 

For the simplest example of the use of the vertex model, we compute the expectation value 
of an unknot, projected to the plane in an obvious way (fig. 22). Summing over the two 
possible labelings of the diagram, we arrive at the result —q^/"^ — q~^/'^ . 

For a slightly less trivial example, we consider the two projections of a single strand to 
the — plane shown in fig. 23. In either (a) or (b), the knot projections shown on the 
left or right can be deformed into one another, so one might expect them to be equivalent. 
But in the present context, this is actually not the case. The twist of (a) relative to (b) 
introduces a factor of —q^^^, which can be evaluated by making use of the weights of the 
vertex model. It is instructive to actually do this; the same factor —q^^^ arises whether the 
magnetic charge of the strand is +1 or —1, but the details of the computation are quite 
different in the two cases. 

Since the factor —q^^^ does not depend on the magnetic charge carried by a given 
strand, it is universal: adding a twist of the type shown in the figure to any strand in an 
arbitrary knot or link multiplies the associated invariant by —q^^^. A similar twist of the 
opposite handedness multiplies the invariant by —q~^^^, for similar reasons. This factor 
means that the invariant associated to a knot (or link) by the quantum field theory depends 
on a choice of "framing." A framed knot is a knot that is slightly thickened into a ribbon. 
One keeps track of how the ribbon is twisted and (in the present context) adding a twist 
multiplies the invariant by a factor of — g'^/^. A knot that is presented with a projection to 
a plane comes with a natural framing, given by a slight thickening in the vertical direction, 
normal to the plane. A little thought (or experimentation with a strip of paper) shows that 
although the pictures on the left and right of fig. 23(a) or (b) are topologically equivalent 
if one ignores the framings, they do differ by one unit of framing. 

In the context of three-dimensional knot invariants that are associated to two-dimensional 
conformal field theory, conformal primary fields in two dimensions are associated to line 
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Figure 23. Use of the vertex model to compare two different projections of a single strand to 
the — plane - with a twist (left) or no twist (right). For either sign of the charge carried by 
the strand, the twist introduces a factor of —q^^'^. The computation is quite different for the two 
possible values of the charge carried by the strand; for charge —1, as shown in (a), the vertex model 
sum has only one nonzero contribution, corresponding to the indicated labeling, but for charge +1, 
there are two possible contributions, shown in (b); they add to the same result. In (b), the second 
contribution involves a charge exchange process in which upper and lower strands exchange charge 
where they cross. A similar twist with undercrossing instead of ovcrcrossing (or with the whole 
picture replaced by a mirror image) leads instead to a factor of —q^^^^, as the reader can verify. 

operators in three dimensions. If a conformal primary has dimension h, then in a unit 
change in framing, the corresponding line operator is multiplied by exp(27ri/i). The fac- 
tor — g^/^ is indeed exp(27ri^), where the conformal dimension of the degenerate Virasoro 
primary field V_f^/2b is 

, k k(k + 2) , , 

V(fc) = ---^^, (6.27) 

and in addition q = exp(— 27rf/6^), and the vertex weights that we have described are for 
the minimum charge case k = 1. 

At this stage, an important detail arises. The vertex weights that we have described 
are appropriate for a certain natural normalization of the Jones polynomial, which has 
been used in the literature, for instance in [49]. However, a shghtly different normalization 
arises in SU (2) Chern-Simons theory. In Chern-Simons theory, the expectation value of 
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an unknot labeled by the two-dimensional representation of SU{2) is g^/^ + which 
differs in sign from what we deduced in fig. 22 using the vertex weights. (The sign is easily 
checked in Chern-Simons theory. Since q = exp(27ri/(fc^ + 2)) in SU{2) Chern-Simons 
theory, where fc^ is the level, the weak coupling limit fc^ — )■ oo corresponds to g = 1; 
for q = 1, the expectation value of a Wilson loop in any representation is simply the 
dimension of the representation, or +2 for the two-dimensional representation.) Similarly, 
the dimension of a primary field related to a representation of SU{2) of spin j = k/2 is 

so that the factor acquired in a unit change of framing is exp(27ri/ics) = Q^^^, without 
the minus sign found in fig. 23. Clearly the discrepancy in sign reflects the fact that 
hv{k) — hcs{k) = —k/2, independent of b and k'^ . 

The comparison with Chern-Simons theory is not necessarily a problem for the present 
paper, in which we have simply started with the four-dimensional gauge theory equations 
(1.1). However, one would like to know the best interpretation of these signs in the context 
of the duality presented in [14] between Wilson operators of Chern-Simons theory and 
singular monopoles at the boundary. We believe that the interpretation may be that the 
dual of a Wilson operator of spin fc/2 in Chern-Simons theory is actually a boundary 't 
Hooft operator that carries angular momentum k/2 and is fermionic when k is odd. (For 
't Hooft operators defined on the boundary of a four-manifold, the relevant rotation group 
is 50(2) or rather its double cover Spin(2); this group is abelian and has one-dimensional 
representations, labeled by the angular momentum k/2.) The fermi statistics for odd k 
would give a minus sign for every crossing (relative to what is presented in fig. 20) and 
a minus sign for every closed loop; including these signs brings the results obtained by 
braiding of Virasoro degenerate fields in agreement with the results obtained by braiding 
in Chern-Simons theory. 

Finally, it is instructive to compare the computation in fig. 23 to an equivalent com- 
putation if the gauge group were simply G = U{1) instead of 50(3). There would be two 
differences. First, the factors in fig. 21 associated to creation and annihilation of a pair 
of strands would simply be 1. (Those factors come entirely from integrating out massive 
degrees of freedom of the 50(3) theory, in reducing to an effective abelian description at 
low energies.) Second, charge exchange processes are absent for G = U{1) (as there are no 
smooth monopoles), so the second contribution in fig. 23(b) would be absent. A look back 
to fig. 20 shows that for G = U{1), evaluation of either fig. 23(a) or (b) gives a simple 
factor of q^^^, instead of —q^^^. Two comments are in order: 

• The minus sign of fig. 23 is absent for G = U{1) (and similarly the minus sign in fig. 
22 is absent). 

• In Chern-Simons theory, the power of q is the quadratic Casimir invariant of the 
relevant representation of G^. The quadratic Casimir of a representation of highest 
weight j is f for O^ = U{1) and j{j + 1) for O^ = SU{2). For j = 1/2, this gives 

Qp g3/4 ^j^g abelian and nonabelian cases, respectively. 
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Figure 24. To make the charge exchange process more visible, we have exchanged the coordinate 
axes relative to fig. 18. The axis now runs vertically while runs horizontally. Charge exchange 
occurs at values of x^ at which two strands differ only in the value of x^, so with the coordinate 
axes aligned as in this picture, the charge exchange invoves a flow of charge in the vertical direction, 
represented by the dotted line. We take this to represent the propagation of a BPS soliton. The 
soliton propagates along the axis of symmetry breaking, so it is described by a solution with real 
symmetry breaking only - in fact, by the bare Miura oper of equtation (3.31) with no 't Hooft 
operator and a single Bethe root. 



6.7.3 Gradient Flow And Strings 

We can give a more concrete physical interpretation to the gradient flows which occur when 
strands of appropriate charge cross. Propagation of magnetic charge from one strand to 
another can be described by motion of a magnetic monopole between the two strands. If 
we rotate the coordinates and think of the direction as "time," we should be able to see 
the relevant object as a time-independent solution in the presence of only real symmetry 
breaking. 

In fact, we have already described precisely the necessary solution: it is associated to 
the bare Miura oper with a single Bethe root and no singular monopole that was described 
in eqn. (3.31). So we visualize the charge exchange process as propagation of an object in 
the direction, as in fig. 24. This picture is oversimplified, as it ignores the existence of 
a fourth dimension, normal to the boundary that contains the knots. An alternative 
picture showing the role of the y direction is given in fig. 25. When propagating in the 
x?' direction, the soliton settles at a value of y that is given by the solution for the bare 
Miura oper. (How to compute this value is explained most fully in section 7.) Of course, 
this description is only good if the soliton propagates far enough in that it has "time" 
to reach the equilibrium value of y. So it is only good if the strands are sufficiently far 
separated, or the symmetry breaking is strong enough. 

We have gained an intuitive picture of the solutions of the four-dimensional BPS 
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Figure 25. An alternative view of the charge exchange process that was pictured in fig. 24, to 
show the role of the y direction (here depicted as the direction normal to the plane that contains 
the two strands). A soliton propagating between two boundary points A and B that are separated 
by a long distance in the direction will bend away from the boundary, to reach a value of y 
corresponding to the solution for the bare Miura oper. 

equations at strong symmetry breaking: they describe smooth monopole configurations 
stretched between the singular monopole strands. The contribution of each configuration 
to the knot invariant is then the product of two types of factor. One type arises from 
the map from the microscopic nonabelian theory to the abelian theory, while the other is 
computed in the abelian theory. Factors of the first type appear where pairs of strands are 
created or annihilated and where a smooth monopole is emitted or absorbed by a boundary 
singular monopole. In the abelian description, the smooth monopoles also look like singular 
monopoles, of charge 2, but not attached to the boundary. The whole configuration looks 
like a web of monopole strands, to which the abelian theory associates an overall power of 

6.7.4 Analog in the Dual Chern-Simons Theory 

It is entertaining to carry this picture all the way back to three-dimensional Chern-Simons 
theory, before all the dualities which brought us to the four-dimensional description studied 
in the present paper. A key step in relating the two pictures is 5-duality, which luckily is 
very transparent at strong symmetry breaking, as it reduces to electric-magnetic duality 
in the abelian gauge theory. We get immediately a sum over configurations of massive 
VF-bosons stretched between Wilson lines, and only carrying electric charge towards the 
positive direction. 



This suggests to look for a gauge condition in nonabelian Chern-Simons theory which 
would have this effect. It is easy to indentify it: one can pick a "partial" axial gauge fixing 
which reduces the A2 component of the gauge field to the Cartan subalgebra. We write 
the nonabelian gauge field as A = Bi^, + W+t+ + W-i-, and impose the = gauge 
condition for W±; it will not be necessary here to make a gauge choice for the diagonal 
gauge field B. The action for charged M^-bosons then reduces to 

^ / d^x W+D2W., D2 = ^ + [B2, ■ ]. (6.29) 



The equation for the propagator is 
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D2Gix,y) = —6Hx-y) (6.30) 



and has a solution 



2w 



exp n ^ B2j 6{x^ - y^)5{x^ - y^)e{x^ - y^), (6.31) 



which only describes propagation of charge towards the positive direction. It is pretty 
clear that such a partial gauge fixing will lead to a version of the vertex model, though 
it may be tricky to compute the precise vertex weights. The necessary computations are 
likely to be somewhat similar to those involved in studying Chern-Simons theory in an 
ordinary axial gauge - for example, see [52, 53] - or in a certain almost axial gauge [54]. 
Somewhat analogous is the use of a complex version of axial gauge to derive the Knizhnik- 
Zamolodchikov equations [55]. 

6.8 Up to Six Dimensions 

With an eye to future categorification, it is useful to lift this abelian picture all the way to 
six dimensions. After all, in the presence of strong symmetry breaking, the six-dimensional 
(0, 2) theory is not that mysterious. It is a theory of self-dual two-forms coupled to heavy 
dynamical BPS strings. The six-dimensional setup which leads to knot invariants involves 
the six-dimensional theory on the product M x M3 x D, where D is a copy of with 
the geometry of a semi-infinite cigar (see eqn. (7.28)). The knot itself is represented by a 
knotted two-dimensional defect placed at the tip of the cigar. The g- grading comes from 
the conserved angular momentum derived from the rotational symmetry of D. 

In the absence of symmetry breaking, the fact that the metric on D is cigar-like rather 
than being the Euclidean metric is important in order to get a well-defined space of states 
for the defect. A defect placed in flat six-dimensional space would be strongly coupled to 
the bulk SCFT. On the other hand, in the presence of symmetry breaking the IR physics 
is free, and we can hope to recover the space of states from bound states of the dynamical 
heavy strings and the defect. The cigar geometry would then not play a significant role, 
and D can be replaced by a flat M^. 

Now, we will specialize to M3 = M^, and for brevity we will take the knot to be time- 
independent. The BPS condition for a time-independent dynamical string is very simple: 
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it must be straight, and aligned with the symmetry breaking (which we will still take to 
be the direction). This is literally the six-dimensional lift of the condition satisfied by 
the smooth monopoles, and leads to the same vertex- model picture when projected on . 
The defect strands have two possible ground states in the abelian low energy description, 
of opposite two-form charge. A junction with a dynamical string allow the two-form charge 
to jump. 

The interesting question is to reproduce the weights of the vertex model from six- 
dimensional considerations. This includes both contributions from the abelian theory of 
self-dual two-forms sourced by the strings, and contributions from the worldvolume theory 
of the strings. We will not compute the former here, but we can give some insight on the 
latter. 

As summarized in fig. 20, the vertex model weight for a charge exchange process 
contributing to the Jones polynomial is 

^±1/4 _ ^T3/4^ (g 32) 

In the six-dimensional picture, the Jones polynomial is supposed to come from a sum 
over BPS states, weighted by 1)^, where P is the conserved charge that corresponds 
to rotation of D, while F, which one might loosely call fermion number, is a certain K- 
symmetry generator. We interpret the relative factor —q^^ between the two terms in (6.32) 
to mean that a BPS string connecting two strands in a knot has two physical states, differing 
by 1 in both angular momentum (to account for the factor of ^^^) and fermion number 
(to account for the minus sign). To understand this, we can focus on the two crossing 
strands, and the string stretched between them. Crucially, a single strand is half-BPS, 
two non-parallel strands are quarter-BPS, but the configuration of two strands exchanging 
a BPS string breaks symmetry further down to eighth-BPS. So the string breaks two 
supercharges, and hence it must have two ground states, exchanged by the action of the 
broken supercharges. This immediately gives the desired difference in quantum numbers. 

Hence the space of approximate ground states for the system can be described as 
follows. Let S be the set of configurations of the vertex model (labelings of strands by -|- 
or — , with smooth strings attached at crossings where a label jumps). For every smooth 
string, introduce a two-dimensional Hilbert space with quantum numbers derived from 
the last paragraph. To each s E 5, introduce a Hilbert space Us that is defined as the 
tensor product of the two-dimensional factors associated to the smooth strings. Then an 
approximation "Hq to the space of BPS states is "Hq = ®s^s^s- The grading of (by 
angular momentum and fermion number) is affected by the abelian field configurations 
sourced by the system of strings; for example, the self-dual abelian tensor fields can carry 
angular momentum. 

In order to compute Khovanov homology, one needs to evaluate the differential acting 
in this space of approximate ground states, by searching for instanton configurations which 
interpolate between different states in the past and future. The BPS condition for a time- 
dependent BPS string is still rather transparent. The worldsheet should be holomorphic in 
complex coordinates x'^-|-ix^ and -l-zx^. Notice that a time- independent string stretched 
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along the direction (and thus parametrized by and x^, with and x^ fixed) indeed 
has a holomorphic worldvolume. 

Thus we expect to be able to build the differential for Khovanov homology from the 
data of holomophic curves in R x M3 which end on the knot. This avenue seems promising 
for future development. 

7 An Effective Superpotential For Monopoles 
7.1 Overview Of Results 

Starting in section 2.4, we interpreted solutions of three-dimensional supersymmetric equa- 
tions in terms of configurations containing smooth BPS monopoles. However, the consider- 
ations were purely qualitative. In this section, we will make the reasoning quantitative. We 
will construct an effective superpotential for smooth BPS monopoles on x M_|. interacting 
with a Nahm pole and singular monopoles on the boundary. This effective superpotential 
will account in a direct way for all qualitative results from sections 2 and 3 about what 
solutions to our equations do or do not exist for various values of the parameters. Also, 
by integrating out some massive fields, we will be able to recover the Yang- Yang function 
(3.52) that has been one of our main tools. 

Let us first consider our underlying supersymmetic equations 

{F - (j) ^(|) + tAA(|))'' = 
{F - (j) A (p-t-^d A(l))~ = 

dA*</' = 0, (7.1) 

on M X (where the first factor is the "time" direction) and ask how the solutions of Bo- 
gomolny equations for smooth monopoles can be embedded as solutions of these equations. 
This is possible precisely if t = 1 (or —1), the only nonzero component of (j) is the time 
component, which here we will call 0j (rather than (pi^ as before), and we also set the time 
component of A to zero. Then the equations (7.1) reduce to the Bogomolny equations 

F = *dA0t. (7.2) 

For the Bogomolny equations to have smooth monopole solutions, the field 4>t must have 
an expectation value at infinity, which means that the real symmetry breaking parameter 
ai of section 2.4 must be nonzero. On the other hand, the complex symmetry breaking 
parameter a of section 2.5 must vanish (or (pt would not be the only nonzero component 
of </.). 

The basic solution of the Bogomolny equations on is the one-monopole solution 
for G = SO{3) or SU{2). It has has magnetic charge m = 2 (in units in which the basic 
singular monopole has charge m = 1). The moduli space of the one-monopole solution 
is P = X S^, where measures the center of mass position of the monopole on the 
spatial manifold M^, and is parametrized by a collective coordinate 1} for the U{1) 
gauge symmetry that is unbroken at infinity. Of course, V is a hyper-Kahler manifold; this 
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follows from the unbroken supersymmetry of the Bogomolny equations on W^. However, 
since we will be considering perturbations that break some of the supersymmetry, for our 
purposes it is more useful to merely look at P as a complex manifold in one of its complex 
structures. We decompose as x M (which we will eventually replace with x M_|_). 
The motion of the smooth monopole along M? is parametrized by a chiral superfield W. 
And the position y of the smooth monopole in the M direction combines with the collective 
coordinate i9 to a second chiral superfield 



Now we want to construct an effective superpotential W(W^, Y) that describes this 
situation. If a = C = 0, there is a smooth monopole solution for every value of W and 
Y, which means that every value of W and 1" is a critical point of W, so W must vanish 
(modulo an irrelevant constant). On the other hand, if either a or is nonzero, then 
there is no supersymmetric monopole solution, meaning that W has no critical point. But 
turning on a and (" preserves the symmetries of adding a constant to W or Y. So W must 
be invariant modulo an additive constant under constant shifts of W or Y; in other words, 
W must be a linear function of W and Y. In fact, the form of W is 



This follows from the following considerations. Invariance under rotations of M? implies 
that a and W can only appear as the product aW. On the other hand, although at ^ = 0, 
the superpotential W is given microscopically by the single- valued function (2.6), as soon 
as C becomes nonzero it is given by the Chern-Simons function (3.6), which is only single- 
valued mod 27rzZ. Since Y is similarly single-valued mod 27riZ (because of the angular 
nature of •&), the effective superpotential (7.4) has a multivaluedness that just matches 
that of the microscopic description. 

Now let us turn off a and ( but replace x M by x M_|., where M+ is the half-line 
y > and we assume the usual Nahm pole at ?/ = 0. We assume symmetry breaking for 
y — )• oo with (pt — )• diag(ai, — ai) while of course A and (at a = 0) 4> vanish for y — )• oo. 
In any solution, the fields approach these asymptotic values exponentially fast for y — oo. 
This reflects the fact that near y = oo, the gauge symmetry is reduced from SU (2) to U (1) 
by the expectation value of (pi, and all charged fields have masses proportional to |ai|. In 
particular, in the absence of singular or smooth monopoles, the solution for the full system 
(7.1) is given by a solution of Nahm's equations (the relevant solution is described in [16]) 
in which the fields approach their vacuum values exponentially fast for y — )■ oo. 

Next, still with a = C = 0, let us add a smooth monopole with positions W, Y. There 
is not an exact solution for the smooth monopole, because the Nahm pole forces charged 
components of (p to have nonzero values that "repel" the monopole to y = oo. However, 
these charged fields are exponentially small for large y, so a smooth monopole located at 
large y is exponentially close to being a solution. This being so, a smooth monopole that 
is located at large y must be governed by an effective superpotential. This superpotential 
must be a single-valued function of Y that vanishes exponentially for y — )• oo. These 



Y = aiy + I'd. 



(7.3) 



W = aW + (Y. 



(7.4) 
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conditions are satisfied by a linear combination of exponentials exp(— ny), n = 1,2,3 ... . 
However, it will soon become clear that the expected qualitative picture emerges if and 
only if at a = = 0, W is linear in exp{—Y): 



W = Aexp(-y) 



(7.5) 



for some constant A. 

In principle, it should be possible to compute this result by evaluating the microscopic 
superpotential (2.6) for an approximate solution consisting of a smooth monopole at large 
y in the presence of a Nahm pole. The exponentially small term should come from W 
boson exchange between the boundary and the monopole. Instead of attempting such a 
computation, we will take a shortcut in this paper: we will consider a representation of the 
smooth monopole and the Nahm pole by a configuration of branes, and in that context the 
exponential superpotential (7.5) will emerge from a simple brane instanton. 

Postponing that analysis to section 7.3, let us discuss the implications of (7.5). First 
we assume that a = ^ = 0, so that (7.5) is the full superpotential. We see at once that W 
has no critical point so there is no super symmetric solution in the presence of the smooth 
monopole. This is in full accord with the analysis in section 2.4. With real symmetry 
breaking only and ^ = 0, a supersymmetric solution is expected to be unique (this remains 
true even in the presence of singular monopoles at y = 0, which we have not yet included 
in W) and does not require smooth monopoles. Now let us see what happens when we 
turn on a and ^. We construct the full superpotential by simply adding the various terms 
that we have found so far: 



The justification for including a and C in this way is the same as before (A may now depend 
on (^). Let us look for critical points. We see at once that there is no critical point unless 
a = and C 7^ 0. If these conditions are satisfied, there is a one-parameter family of critical 
points, parametrized by W, with 



(The uniqueness of the solution for exp(— y) holds precisely because in (7.5) we took W 
to be linear in exp(— y).) We found a similar result in section 3.5.1 where we described in 
eqn. (3.31) a family of solutions that exist precisely if a = and C 7^ 0- The solution in 
question corresponded to a Miura oper with only real symmetry breaking, a single Bethe 
root at an arbitrary point i(j G M^, and no singular monopoles. We interpret the Bethe root 
w as the value of the superfield 11^ - in other words, the position of the smooth monopole in 
- while the position of the smooth monopole in the M4. direction is determined in (7.7). 
The reason that this description makes sense is that for small the smooth monopole is 
located at large y, where the effective superpotential W is valid. 

The next step is to include singular monopoles. As in section 2.2.1, in the presence of 
singular monopoles of charge ka located at positions Za, a = 1, . . . , d in the complex plane. 



W = aW + Cy + Aexp(-y). 



(7.6) 



exp(-y) = ^. 



(7.7) 
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it is convenient to introduce the polynomial K{z) = Y[t=ii'^ ~ ^a)^"" ■ In section 7.2, we will 
argue that the only effect of the singular monopoles is to multiply the exponential term by 
KiW), so that the superpotential becomes 

W = aW + CY + KK{W)e^v{-y)- (7.8) 

Now let us examine the implications of this formula. Suppose first that C = 0- Then the 
conditions for a critical point are 

K{W) = ^ = a + kK'{W)e^v{-y)- (7.9) 

The first condition says that W must equal one of the Za- For a = and all /cq = 1, the 
second condition then has no solutions. If /cq > 1 for some a, still with a = 0, the second 
condition is satisfied for arbitrary Y . All this is in keeping with what we found in sections 
2.3 and 2.4. Now suppose that a 7^ 0. li ka = 1, the second condition in (7.9) determines 
Y uniquely, and if a is small, then Y is large so that the analysis is valid. For /cq > 1, the 
second condition in (7.9) cannot be satisfied. All these statements match what was found 
in section 2.5 from a quite different point of view.^^ 

Now let us consider the case that C 7^ 0- The most illuminating way to proceed is to 
integrate out the massive field Y to generate an effective superpotential for W . For fixed 
W, the condition dW /dY = has the unique solution exp(— y) = C,/kK{W). Setting Y 
to this value and evaluating W, we find (modulo an irrelevant constant) 

W = aH^ + C log K{W) =aW + C^ka \og{W - Za). (7.10) 

a 

But this is the Yang- Yang function (3.52) for the case of a single Bethe root W = 
modulo terms that depend only on the Za and not on W] the present derivation is not 
sensitive to those terms. 

At this point, the reader hopefully would like to see a similar derivation leading to 
the general Yang- Yang function with any number of Bethe roots. For a general case 
with q smooth monopoles, we describe the positions of the i^^ smooth monopole by chiral 
superfields Wj, 1^, i = 1,. ..,(;. The definition of these fields is somewhat subtle and is 
discussed in section 7.2. The expectation values of the Wi will turn out to be the Bethe 
roots Wi of section 3.4. As in that discussion, it is convenient to introduce the polynomial 
Q{z) = Jli'=i(-^ ~ Wi). It turns out that the generalization of the superpotential (7.8) to 
an arbitrary number of smooth monopoles is 

w = «E^^ + cE^«+^E §St -m-y^). (7.11) 



In section 2.5, we found, for general ka and a 7^ 0, a solution with no smooth monopoles at W = Za 
and a solution with ka of them. (These correspond to the two ways of solving PQ = {z — Za)'''^ such that 
P and Q have no common zero at z — Za-) Since (7.8) is the superpotential for just one smooth monopole, 
it describes a solution with ka smooth monopoles only if fca = 1. The general analysis for arbitrary ka can 
be made and matched to section 2.5 using the superpotential (7.11) for an arbitrary number of smooth 
monopoles. 
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To recover the qualitative results of section 2, we first set C = 0, a 7^ 0. To find a 
critical point, the Wi must each equal zeroes Za oi K. Assuming that the charges ka are all 
1, no more than one of the Wi may equal the same Za (otherwise a zero of Q'{Wi) cancels 
a zero of K{Wi) and the condition dW/dYi = is not obeyed). Summing over all values 
of q, there are a total of 2"^ solutions - each Za may or may not be equal to one of the Wi. 
If we take a — )• 0, then all but one of these solutions (the one with no smooth monopoles 
at all) disappear, with Yi ~ log(l/a). 

For C / 0, just as in the derivation of (7.10), it is convenient to integrate out the 
massive fields Yi. Modulo terms that do not depend on the Wi, the superpotential that we 
arrive at is precisely the Yang- Yang function: 

W = Y,aWi+cY.^og{K{Wi)/Q'{Wi)) = Y,o-Wi+Cj2log{Wi-Za)-2CY,^og{W^-Wj). 

i i i i,a i<j 

(7.12) 

The attentive reader may notice one gap in what we have said. In the case = 0, 
we have not analyzed the problem with ka > 1 for some a. To do this, it is important 
to consider the case that Wi = Wj = Za for some but the coordinates that we have 
used to describe the monopole moduli space are actually not adequate when Wi = Wj. We 
explain a better description momentarily. 

7.2 Coordinates for Monopoles 

The moduli space of several smooth BPS monopoles on is the subject of a rich math- 
ematical theory [56]. From this theory we only need a small part: when the monopoles 
are widely separated in space (or equivalently when the symmetry breaking is strong), 
an effective abelian description of the moduli space is possible. In this description, the 
monopoles are regarded as "point" Dirac monopoles that interact with each other via the 
abelian gauge multiplet. Each monopole has a position in and an angular coordinate 
■d that is a collective coordinate for charge rotations. As in section 7.1, once we pick a 
particular complex structure on the moduli space, the position and angular coordinate of 
each monopole combine to a pair of chiral superfields W , y. However, there is a subtlety 
in the definition of 3^, which is the reason that we have changed our notation from section 
7.1. 

The angular part of the coordinate y parametrizes the freedom to do a U{\) gauge 
transformation on the smooth monopole solution before "gluing" it to the abelian solution. 
Hence is an holomorphic section of the U{1) gauge bundle over the VF-plane. This is 
the reason that in the presence of boundary singular monopoles, the exponential superpo- 
tential exp(— 3^) for a single smooth monopole needs a prefactor K{W): the superpotential 
should be a function, but the singular monopoles make into the section of a bundle 
®aO{za)~^'' , where the exponents are the charges of the singular monopoles in the abelian 
effective field theory. So we compensate for this by multiplying by K{W). A more intuitive 
explanation is that the superpotential encodes the interaction of the BPS monopole with 
the off-diagonal part of the complex Higgs field which has a zero of order ka at Za. So 
the superpotential acquires a factor {z — Za)'^°- . 
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On the other hand, in the abehan theory, the BPS monopole behaves hke a singular 
monopole. At any given value of y, one can restrict the U{1) gauge bundle of the low 
energy description to the W plane. As one increases y so that one passes the location of a 
BPS monopole, the U{1) gauge bundle on the W plane jumps. In other words, a second 
BPS monopole to the right of the first will feel the presence of a singular monopole of 
charge +2 at the location Wi of the first monopole. Hence we expect a superpotential 

K{W2) 
{W2 - Wif 

and similarly for several monopoles with increasing values of the real parts of 3^j: 



W = K{Wi) exp(-3^i) + ^ Z N2 ^M-y2) (7.13) 



This superpotential is equivalent to the relevant part of (7.11) as long as the Wi are 
distinct. The two are related by the change of variables 

This re-definition does not affect the part of the superpotential linear in a and Q, since 
^j3^j = Yli^i- However, it is the 3^i, not the Yi whose real parts are the actual positions 
of the monopoles in the y direction; moreover, the difference between the and the Yi is 
divergent when two or more Wi coincide. The superpotential expressed in terms of the J^j 
reproduces correctly the counting of solutions at ^ = for arbitrary values of the ka- For 
a = 0, if K{z) = 2;^° the superpotential is 

This function is extremized for arbitrary if the Wi are all zero and the number q of 
smooth monopoles is no greater than A;/2, because the prefactors have a zero of order at 
least 2 when the Wi are all zero. This reproduces what we found in section 2.3. If we do 
the same computation in terms of the 1^, we would seem to get solutions even when the 
number of monopoles at the origin is greater then ka/2, but the Yi are not good variables 
when the Wi coincide. We will leave the case a 7^ 0, C = to the reader. For C 7^ 0? the 
difference between the Yi and the is not important. 

It is interesting to match the coordinates in the low energy description to the exact 
nonabelian description of the monopole moduli space. The exact monopole moduli space is 
parametrized [57] by a scattering matrix for the operator Vy = Dy + i[(j)t, • ]. The scattering 
matrix takes the form 

where Q is a monic polynomial of order P and P are polynomials of degree up to (7 — 1, 
and i? is a polynomial of degree up to g' — 2. P and Q are necessarily coprime (this follows 
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from the condition QR — PP = 1), and both P and R are uniquely determined by Q and 
P. One can define coordinates Wj, Yi on the monopole moduh space by QiWi) = and 
PiyVi) = expYJ. However, this definition does not work well if the Wi are not distinct. 

To find a parametrization that works better as long as the symmetry breaking is strong, 
we can proceed as follows. For a single monopole, the scattering matrix takes the form. 



{z-Wi) e^i 
-e-^i 



(7.18) 



This expression can be matched naturally to the fact that at low energies, the smooth 
monopole can be approximated as a singular Dirac monopole. It tells us that there are two 
solutions il)± of the equation Vyif) = that behave as 

y « 

(z-t^i)e"iJ^/2 |lj jOj y»Q (7.19) 




^_ ~ e-^^y'"^ ^Jj y » 

~ (z - Vt^i)e-"i^/2 ^ g-yi+aij//2 y « (7.20) 

Here ("0-) is the unique solution which is small for y<<0(y>>0). A singular 
monopole solution in the abelian theory would have given the same exponential growth for 
y >> (y << 0), and the subexponential correction is due to the exponentially decreasing 
corrections for the smooth monopole solution. 

For a configuration of many well-separated smooth monopoles, the scattering matrix 
is a product 

5 = 5,5,_i---5i. (7.21) 

where Sa is the scattering matrix due to the a*'* monopole and the monopoles are taken to 
be ordered in the y direction. It is natural to parametrize the moduli space by 

This enables us to write P and Q in terms of the Wi and J^j, and finally, using PiWi) = 
exp(yj), to express the Yi in terms of the Wi and the 3^^. 

7.3 Realization Via M-Theory And Branes 

Here we will explain an M-theory approach to understanding the exponential superpoten- 
tial (7.5). We begin with some preliminaries. 
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7.3.1 A/-Theory Preliminaries 

The six-dimensional (0, 2) model of type Ai can be realized on a pair of parallel M5- 
branes. Thus, we begin on R^"*^ with coordinates x^, . . . , x^^, and we consider two M5-branes 
parametrized by x^, . . . ,x^ and located at = • • • = x^^ = 0. This system preserves 16 
global supersymmetries. Their generators can be understood as eleven-dimensional spinors 
£ that obey 

rori---r5e = r6r7---rioe = e. (7.23) 

Here the gamma matrices obey {r^,r^} = 2g^y. 

To simplify the picture, we introduce symmetry breaking, separating the two M5- 
branes in, say, the direction. So now we place one at j;^ = and the other at x^ = L, 
for some L. At low energies, this system is described by a pair of abelian tensor multiplets, 
coupled to BPS strings. The strings arise from M2-branes stretched between the two 
M5-branes. The string tension is T = Tm2L, where Tm2 is the M2-brane tension. 

We will consider a string whose world- volume is parametrized by x^ and x^, and that 
is located at specified values of x^, . . . The string is of course represented by an M2- 
brane that stretches from x^ = to x^ = L, though this direction will just factor out of 
the following analysis. 

The string described in the last paragraph preserves those supersymmetries whose 
generator obeys 

roTirge = e (7.24) 

as well as (7.23). Altogether, there are eight unbroken supersymmetries, corresponding 
to A/" = 4 supersymmetry in the two-dimensional sense. There is an 50(4) symmetry 
group rotating the x'^,...,x^ coordinates, and a second 50(4) symmetry, which we will 
call SO(4:)ji, that rotates x'^, . . . , x^°. 

However, we will soon modify the construction in a way that will break half of the 
supersymmetry and also reduce 50(4) x 5*0(4)^^ to a maximal torus. So it will help to focus 
on the relevant symmetries to begin with. We consider the M = 2 subalgebra consisting 
of supersymmetries that (in addition to the conditions already given) are invariant under 
a combined rotation of (say) the x^ — x^ plane together with an i?-symmetry rotation of 
the x^ — x^'^ plane: 

(r4r5 + r9rio)e = o. (7.25) 

Given (7.24), this is equivalent to 

(TsTs + FrFs) e = 0. (7.26) 

We note that the equations (7.25) and (7.26) are exchanged if we exchange W = x"^ + ix^ 
with Z = x'^ + ix^, and similarly exchange x'^,x^ with x^,x^'^. So in particular, our 
conditions on e are symmetrical between W and Z. 

The M = 2 supersymmetry algebra singled out by the above conditions has a U{1) x 
C/(l) group of i?-symmetries generated by J = r7r8 and J' = rgTiQ. We want to compare 
J and J' to i?-symmetry generators that we will call J+ and J_ that only act, respectively, 
on supersymmetry generators that have positive or negative two-dimensional chirality, in 
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other words, that obey xe = ±e, where x = ToFi is the two-dimensional chirahty. The 
conditions given above can be combined to give 

xe = JJ'e (7.27) 

and this implies that (with a suitable choice of sign for J+ and J_, and normalizing them 
so that they square to 1 on supersymmetry generators of the appropriate chirality) J and 
J' can be expressed as J = J+ + J_ , J' = J+ — J^. 

In particular, an exchange J ^ J' amounts to J± — )• ±J±, an operation known as the 
mirror symmetry automorphism of the Af = 2 algebra. The automorphism of the above- 
described N = 2 algebra that exchanges W and Z also exchanges J and J', so it is a mirror 
symmetry. Hence, if we view W as a, chiral superfield in the two-dimensional worldsheet 
theory of the string, we must view Z as a twisted chiral superfield, in the sense of [58]. 

7.3.2 Reduction To Gauge Theory 

So far we have half-BPS strings, but no gauge theory description of them. 

To get a gauge theory description, we compactify one direction, say the direction, 
on a circle of radius R. M-theory on a circle reduces at long distances to Type IIA 
superstring theory. The M5-branes become D4-branes and the theory on the D4-branes is 
at long distances a U{2) gauge theory, broken to U{1) x U{1) by the separation between 
the D4-branes. What is relevant to us is the SU (2) subgroup, broken at low energies to 
U (1). The string that was originally described via a stretched M2-brane is now represented 
by a D2-brane stretched between the two D4-branes. 

This situation has been studied in [59]. The D2-brane stretched between two D4- 
branes carries magnetic charge and corresponds to a smooth BPS monopole in the low 
energy SU{2) gauge theory. 

Before compactifying the direction, the low energy theory along the string was 
described by the chiral superfield W = x"^ + ix^ and the twisted chiral superfield Z = x"^ + 
ix^. After the compactification, we can replace Z by the single- valued field exjp{Z/R), which 
is still a twisted chiral superfield. However, for matching to the theory of BPS monopoles, 
another variable is more useful. The moduli of the BPS monopole corresponding to the 
D2-brane are the positions of the underlying M2-brane in x'^,x'^, and x^, and the dual of its 
position in x^. The duality in question is a T-duality in the two-dimensional effective field 
theory governing the string. This T-duality is a mirror symmetry in the two-dimensional 
sense. It replaces x^ /R with a new angular coordinate So while Z/R = x^/R + ix^ /R 
is a twisted chiral superfield, Y = x'^/R + i^D is an ordinary chiral superfield, just like W . 
Of course, the single- valued chiral superfield is not Y but e^. 

We conclude with two comments: 

• The fact that the angular coordinate "Q of the BPS monopole is T-dual to the angular 
position x^ is part of the relation between M5-branes on a circle and D4-brane gauge 
theory. The symmetry that rotates x^ becomes instanton number in the 4+1- 
dimensional gauge theory, while the symmetry that rotates •& is electric charge. 
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• The single- valued field 0, = exp{Z/R) can be understood as a map to C*; it can be 
neither zero nor infinity. In section 7.3.3, we modify the problem to make it possible 
to have = 0. 

7.3.3 Reducing On A Half Space 

To get a non-trivial superpotential, we will have to break some of the translation symmetries 
of the problem. In fact, we are interested in gauge theory on a half-space, so we want to 
restrict y = x'^ to be non-negative. 

The gauge theory problem studied in the present paper arises if and x^ parametrize 
not M X S""*^, as is the case in our presentation so far, but rather a copy of with a cigar-like 
metric 

ds2 = d?/2 + /(y)(dx5)^ (7.28) 

Here f{y) ~ y'^/R'^ for y — )• and f{y) — )■ 1 for y — )• oo. In fact, this was the starting point 
in the derivation in [14]. 

As a complex manifold, M? is the same as C, so we can now parametrize the x^ and 
x^ directions by a C-valued chiral superfield 0, which asymptotically at large y (but only 
there) can be written 

n = exp{Z/R), Z = x^ + ix^. (7.29) 

In contrast to the concluding remark of section 7.3.2, ^2 is now C-valued rather than C*- 
valued, and in particular there is no problem in having = 0. 

We will make use of this shortly in generating a superpotential. 

7.3.4 The Instanton 

We now want to describe an M2-brane instanton that will generate the superpotential that 
we are looking for. 

The instanton is supposed to correct the physics of a string that is parametrized by a 
worldsheet coordinate X = x^ + ix^. The string is located at definite values of W and 0, 
say W = Wo,n = Qq. 

We can understand qualitatively what sort of instanton can generate a superpotential. 
We consider a two-dimensional model with N = 2 supersymmetry whose chiral ring is 
generated at the classical level by the chiral superfield W and whose twisted chiral ring is 
generated by the twisted chiral superfield 0. The chiral ring is the ring of observables of a 
twisted -B-model, and the superpotential that we want to generate will give a deformation 
of this chiral ring. A hypothetical superpotential must be generated by configurations 
that preserve the i?-model supersymmetry. These are configuations in which the chiral 
superfield W is constant, while the twisted chiral superfield O is holomorphic as a function 
of the worldsheet coordinate. (It is a familiar fact that i?-model supersymmetry requires 
a chiral superfield to be constant. The fact that S-model supersymmetry allows a twisted 
chiral superfield to be holomorphic is mirror to the perhaps more familiar fact that A-model 
supersymmetry allows a chiral superfield to be holomorphic.) 
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The instanton that generates the superpotential is accordingly given hy W = Wq while 
O is a nontrivial but simple holomorphic function of X: 

^ = X-Xo, (7.30) 

(Note that this formula only makes sense because il. is allowed to vanish.) Here Xq is a 
constant, which we interpret as the instanton position; as always in instanton physics, to 
calculate physical amplitudes, it is necessary to integrate over the instanton moduli, which 
here mean Xq as well as some fermionic moduli associated to the supercharges under which 
the instanton solution is not invariant. 

The worldvolume of an M2-brane instanton is supposed to be a three-manifold. The 
three-manifold we want is just the product of the two-manifold S that was defined in (7.30) 
with the one-manifold < < L (all at x"^ = • • • = x^^ = 0). 

Since it is invariant under i?-model super symmetry, and has no moduli except what 
follows from translation invariance and supersymmetry (the parameter Qq corresponds 
roughly to a constant value that Q would have in the absence of the instanton), this sort of 
instanton will generate a superpotential. To understand just what superpotential will be 
generated, we use the asymptotic formula (7.29) and look at the disturbance in the string 
that is generated by the instanton, at great distances. For large values of y = x^, we can 
write 

exp{{y + ix^)/R) = no{X -Xq). (7.31) 

We see that as X circles once around Xq in the clockwise direction (at large values of 
\X — Xq\ so that the formula (7.31) is valid), x^ increases by 27rR. To produce this effect, 
the operator inserted at X = Xq must be a twist field. As the T-dual of x^ is the angular 
variable ??, a twist field is exp(— i??), and this must be the t?-dependence of a superpotential 
that captures the effects of the instanton. The holomorphic expression must therefore be 
W = Aexp(— y), where Y = y + I'd, for some constant A. This is precisely the result 
claimed in (7.5). 

For further confirmation, and also to check the sign in the exponent of W, let us 
consider the behavior of the field y at large distances, far from X = Xq. At long distances, 
the fiuctuations in y are described by a free-field path integral 

J Dy exp j dx° dx^ [Vyp j . (7.32) 

When the operator exp{—y/R) is inserted in such a path integral at a point X = Xq, the 
result is that at large distances, y/R grows as | log(X — Xq)\. But this is exactly what we 
see in (7.31). 

A final comment is that if the worldvolume dimension of the string were bigger than 
2, we would have considered the instanton as a fiuctuation around a vacuum defined by 
a limiting value of y (and all the other worldvolume fields) for X ^ oo, and we would 
have asked for the instanton to approach this limiting value at infinity. In two dimensions, 
because of the usual infrared divergences - which appear, for instance, in the logarithmic 
growth mentioned in the last paragraph ~ such a formulation is not valid. 
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8 Opers And Branes 



The purpose of the present section is to place some of the ingredients that have appeared 
in the present paper in a wider context. 

We continue to study = 4 super Yang-Mills theory, with a twist that preserves half 
the supersymmetry, on the four-manifold M4 = M x C x and with the usual Nahm 
pole boundary condition at the finite end of 1R+ . The novelty, compared to what has been 
said so far, is that we will view the problem from the point of view of compactification 
on C from four to two dimensions. In general [60, 61], assuming for simplicity that C has 
genus at least 2 (we relax this condition in section 8.4), compactification on C gives at low 
energies a two-dimensional sigma- model in which the target is Ain, the moduli space of 
solutions of Hitchin's equations [20]. The Nahm pole boundary condition must reduce at 
low energies to a brane in this sigma-model, and this brane must be half-BPS because the 
Nahm pole boundary condition is half-BPS in four dimensions. 

8.1 Back to t = 1 

We begin by analyzing the case t = 1. For simplicity, we take G to be SU{2) or SO{3). 
In section 2.2, we found that at t = 1, the Nahm pole boundary condition (in the absence 
of singular monopoles) describes a Higgs bundle {E, ip) ^ C that is endowed with a line 
sub-bundle L C E that is nowhere c/p-invariant. Viewing E as a rank 2 complex bundle of 
trivial determinant, the inclusion L C E is part of an exact sequence: 



Here we use the fact that, as E has trivial determinant, the quotient E/L must be isomor- 
phic to L~^. 

We view 99 as a holomorphic map E — >• E K, where K is the canonical bundle of 
C. We can restrict 93 to L, to get a holomorphic map L — )• E ® and then using the 
projection E — )• L~^, we get a map if : L ^ (X" K. The condition that L is nowhere 
(^-invariant means precisely that the map ip : L ^ ® K is everywhere nonzero. In 
other words, this map is an isomorphism. 

Tensoring with L, we learn that is isomorphic to K, so that L is a square root K^^"^ 
of K. If G = SU{2), a solution of the Nahm pole boundary condition involves a choice of 
^1/2, while if G = 50(3), since we really should be working with the adjoint bundle ad(ii^) 
rather than E, the choice of K^^^ does not matter. In what follows, we assume that either 
G = S0{3) or we have picked a particular square root of K. 

It is possible to make a non-trivial extension — t- K^^'^ — > — ?• K~^/'^ — t- 0, and we 
will exploit this fact in section 8.2. However, for Higgs bundles, we want £^ to be a direct 
sum since in the case of a non-trivial extension, the Higgs fields that we are 

about to write would not exist. If we write E in column form 







L-^ E ^ L 



1 



^ 0. 



(8.1) 




(8.2) 
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then up to an automorphism of E, a possible Higgs field tp takes the form 




(8.3) 



where g is a quadratic differential. To be more exact, we assume the upper right matrix 
element of (p to be nonzero as otherwise L would be (^-invariant (and the Higgs bundle 
{E,ip) would be unstable, as explained in [20]). Given this, by a bundle automorphism 
diag(A, A~^), we can take the upper right matrix element to be 1, and by a lower triangular 
bundle automorphism, we can make the diagonal matrix elements of ip vanish. Finally, for 
E as in (8.2), the lower left matrix element of is a quadratic differential (an element of 
H°{C,K'^)), which we call q. 

Let T C Mh be the submanifold parametrizing the Higgs bundles {E,ip) described 
in the last paragraph. At t = 1, the brane in Mh defined by the Nahm pole is supported 
on T. What sort of subvariety is T? As in [20], let / be the complex structure on A^/f in 
which it parametrizes Higgs bundles, J the complex structure in which Ain parametrizes 
fiat bundles with connection A = A + and K = IJ. The Hitchin fibration is the map 
from A4h the space of quadratic differentials that maps {E,ip) toTi(p^. This map IS 
holomorphic in complex structure I. For the Higgs field in (8.2), we have Tnp'^ = 2q, so 
there is a unique such ip for every desired value of Tr 99^ . Accordingly, T is a holomorphic 
section of the Hitchin fibration; in fact it is the holomorphic section constructed in [20]. 
Actually, T is complex Lagrangian from the point of view of complex structure /. That 
assertion means that the complex symplectic form 



vanishes when restricted to T- This is the case since, as the holomorphic type of E is fixed 
for all Higgs bundles that represent points in T, 6Az is zero (up to a gauge transformation) 
when restricted to T- 

Since T is complex Lagrangian in complex structure /, we can identify as follows the 
supersymmetry of the half-BPS brane produced by the Nahm pole. This is a brane of type 
{B,A,A), that is, it is a S-brane in complex structure /, but an j4-brane from the point 
of view of J or i^. 

8.2 General t 

At general t, we are dealing with a fiat bundle rather than a Higgs bundle. The Nahm pole 
still gives a line sub-bundle L G E, so we still have an exact sequence 



as in (8.1). The covariant derivative D/Dz now gives a holomorphic map E ^ E <^ K . 
We can still restrict this map to L and project the image to (g) K, to get a linear map 
V/Vz : L —7- (g) K. The condition that L is nowhere invariant under V/Vz implies, 
just as in section 8.1, that this map is an isomorphism from L to K, and again we 
conclude that L = K^l'^ . 




(8.4) 




(8.5) 
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The difference from section 8.1 is that now the bundle E is not a direct sum K^^"^ © 
but a non-trivial extension. Indeed, as we assume that the genus of C exceeds 
1, a bundle that holomorphically is a direct sum K^^"^ © K~'^/'^ would not admit a flat 
connection. 

Non-trivial extensions of K~^/'^ by K^^"^ are all isomorphic; this is so because such 
an extension is determined by an element of H^{C,K) = C, and the choice of a nonzero 
element does not matter, up to a bundle automorphism. 

The simplest example of a flat bundle that from a holomorphic point of view is the 
extension described in the last paragraph can be found by placing on C a Kahler metric of 
scalar curvature R = —1. Let u be the spin connection of such a metric and e the vierbein. 
The flat connection is 

A = L0t3 + e^t_+e;,t+. (8.6) 

In differential geometry, since co is the spin connection, the flat bundle E is the spin bundle 
of C, or more precisely the direct sum 

© /s:-i/2 of the two spin bundles of opposite 
chirality. But in this basis, the complex structure of E is defined by the (0, 1) part of A, 
which is Az = oJzH + ^zi- this is lower triangular, but not diagonal, so E is an extension 
rather than a direct sum. 

Having found a single fiat connection A on the bundle E, it is straightforward to find 
them all. We do not want to change A-z (since the holomorphic structure of E is supposed 
to be unchanged), but we can change Az by Az Az + ^z, where (to preserve fiatness) 
Xz is annihilated by D^. For E as described in the last paragraph, the relevant choice is 
= where ^ is a quadratic differential. 

Mathematically, a flat bundle E ^ C that from a holomorphic point of view fits in 
a non-split exact sequence (8.5) is called an oper; see [22] for a detailed explanation. We 
have learned that, at general t, the brane defined by the Nahm pole boundary condition 
is supported on the variety of opers. Actually, we should be more precise, because the 
Nahm pole boundary condition depends in general on a parameter ^ that was introduced 
in section 3.1, and the complex connection A that obeys the oper condition is in general 
not A-ir icj) but the more general connection A'^ defined in eqn. (3.3). We write for the 
subvariety of A4h defined by requiring that A"^ obeys the oper condition. 

When restricted to V^, A^ is fixed, up to a gauge transformation, so the complex 
symplectic form 

ni=^fdzdzTi5Ai5Ai (8.7) 

vanishes. Accordingly, the brane is a complex Lagrangian brane, just as in section 8.1, 
but now in a rotated complex structure. In the context of the present paper, the rotated 
complex structure is defined in section 3.1. might be called a brane of type {B, A, A)^^, 
being related to complex structure as a {B,A,A) brane is to complex structure /. 

As long as ( ^ 0, oo, the complex structures are all equivalent. If we simply set 
C = i, we get the usual variety of opers for complex structure J; alternatively, in the limit 
C — 7- 0, reduces to the holomorphic section V of the Hitchin fibration, described in 
section 8.1. 
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8.3 5'-Duality 

A particularly simple boundary condition in A/" = 4 super Yang-Mills theory is the Neumann 
boundary condition for gauge fields, extended to the whole supermultiplet in a half-BPS 
fashion. In terms of branes, this is the boundary condition for a family of D3-branes ending 
on a single NS5-brane in the absence of a gauge theory 0-angle. 

Upon compactification on C and reduction to two dimensions, this boundary condition 
gives a brane i^Nss OTi corresponding to a trivial fiat line bundle over Mh- In other 
words, the support of the brane ^ns5 is all of and its Chan-Paton connection is 

trivial. The brane Sns5 is of type {B, B, B), meaning that it is a S-brane in every complex 
structure. This reflects the fact that the trivial bundle over A4h is holomorphic in every 
complex structure. 

Under S-duality or electric-magnetic duality, the D3-NS5 boundary condition is con- 
verted to a D3-D5 boundary condition, still with 9 = 0. On the other hand, S-duality 
acts in the dimensionally reduced theory as T-duality on the fibers of the Hitchin fibration 
[60, 61]. Hence the brane 5ns5 must be mapped by 5-duality to a brane 5d5 supported on 
a section of the Hitchin fibration. Moreover, the S'-dual of a brane of type {B, B, B) is of 
type {B, A, A) (this is shown in [15]). For a middle-dimensional brane to be an A-brane, its 
support must be a Lagrangian submanifold, and its Chan-Paton bundle must be fiat. So 
the section of the Hitchin fibration on which Bd^ is supported must be complex Lagrangian 
from the point of view of complex structure /. On the other hand, concretely the D3-D5 
boundary condition, for the case of a single D5-brane, is described by the Nahm pole [19] 
at t = 1. (This value of t corresponds to unbroken supersymmetry of type {B,A,A).) 
Our analysis above determines the section of the Hitchin fibration that corresponds to the 
Nahm pole; it corresponds to the family of Higgs bundles described in (8.2) and (8.3). 

The D3-NS5 boundary condition can be deformed by turning on and more gener- 
ally by turning on a U{1) gauge field on the NS5-brane. These deformations, which are 
described in [19], preserve the half-BPS nature of the boundary condition but rotate the 
unbroken supersymmetry. A particular deformation in this family, described in section 12 
of [15], gives a brane - the canonical coisotropic brane Bcc ~ that is important in the gauge 
theory approach to the geometric Langlands correspondence. Bcc is a rank one brane sup- 
ported on all oi Mh, with a Chan-Paton bundle whose curvature is a linear combination 
of the Kahler forms oi Mh- The precise combination depends on a parameter analogous 
to our ^. With a choice that is convenient for geometric Langlands, the curvature of the 
Chan-Paton bundle is a multiple of uj (the Kahler form for complex structure J) and then 
Bcc is of type {A,B,A). 

The S'-dual of the deformation of ;Bns5 that gives Bcc is a deformation of Bd^ that is 
obtained by rotating the Nahm pole in the space of fields A and (p. This is analyzed in [19], 
and the appropriate type of "rotation" was briefly described in section 3.1. As we have 
seen, the rotated Nahm pole boundary condition leads to a brane Boper that is supported 
on the variety of opers, that is on for some If Sec is defined in the standard fashion 
as a brane of type {A,B,A), then its 5'-dual must have the same supersymmetry. In that 
case, BopcT is supported on the ordinary variety of opers, with Q = i. (In the present paper. 
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it is more natural for to be real.) 

That the S-dual of the brane Bcc is the brane i3opor supported on the variety of opers is 
important in mathematical treatments of the geometric Langlands correspondence. See for 
example [62] for an explanation of the role of opers in the geometric Langlands correspon- 
dence. The facts that we have just described give a gauge theory way to understand the 
5-duality between Bcc and .Soper- It has been argued [63] that the iS-duality between these 
two branes is important in understanding the AGT correspondence [64] as well as recent 
developments relating supersymmetric gauge theory and integrable systems [65]. The role 
of the S-duality between Bcc and Boper is more explicit in [66]. 

8.4 Monodromy Defects 

The concept of a Higgs bundle can be generalized by allowing singularities at isolated points 
Pi G C. For what follows, the case of interest will be a regular singularity. To introduce a 
regular singularity near a point p in C, we pick a local complex coordinate z near C and 
then we introduce polar coordinates r,6 with z = re*^. We select elements a, /3,7 in the 
Lie algebra t of a maximal torus T C G, and look for solutions of Hitchin's equations with 
a singularity at p of the form 

A = ad9 + ... 

dr 

(^ = /3 -fde + ... (8.8) 

r 

where the ellipses refer to additional terms that are less singular than 1/r. We call this sort 
of codimension two singularity a monodromy defect. The general theory of Hitchin's equa- 
tions adapts well to this situation [67] and one can define a moduli space AiH{C;p, a, /?, 7) 
of solutions that is a hyper-Kahler manifold with properties rather similar to what one has 
in the absence of the monodromy defect. Everything we will say generalizes in an obvious 
way to the case of any number of monodromy defects. 

Once we introduce monodromy defects, the limitation of some of the above statements 
to the case that the genus of C is at least 2 can be dropped. All above statements hold for C 
of any genus in the presence of a sufficient number of monodromy defects (for G = SU{2), 
the required number is 3 if C has genus 0, and is 1 if C has genus 1). 

In the context of = 4 super Yang-Mills on S x C, where S is another two-manifold, 
one can consider a monodromy defect supported on S x p, with p £ C. The singular 
solution (8.8) of Hitchin's equations embeds naturally as a solution of the four-dimensional 
equations (1.1). In the limit that G is small compared to S, the = 4 theory on S x C 
reduces to a sigma-model on S with target J^H{G;p,a, f3,^). In this description, there is 
an additional parameter rj that arises [68] as a theta-angle for the abelian subgroup of G 
that is unbroken along Exp. So quantum mechanically, a monodromy defect really has 
four parameters a, /3, 7, r]. Under S-duality, a monodromy defect of the above-described 
type in G gauge theory is mapped to a similar monodromy defect in G"^ gauge theory. The 
transformation of the parameters under 5-duality is (a, rj) — )• (r/, —a) while /3 and 7 are 
rescaled (for more detail see section 2.4 of [68]). 
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If we drop the subleading terms represented by the elhpses in (8.8), we find that the 
monodromy of the complex flat connection A = A + i(p is U = exp(— 27r(a — i'j))- The 
subleading terms do not modify the monodromy as long as U is regular - meaning that 
the subgroup of G that commutes with U has dimension equal to r, the rank of G. If U 
is not regular, there is an important subtlety, explained in detail in [68], section 3.3. For 
brevity, we will here consider only the case that G = SU{2), so that the nonregular values 
of U are only ±1. If C7 = ±1, the monodromy F of a connection of the form (8.8) is not 
necessarily conjugate to C/; on the contrary, generically it is in the "unipotent" conjugacy 
class containing the element 



A general element V of this conjugacy class is ±1 plus an arbitrary nilpotent matrix: 



The equation x"^ + yz = describes an Ai singularity C^/Z2. The singular point is located 
atx = y = z = where the monodromy V is precisely ±1; in other words, this is the 
case that the subleading terms in (8.8) do not correct the monodromy. In setting a and 
7 to special values at which U = ±1, we will assume that /3 remains generic. In this case, 
even if [/ = ±1, the solution has a "symmetry breaking direction" built in, given by the 
singular term in the connection proportional to /?. The effect of this is to blow up the Ai 
singularity, replacing C^/Z,2 with r*CP^. This important fact is established in [67]. 

The precise meaning of this T*C¥^ is that if U = ±1 and /3 = 0, then A4h{C; a, /3, 7) 
has a locus of Ai singularities, which parametrizes Higgs bundles for which the monodromy 
around p is precisely ±1. But if /3 / with U still equal to ±1, then this singular locus is 
blown up, replacing the singularities by a family of CP^'s. 

The moduli space A^//(C;p, q,/3,7) is invariant under shifting a by a cocharacter - 
for G = SU{2), this means that it is invariant under a — )• a + diag(i, —i), which is a shift 
that can be induced by a gauge transformation that has a singularity at p. However, we 
will be interested in brane constructions that are not invariant under such shifts of a, and 
for this reason, it will be best for our purposes not to view a as a periodic variable. 

The first brane that we want to consider is the oper brane. It is defined as usual by the 
Nahm pole boundary condition. For this, we take S to be M x M+, where is the usual 
half-line y > 0, and we impose the Nahm pole boundary condition at y = 0. Suppressing 
the M or time direction, the picture on C x R_|_ is sketched in fig. 26: there are monodromy 
defects supported on pi x R+, with respective parameters {oi, I3i,^i,r]i), and a Nahm pole 
boundary condition at y = 0. Of course, we need to explain what sort of singularity we 
want where the monodromy defect ends on a boundary with the Nahm pole. As usual, this 
kind of question is answered by finding a model solution with the desired singularity. For 
the present case, this has been done in section 3.6 of [14], and in greater generality in [69]. 

Let us set a — = Adiag(z,— i), for a complex parameter A, and consider a Higgs 
bundle E with a singularity of this type at, say, z = 0. If E (viewed in complex structure 




(8.9) 




x'^ + yz = 0. 



(8.10) 
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Figure 26. Monodromy defects in C x K_|_. supported on pi x IR+ where the pi are points in C. C 
is represented by the reetangle. We assume a Nahm pole boundary condition sX y ~ Q. 



Iq) is also an oper, then it can be described by the classical stress tensor 

* = -^ + .... (8.11) 

where we have omitted less singular terms. This formula is just like (3.48), with ja = ka/2 
replaced by A. Flat sections of E correspond to holomorphic solutions of the differential 
equation 

9^2 



— + t]i) = Q. (8.12) 



For generic A, one can find two linearly independent solutions with tpi = z ■^(l + X^i^iQ-z*), 
ip2 = z^~^^ (l + X^i^i Ci2;*). This means that, as expected, the monodromy is 

U = diag(exp(-27riA),exp(27ra)). (8.13) 

What happens if instead \ = k/2 with k G Z? There is always a solution ^2 = 
+ Cj-z*), but if we look for a solution with Vi = z~^/'^{l + ciz + . . .), we 
find that generically when we carry this expansion to order 2*^/2+1^ need logarithmic 
terms of order 2:'^'/2+i log z + . . . . The logarithmic terms are simply a multiple of (log z)iIj2- 
Accordingly, the monodromy around 2: = is actually generically of the unipotent form 

(8.14) 

for some complex constant s. 

So if we want the monodromy around z = to be trivial, we need to impose one 
condition on the subleading coefficients in the stress tensor (8.11), so as to get s = 0. This 
means that having trivial monodromy around z = is a middle-dimensional condition. 
Indeed, without this condition, a monodromy defect for G = SU{2) increases the complex 
dimension oi Mh by 2, but the trivial monodromy condition fixes 1 of the 2 parameters. 
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As this point is important, we will dwell on it a bit. Generically, the monodromy 
around the defect is an element of SL{2, C) (complex dimension 3) that obeys 1 constraint 
specifying its conjugacy class, leaving 2 complex parameters. For example, when a = 7 = 0, 
the conjugacy class is two-dimensional, as exhibited explicitly in (8.10). The condition of 
trivial monodromy (which is defined only when U = exp(— 27r(a — i'j)) equals ±1, and has 
no analog for other values) fixes 1 of the 2 parameters associated to the defect, so it leaves 
1 parameter. One can think of this 1 parameter as the direction of symmetry breaking 
associated to the term f5dr/r in eqn. (8.8). The choice of a symmetry-breaking direction 
determines a point in a copy of CP^; this CP^ is the projectivization of the fiber of E at the 
point p £ C where the monodromy defect lives. A more detailed explanation of the origin 
of this CP"*^ is as follows. First of all, because of the equation + yz = 0, the unipotent 
conjugacy class described in eqn. (8.10) is explicitly isomorphic as a complex manifold to 
C^/Z2, with an Ai singularity at x = y = z = 0. The singularity is precisely the point 
at which the group element V in (8.10) equals ±1. In the context of the construction of 
Mh as a hyper-Kahler manifold, the f3 parameter is a Kahler parameter that blows up the 
Ai singularity, replacing the conjugacy class C^/Z2 by its resolution, the Eguchi-Hansen 
manifold T*CP^. In the blowup, the singular point at the origin is replaced by a copy of 
CP^. See [67] for the interpretation of /3 as a blowup parameter, and [68] for a leisurely 
explanation of some of these matters. 

We have essentially already run into the fact that in this situation, vanishing mon- 
odromy is a middle-dimensional condition. Let us specialize to the case that C = CPI (we 
could similarly treat the case that C = C with an irregular singularity at infinity). We 
know from section 3.4 that for a given set of singular points Za and charges ka, a = 1, . . . ,d, 
there are finitely many opers with monodromy- free singularities. The condition that a fiat 
Gc bundle should be an oper is a middle-dimensional condition. To reduce to a finite set of 
opers with monodromy-free singularity, the condition of vanishing monodromy must also 
be middle-dimensional. (This assertion tacitly assumes that the two conditions intersect 
in a transverse fashion, which is in fact the case.) 

In fact, dropping the oper condition, we can explicitly describe the moduli space of 
solutions of Hitchin's equations on C, with monodromy defects characterized by = ka/2, 
for which the complex connection A has trivial monodromy around those points. As C 
is simply connected, a flat bundle on C with no monodromy around the points Pa is 
completely trivial as a flat bundle. The only possible moduli arise because the symmetry 
breaking associated to the parameters (3a (which we assume to be all nonzero) generates a 
copy of CP^ at each singular point pa- To get the moduli space, we must divide the product 
of these CP^'s by the automorphism group of the trivial flat bundle E; this is a copy of 
SL{2,C). So finally the locus U of solutions of Hitchin's equations corresponding to fiat 
bundles with trivial monodromy at each singular point is isomorphic to {CF^)'^ / SL{2,C). 
This is a complex submanifold of Mh in complex structure J (it is defined by a condition 
on the monodromies, which are holomorphic in that complex structure). Its dimension is 
d — 3, which is one-half the dimension of Ain- In fact, U is complex Lagrangian from the 
point of view of complex structure J; this is true roughly because each CP^ is complex 
Lag rangian in r*CP^. So the brane Btm supported on U with trivial Chan-Paton bundle 
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Figure 27. This figure differs from fig. 26 only in that each monodromy defect line ends on a 
singular monopole, indicated by a black dot on the right. Since the defect lines themselves are 
selfdual (with a suitable transformation of their parameters), the brane defined by ending of the 
defect lines is also selfdual. 

is a half-BPS brane of type {A, B,A). 

This gives us a new way to think about opers of trivial monodromy. They are intersec- 
tion points of two Lagrangian submanifolds of type {A,B,A) - one is the variety of opers 
and one parametrizes bundles with trivial monodromy. So the opers of trivial monodromy 
give a basis for the space of supersymmetric open strings stretching between the brane 
i3oper and the brane ^Striv We call this the space of {Boper, IStm) strings. Technically here 
we want the space of {Boper, Btm) strings in the S-model of type J. 

We can study this space of supersymmetric string states using S'-duality, which converts 
the i?-model of type J to the A- model of type ujk- S-duality converts the brane ;Bopcr to 
the canonical coisotropic brane Bcc, as we learned in section 8.3. It turns out that, as we 
describe shortly, Btm is mapped to itself by 5-duality (with the usual transformation of the 
monodromy defect parameters {oa, f3a,la,i]a))- So the S-dual of the space of {Boper,Btriv) 
strings is the space of {BccBtriv) strings, now viewed in the ^4- model of type ujk- The 
key aspect of this problem is that although the support of Btm is Lagrangian for ujk, 
it is actually symplectic for ujj - indeed, the support of fitriv is a complex submanifold 
in complex structure J, and accordingly has ujj as a Kahler form. This being the case, 
the problem of describing the space of (-Sec, -Striv) strings is governed by the analysis of 
quantization and branes in [70]. The space of (i3cc, i3triv) strings is obtained by quantizing 
the support U of Btriv] here U is viewed as a symplectic manifold with symplectic structure 

UJ. 

8.5 A Selfdual Brane 

There is a simple gauge theory explanation of why Btm is selfdual. Forgetting about su- 
persymmetry for a moment, we can think of a monodromy defect line as the Dirac string 
associated to a magnetic monopole that may have been improperly quantized. Hence a 
monodromy defect line can end on a singular magnetic monopole (fig. 27). Since mon- 
odromy defects are mapped to themselves by S'-duality (with some transformation of the 
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parameters), pictures in which the monodromy defects end on singular monopoles are sim- 
ilarly mapped to themselves by duality. 

Supersymmetry imposes some constraints on the values of the parameters at which 
such pictures exist. In the context of the i?-model of type J, the monodromy around a 
given defect line must be trivial if the defect line is going to end. This means that, in 
this -B-model, the picture of fig. 27 only exists if a = 7 = (here we will view a and r] 
as periodic variables). Of course, that is anyway the only case that the brane Btm can 
be defined. Dually, in the ^- model of type ujk, a picture like that of fig. 27 only exists 
if 7 = ?7 = 0. (For example, r] must vanish because the worldsheet theta-angle rj fails to 
preserve the topological supersymmetry of the A-model if the support of the monodromy 
defect ends at a place where the U{1) bundle along the monodromy defect is not trivialized.) 

In the context of the present paper, opers with trivial monodromy arise most directly 
from singular monopoles at y = 0. However, without changing anything essential, we can 
move the singular monopoles away from the boundary as long as we connect them to the 
boundary via monodromy defects, as in fig. 27. This has the advantage of making it 
obvious that opers with trivial monodromy are intersection points of two branes, and also 
making clear the selfduality of one of these branes. 

In the general context of a defect line ending on a singular monopole, the monopole 
may be incorrectly quantized. However, for A = k/2, which is equivalent to 7 = 0, a = 
(A;/2)diag(i, — i), the monopole at the end of the string obeys Dirac quantization, but the 
string is observable because we assume /3 7^ 0. 

8.6 Application To The Gaudin Model 

The selfduality of the brane Btriv provides a gauge theory explanation of the main result of 
[22, 25]: opers on CP^ with trivial monodromy correspond to simultaneous eigenvectors of 
the commuting Hamiltonians of the Gaudin model. Let us consider the duality between the 
space of (yScc^triv) strings and the space of (^Soper, -^triv) strings. The following discussion 
assumes familiarity with the framework of [70]. 

To construct the space of (^Scc, ^Btriv) strings, we have to quantize a moduli space 
(nf=iCP^)/5'-L(2,C), where CP^ is a copy of CP^ attached to the monodromy defect 
at z = Za- Quantization of CP^ gives an irreducible representation Ra of SU{2) of spin 
ja = ka/2, and quantization of (nf=i 'CPa)/5'L(2, C) gives a quantum Hilbert space that 
is the 5C/(2)-invariant part of ®aRa-, 

^ = {^aRaf""^^^ . (8.15) 

The classical commuting Hamiltonians of Hitchin's integrable systems can be interpreted 
(in the A-model of type uk) as {Bcc, -Bcc) strings. So they act on the space of (^cci B^m) 
strings. In fact, the Hitchin Hamiltonians become the commuting Hamiltonians of the 
Gaudin model. To demonstrate the last statement, one interprets the generators of the 
SU{2) action on Ra as arising from first order differential operators on CP^, whence the 
Gaudin Hamiltonians (3.23) become second order differential operators. The "symbols" 
(or coefficients of the leading terms) of these operators are functions on the base of the 
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Figure 28. Another selfdual brane can be constructed by replacing tlie ends of monodromy defects, 
whiclr we used in fig. 27, with junctions of monodromy defects, as depicted here. 

Hitchin fibration that are precisely the Hitchin Hamiltonians. So, reading this in reverse, 
the Gaudin Hamiltonians represent a quantization of the Hitchin Hamiltonians (and this 
quantization is unique, given the commutativity of the Hitchin Hamiltonians, modulo the 
possibility of adding c-numbers). 

To understand the eigenvectors and eigenvalues of the commuting Hamiltonians, we 
use the equivalence of (^cci -Striv) strings in the ^4- model of type ujk to (Bopcn Striv) strings 
in the i?-model of type J. The latter strings simply correspond to intersection points 
of the classical branes Bcc and .Soper- So opers with trivial monodromy give a basis for 
the quantum Hilbert space S) of the Gaudin model. In the i?-model description, the 
commuting Hamiltonians simply become functions on the variety V of opers, which is the 
support of the brane 6oper- Hitchin's classical Hamiltonians are holomorphic functions on 
the space of quadratic differentials on C (with poles of prescribed type at the positions Za 
of the monodromy defects). The support of V is the space of stress tensors on C (with 
prescribed poles at the Za)- The space of stress tensors differs from the space of quadratic 
differentials only because of the c-number conformal anomaly. This matches the additive 
c-number ambiguity in the quantization of the Hitchin Hamiltonians. The eigenvalues of 
the quantized Hamiltonians corresponding to a given oper are simply given by the stress 
tensor associated to that oper. 

Recently [66], a "noncompact" version of the Gaudin model has been described in 
which the finite-dimensional representations Ra are replaced by infinite-dimensional ones. 
The eigenvectors of the commuting Hamiltonians are again expressed as opers, now with 
certain conditions on their monodromies. It is natural to suspect that this construction 
again reflects the existence of a selfdual brane. There actually is a good candidate - a 
selfdual brane that is constructed by replacing ends of monodromy defects, as in fig. 27, 
by junctions of such defects, as in fig. 28. Such a junction is defined by a solution of 
Hitchin's equations on a small two-sphere S linking the junction with singularities (of a 
type depending on the parameters ota-, (iaila-ifla) at the intersection points of S with the 
monodromy defects. 
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A Three-Dimensional BPS Equations From Six Dimensions 

The time-independent configurations we consider in section 2 are solutions of the 3d BPS 
equations 

3 

[VuVj]=Q (A.l) 

i=l 

for a 3d connection together with three adjoint scalar fields, packaged together in the 
operators Vi as 

^ D . D 

^1 = + «- 



T^2 = [(t>2 - i(l)3, •] 

'^3 = §^-^[^l,■] (A.2) 



Or, in a complex notation. 



D 

=2— 
Dz 

V, = (A.3) 

This system of equations can be generalized to a one-real-parameter family of 3d BPS 
equations, which can be written as in A.l, but with a different choice of operators Pj: 

Pi = 2^ + 2C[^, •] 
Dz 

V2 = -2C^ + 2[ip, ■] 
V 

P3 = — . (A.4) 

The generalization was studied in section 3. 

This family of 3d equations can be usefully derived from six dimensions. We start in 
with coordinates x*^, x""'"'^, a = 1,2,3. Then we constrain a gauge field by requiring 
that the field strength, seen as an element of the SO{Q) Lie algebra, lies in a specified 
SU{?>) subgroup. As long as one is in six dimensions, the choice of a subgroup does not 
matter; it just amounts to the choice of an identification of with C^. But if we require 
that the fields actually only depend on the first three coordinates x", and are invariant 
under constant shifts of x""*"^, then the choice of an SU{3) subgroup does matter. So 
after dimensional reduction to three dimensions, one can obtain a family of inequivalent 
three-dimensional equations depending on a parameter. 
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A simple way to show that the family is of real dimension one, modulo equivalences, is 
as follows. First, a choice of embedding of SU{3) in SO{6), parametrized by SO{6)/U{3) ~ 
SU{4) /U{3) ~ CP^, is equivalent to the choice of a complex line in the space of 6d spinors of 
positive chirality. After dimensional reduction to 3d, the inequivalent sets of 3d equations 
are parametrized by such a choice modulo the 50(3) x SO{3) group of space rotations and 
rotations of the three scalars (fti. Although this group is six-dimensional, just like CP"^, it 
does not act freely on CP'^; rather, a generic point in CP^ preserves an SO{2) subgroup 
of SO{3) X 50(3). For example, all the 3d equations parametrized by ^ are invariant 
under a simultaneous phase rotation of D/Dz and (pz- In general, 50(3) x 50(3) acts as 
50(4) C 5C/(4) on the space of 6d spinors. So a complex spinor of 5f/(4) is a complex 
vector of 50(4), and its real and imaginary parts break 50(4) to 50(2). Hence the family 
of 3d equations obtained as dimensional reduction of the 6d equations is of real dimension 
one, and (A. 4) is a generic representative. 

A consequence of this picture is that we can change C by an 50(6) rotation. Indeed, 
we can change ^ as desired by acting with an appropriate element of a group that we will 
call 50(2)^;, which rotates D/Dx^ and [(/>2, •] into each other, and also rotates D/Dx^ 
and [(ps, ■] into each other. To be precise, the 50(2)^ rotation acting on the T>i defined 
at = will give a slightly rescaled version of the Vi, with a prefactor (1 + C^)""^/^. This 
prefactor can be absorbed by a simple rescaling of the z coordinate. 

Now let us discuss the Nahm pole boundary condition that has been so important in 
the present paper. If we assume a dependence on y only, and further assume that A-^ = 
(so that we can disregard T>i), the equations (A.l) with the T>i defined as in (A. 3) reduce 
to Nahm's equations. The Nahm pole boundary condition is defined by requiring that for 
y — 7- 0, the fields can be approximated by a certain singular solution of Nahm's equations. 

There is a similar boundary condition for the 3d BPS equations at generic Indeed, 
50(2)^; maps a solution of (A.l) which only depends on y to a solution of (A. 4) which only 
depends on y. More explicitly, taking the general form of the Pj in (A. 4), we can look 
for solutions which depend on y only, and such that Di reduces to 2d/dz, i.e. = —OP- 
Then ^2 = 2(1 + C,'^)[p, ■] and hence we can embed solutions of the Nahm equations as 
solutions of the general 3d BPS equations, at the price of a rescaling of the complex scalar 
93 by 1 + C^- This leads to the rotated Nahm pole boundary condition which we found 
useful in this paper. 

Of course, what we have just described is not the only embedding of the Nahm pole 
which would be compatible with the general 3d BPS equations. For example, at ^ 0, 
we could have chosen to look for an embedding in which P2 rather than Vi is trivial; this 
would lead to what we might call anti-opers - flat bundles with an oper-like constraint on 
their antiholomorphic structure, rather than on their holomorphic structure. Any rotation 
of our choice of Nahm pole by the ^7(2) subgroup of 5C/(3) which preserves P3 would 
produce a possible boundary condition, but we will generally stick to the "oper" Nahm 
pole. 

We will conclude with an alternative explanation of the meaning of the parameter C,. 
For flnite, non-zero C,, the can be rescaled and interpreted as a generic complex 3d 
connection. In Cartesian coordinates, we can denote the components of the connection as 
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T>a- The complex equations tell us that the connection is flat. Then we have a moment 
map constraint, which set to zero a certain constant linear combination of the commutators 

From this point of view, C only appears in the choice of moment map equation. A 
generic linear combination of the commutators is described by a 3 x 3 matrix of coefficients 

, ,ab 

j;^«^[P„,P,t] = 0. (A.5) 

a,b 

From eqn. (A. 4), we have a useful relation: oj^^ /oj^^ = C^. If = 1, uj"'' is symmetric, but 
in general that is not so. 

Generically, under linear coordinate redefinitions, there is a one-dimensional parameter 
space of possible w"''. For example, if the symmetric part of Lo"^ is positive definite, as it 
is for (A. 4), we can make it into the identity matrix 6"'^. Then the antisymmetric part i?"* 
can be rotated to live in the z,'z plane, and its magnitude is controlled by a single real 
parameter, which we can identify with 



B Three-Dimensional BPS Equations From Four And Eight Dimensions 

In this appendix, we will discuss how the 3d BPS equations of parameter can arise from 
time-independent solutions of the four-dimensional BPS equations (1.1). We will generalize 
the statement that the C = equations in three dimensions arise from the 4d equations at 
t = 1 if we drop the dependence on one coordinate, say x^, and also set Ai = cpy = 0. 

First we will show that this is not a feature of a specific choice. We can start with any 
choice of t, set d/dx^ = and set Ai and (jiy to any two linear combinations of the other 
three components of 4>, and the resulting 3d equations will be equivalent to the 3d BPS 
equations discussed in the last appendix for some value of the parameter ( in (A. 4). 

For that purpose, it is rather convenient to rewrite the 4d BPS equations in a compact 
form, as a dimensional reduction of BPS equations in eight-dimensional Yang-Mills theory. 
A succinct way to describe the desired eight-dimensional equations is to pick a Spin(8) 
spinor e of definite chirality and require 

FijT^-^e = 0. (B.l) 

If the curvature Fjj is understood as an element of the Lie algebra of SO{8), then the 
equations restrict the curvature to a Spin(7) subalgebra of SO{8). These are really 7 
equations, because of the obvious relation 

e^FjjF^'^e = 0, (B.2) 

as are antisymmetric. 

Dimensional reduction to four dimensions breaks SO (8) to a subgroup that we will 
call S0{4:)s X S0{4:)fi, acting respectively on the first four and last four coordinates. The 
spinor e decomposes into a piece ei which is left chiral under both S'0(4)s and 80(4) 
and a piece eR which is right chiral under both SO{4:)s and 5(9 (4) /j. If both and e/j are 
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non-zero, they fix a choice of a twisted <S'0(4)'^ diagonally embedded in 50(4)^ x 50(4)/?, 
such that eL and are 50(4)^ scalars. Then the 7 equations decompose under S'0(4)'j, 
into a triplet of self-dual two-forms, a triplet of anti-self-dual forms and a scalar equation. 
This is the form familiar from (1.1). 

We write a' as an abbreviation for a + 4 and adopt a complex notation with a as an 
abbreviation for a + ia' and a as an abbreviation for a — ia' . In order to bring the 8d 
equations explicitly to the form (1.1), it is useful to combine the T matrices to raising 
operators 

7a = Ta + iTa+4: (B.3) 

and lowering operators 

Ja = ra- iVa+A (B.4) 

with a = 1, ... ,4. We write \ for a state annihilated by the lowering operators, and \13) 
for its complex conjugate, a state annihilated by the raising operators. Being invariant 
under SO (4)'^, e is a linear combination of and \U)] being real, it is actually e = 
e~*°|il) + e*"|l5), for some real a. 

Then the 8d equations can be written in terms of the (2, 0), (1, 1) and (0, 2) components 
of the curvature Fah, F^i and F-p 

e-'^F^ + e'^'^eab'''F--^ = 

Y.^s.a = 0. (B.5) 

a 

When we reduce to 4d, the first equation tells us that the selfdual part of Re (e~*"Fafe) 
vanishes, as does the anti-selfdual part of Im (e"*"^^^). With (j) = dx"", we recover 

the familiar 4d equations 

(F-0A0 + tdA0)+ =0 
(F - (j) A (p-t-^d A<l))~ =0 

d^*<^ = 0, (B.6) 

with t = tan a. 

If we start from the 8d form of the equations, it is clear that solutions which are inde- 
pendent of some of the eight directions preserve additional supersymmetry. For example, 
any solution such that Fjs = for some / also satisfies 

F/jF^r^e = (B.7) 

and hence preserves the supersymmetry generated by the real anti-chiral spinor T^e of 
SO{8). The 7 equations remain independent, and describe a reduction of SO{7) to G2 
preserving a 7d spinor ey. 

Solutions that satisfy Fjg = and Fjy = preserve generically four spinors: e, r'''e, 
r^e, r''^e. The 7 equations then describe the reduction of 5*0(6) to SU{3) preserving the 
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super symmetries generated by a 6d complex spinor ee and its complex conjugate. They 
decompose into 3 complex equations and a real moment map condition 

[V,,V,] = j;[A,P|]=0, (B.8) 

i 

as discussed in Appendix A. 

This is exactly the situation we are in whenever in the four-dimensional equations 
(B.6), for any value of t, we set d/dx^ = and set Ai, (/)y to any two linear combinations of 
the remaining three scalar fields (j) in 0- Any such choices will produce a 3d reduction of the 
6d BPS equations, and hence, according to the analysis in Appendix A, will be equivalent 
to the standard 3d BPS equations for some C- The 3d BPS equations admit the oper-Nahm 
pole boundary condition. This will induce a boundary condition in the original 4d BPS 
equations, which will be some deformation of the standard Nahm pole boundary condition. 
Vice-versa, with this boundary condition, the usual vanishing theorems will guarantee that 
time- independent solutions arise from solutions of the corresponding 3d BPS equations. 

Finally, we will describe a simple explicit choice of reduction from 4d to 3d which gives 
whatever ^ we wish. Starting from q = and the standard reduction with Ai = (j)y = 0, 
we make simultaneous 50(2) rotations in the (a, a + 4) planes for a = 1, 2, 3, i.e. rotations 
of D/Dx'^ and [(f)a, •] into each other by angles 9a- (We do not make such a rotation for 
= y, as this would not behave well when we introduce a boundary a.t y = 0.) The 
rotation multiplies the creation and destruction operators by phases g^*^"/^, and hence the 
vacuum by the phase e"*^"^"/^. Hence it shifts the angle a by X](j^a/2, and acts 
correspondingly on the t parameter. 

In order to preserve the SO {2) symmetry that rotates and x^, it is natural to keep 
02 = 03- Given how the rotation transforms T>2 and P3 in (A. 3), we will then have clearly 
= tan ^2- 

Concerning the relation between 9i and 82, there are two particularly natural choices. 
If we want to keep three-dimensional topological symmetry along the boundary, we should 
keep 01 = 62 = 03 = 9- A rotation by these angles will change t to tan(30/2 -|- vr/4), and 
set C to tan0. On the other hand, if we content ourselves with two-dimensional symmetry, 
we can keep t = 1, by setting 61 = —202 = —203 = —20. Again, C, will be tan0. With 
this second choice, we deform only the Nahm pole boundary condition, and not the four- 
dimensional equations. 

C On Boundary Conditions And A Special Solution Of The BPS Equa- 
tions 

Here we will describe the Nahm pole boundary condition for the 3d BPS equations with 
generic allowing for singular monopoles on the boundary, and describe explicitly the 
model solution for the case of just one singular monopole. We work throughout on ]R2 X M 
(the generalization to C x ]R_|_ is straightforward). 

We will write the BPS equations simply as a flatness condition 

[A,Pj]=0 (C.l) 
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for a complex 3d connection T>i = di + [Ai, ■] together with a moment map constraint. 
Just as in eqn. (A. 4), the indices i = 1, 2, 3 refer to z, z, and y. The definition of the T>i 
differs from eqn. (A. 4) by a rescahng of T>2- 

This affects the relative normalization of the [Pi,!Dj^] terms in the moment map con- 
straint. Of course, we can always rescale the y coordinate with respect to z, z. If we write 
the moment map constraint as 

Y,uj'^[T>i,Vj^ = 0, (C.2) 
hi 

for a constant diagonal matrix w*-', the statement invariant under scaling is that = 
C^w^^. We will find it convenient to set uj^^ = C~^, <^^^ = 1, f^^^ = 1- If = then 

cj'^did] = dl + 29,^ (C.3) 

is the Laplace operator for a Euclidean metric on the half-space x that is normalized 
in a slightly unconventional way 

ds^ = dy2 + 2|dz|2. (C.4) 

This normalization will be useful later. 

The flatness condition (C.l) tells us that Pj = gdig~^ for a complex gauge trans- 
formation g, that is, a map from x to Gc. The moment map condition (C.2) is 
invariant under unitary (G-valued) gauge transformations g — )• Ug. We can eliminate the 
gauge-invariance by introducing the gauge- invariant hermitean matrix h = g^g. Then the 
moment map equation can be conjugated to 

uj'^di{h~^dih) = (C.5) 



or 



oj'^d^djh = J^{dih)h-^{d]h) (C.6) 



When = 1 , this equation says that the map h from the half-space x M+ to the quotient 
space G\Gc (endowed with its natural Gc-invariant metric Tr (/i~^d/i)^/2) is harmonic. 
Problems of this type are much-studied, but usually (for example, see [71]) in the context 
of a hyperbolic metric on X M , , rather than a Euclidean metric, as in our case. 

For simplicity, we will specialize to the case G = SU{2), so that G\Gc is a copy of 
hyperbolic threespace or AdSs. We can write 



9= . ..1/2 L . (C-7) 




for a real function Y and a complex function S. This is a general parametrization, in the 
sense that every g E SL{2, C) can be uniquely written in this form, modulo a unitary gauge 
transformation g ^ Ug. With this parametrization, we have 
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In these coordinates, the natural metric on takes a famihar form 



In general, in terms of the variables Y and S, the equations for h become 



OJ 



3 'J 



J^d^ (I'^^sjsj = 0. (c.io) 

In the framework of section 3.2, we want a boundary condition that is determined by 
the properties of the "small section." If we write s for the small section in the complex 
gauge Ai = 0, then in the unitary gauge with = gdig~^, the small section becomes gs. 
We must require gs to go as y^^^ as y — >• 0, while g itself diverges as This means 

that h diverges as while hs and s^h are finite and s^hs goes as y. The standard Nahm 

pole solution corresponds to s = f j and 




This formula, which is familiar from eqn. (3.8), is equivalent to 

h=( , V (C.12) 

\-zy \z\^y ' + yj 

Comparing to the general parametrization (C.8), we see that the standard Nahm pole 
solution is y = y, S = z. (The normalization uj^^ = C~'^i "^^^ = li = 1 was chosen 
to ensure that this is a solution for all C-) In other words, this solution is the "identity" 
map from the half-space x ]R_|. endowed with the Euclidean metric (C.4) to the half- 
space endowed with the hyperbolic metric (C.9). For = 1, the assertion that this gives 
a solution is simply the statement that the "identity" map between half-spaces endowed 
with these two metrics is harmonic. 

In general, if s = I q I ' want to require that y ~ y and P — QT, ~ y as y — ?■ 0. 

The last statement means that if we set cr{z) = P/Q, then S = o" at y = 0. The fact that 
y — 7- for y — 7- means that the boundary y = of the half-space x is mapped to 
the conformal boundary at infinity of the hyperbolic space H"^. That conformal boundary 
is a copy of CP^. By adjoing CP^ to H^, one makes the usual conformal compactification 
H of iJ^ The choice of an oper without monodromy determines a holomorphic map cr{z) 
from R2 ^ C to CP^ and the condition S| y=Q = a means that, as a map of the boundary of 
the half-space to CP^, h coincides with a. So our problem is this: given a holomorphic map 
a from the boundary of the half-space to the conformal boundary of the hyperbolic space, 
we want to extend u to a map h -.IS? x R_|_ — )• H that obeys (C.5) when restricted to y > 0. 
For = 1, we are simply trying to extend the given map cr to a harmonic map from the 
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half-space x ]R_|_ to H . (Technically, we assume that the map a has only polynomial 
growth so that it extends to a holomorphic map from the one-point compactification of C 
to CP^, and we similarly require that h extends to a continuous map from the one-point 

3 

compactification of the half-space to H .) 

For any cr{z), at least away from the branch points of the map a - in other words, 
the zeroes of da/dz - it is not difficult to write a systematic expansion of Y and Z for 
y — )• 0, involving powers of y and powers of logy. The expansion roughly starts with 
Y = y\a'{z)\ + ■■■ and S = cr{z) + and the coefficients are rational functions in 
derivatives of cr{z), and of three undetermined real functions of z and z. The denominators 
of these rational functions are powers of (t'{z) and its complex conjugate. So, away from the 
zeroes of cr'(z), boundary condition behaves well, and cuts in half the degrees of freedom 
of a solution. The branch points are precisely the points with PQ' — QP' = - in other 
words, the points at which there are singular monopoles. 

We still need to show that it is possible for a solution to be smooth away from the 
boundary in the presence of branch points or in other words singular monopoles on the 
boundary. The basic problem is to find a model solution in the presence of just one singular 
monopole; we then ask for the behavior near every singular monopole to match the model 
solution. In order to describe a singular monopole of charge k, we consider the special 



of the invariances of the BPS equations. The equations (C.IO) are invariant under scale 
transformations y — ?• Ay, z — ?• Az with real A, and under rotations z — )• e^^z. They are also 
invariant under reflections z — >• z of M^, accompanied, if ("^ 7^ 1, by S — t- S. 

The boundary conditions y ~ 1/y for y — )■ 0, S|y=o = z^^^ /{k + 1) are invariant under 
all these symmetries, accompanied by obvious rescalings of Y and S (which correspond to 
SL{2,C) transformations of the hyperbolic space). We expect the solution of the moment 
map condition that obeys the boundary condition to be unique, so it must be invariant 
under all these symmetries. Hence we require Y to be of the form y\z\^e^^^'^ and S to be 
of the form z^~^^ e^^P^ v{p) / {k + 1), for real functions u and w of p = y/\z\. 

Then the equations for h turn into two unfortunately rather complicated-looking non- 
linear PDEs: 



v{p) {e ((4fc + 3)p2 - 8) - p") u'{p) - Ak{k + l)pev{p) 
+p [p" + + 4) C^) v{p)u"[p) - (p3 + (p2 + 4) p(-2^ v{p)u'{pf 

+P {p" + {p' + 4) C^) v"{p) + [-p' - + 8) C^) v\p) = 
(A: + 1)V {p'+{p' + ^)e)n"{p) 
+v'{p) (2 (p2 + (p2 + 4) ^2) yi^p)u'[p) - A{k + l)pev{p)) 
+{k + l)pu'{p) {{k + (^2 + 1) - ACMp?) + <k + Ife {v{pf - 1) 

+ + (p2 + 4) ^2) ^(p)2„/(^)2 + (^2 + (^2 + 4) ^2) ^,(^)2 ^ (^.13) 



These equations involve v and the first two derivatives of u and v, but not u itself. Indeed, 
a constant shift of is a symmetry of the equations, though not of the desired boundary 
conditions for y — t- 0. 




other words (t(z) 



= z 



,k+l 



/{k + 1). We also make use 



- 113 - 
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Figure 29. The numerical solutions as 5 is varied across the critical value, for fc = 1. 

As PDEs for v and the derivative u', these equations have a space of solutions which 
is locally three-dimensional. The requirement that the solution should be smooth as 2 — )■ 
poses two constraints. It turns out that at large p, u behaves as klogp, so that Y ~ y^"^"*^, 
while v{p) scales as p^^, so that S ~ 6z^~^^ for some constant 6. The solution admits for 
large p a convergent power series expansion in 1/p, which depends on 6. 

On the other hand, the boundary condition at /3 — )■ is more forgiving, and only im- 
poses a single further constraint on the solution, which basically reduces to the requirement 
that t) —7- 1 as /o —7- 0. It is not difficult to check numerically that 6 can be tuned so that 
the solution satisfies the constraint, and it is hopefully possible to prove this rigorously for 
any non-zero finite C- As 6 is tuned, given the behavior for large p imposed in the last 
paragraph, there are two possible behaviors for v{p) as p becomes small. If 5 is small, v{p) 
does not reach 1, and goes to zero as p — 0. If 5 is large, it crosses 1 at some finite p, and 
then blows up before reaching p = 0. The solution we are after corresponds to the critical 
value of 6 which separates these two behaviors. 
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